Какой стиль программирования наиболее подходит для численных расчётов?

Можно подробнее

0

Функциональный (как слышится, так и пишется).

1

Сообщение от Shamil1 Функциональный (как слышится, так и пишется).

Полностью согласен. Лично меня зацепил язык APL, "для себя" пишу именно на нём, но не волшебник, только учусь.

0

Shamil1, а можете объяснить почему? Например, в каком-нибудь сравнении. Ну так, коротко.

0

1. Численные расчёты - это чистые функции.
2. ФП - лаконичен и код очень близок к математическому описанию

0

Сообщение от Shamil1 Численные расчёты - это чистые функции.

Такой аргумент подходит и для обычного процедурного программирования. В интеловской MKL и в прочих реализациях интерфейса https://ru.wikipedia.org/wiki/... ubprograms - чистые функции.

0

https://vimeo.com/299593868

0

Сообщение от VTsaregorodtsev Такой аргумент подходит и для обычного процедурного программирования.

В данном случае это аргумент в обратную сторону. То есть, есть ряд задач, для которых ФП не очень удобно использовать, но численные расчёты к ним не относятся (так как численные расчёты - это чистые функции).

0

В императивщине с памятью понятней.

Сообщение от Вадим Тукаев APL

Во, векторное программирование! Была мысль взять какой-нибудь задачник по числякам и попробовать порешать на J. Но прошла.

1

Сообщение от helter Во, векторное программирование! Была мысль взять какой-нибудь задачник по числякам и попробовать порешать на J. Но прошла.

А я купил все (или почти все) книги про APL на русском языке. В первую очередь полезен "Курс АПЛ" Гилмана и Роуза. Кстати, хочу его купить в оригинале, там более новое издание. Но это уже скорее для коллекции, а пока и на русском сойдёт. Просто берём и решаем все задачки. Если идти постепенно, то всё легко и просто. Благо есть GNU APL, бесплатная версия Dyalog APL и даже онлайн-интерпретатор tryapl.org Подозреваю, что Вы решили изучить J методом кавалерийского наскока.

0

Сообщение от Вадим Тукаев Подозреваю, что Вы решили изучить J методом кавалерийского наскока.

Не уверен, что понял. В любом случае, я у меня скорее была идея познакомиться, чем изучать. J надо знать хотя бы чтобы ворчать на numpy там, где он отличается от J.

Добавлено через 12 минут
Любопытно, как в J реализованы массивы. Как так вышло, что массивы неизменяемые, а работает эффективно? В принципе, в заметках о реализации написано, если не вру, что, дескать, при вычислении выражения память выделяется malloc-ом, а после вычисления освобождается, и никакой сборщик мусора не нужен. Но, вероятно, есть и трюки. Например, в книге сказано, что изменение одного элемента массива такого вида:
не сопровождается массивным копированием данных. Видимо, J понимает, когда массив можно повторно использовать.

Хитрая там логика под капотом или не хитрая - всё равно магия. В numpy, например, такой магии нет. В документации пишут, когда временные массивы создаются, когда - нет; когда возвращается view в массив, когда - новый массив. Поэтому, по моим ощущениям, с императивщиной в общем спокойней, чем с функциональщиной, а векторный язык удобней скалярного.

0

Сообщение от helter Как так вышло, что массивы неизменяемые, а работает эффективно?

Есть структура данных persistent vector. Обычно реализована на основе 32-дерева, где все данные хранятся в листьях. Все операции, включая копирование - О(1).

0

Не уверен, что конкретно J использует хитрые функциональные структуры данных. Тем более, за персистенцию надо платить. По идее, под капотом сишные массивы и циклы по ним. Буквально из сишных структур данных реализация использует только одну: массив, содержащий немного метаданных. Скорее, он просто следит, когда можно сделать присваивание.

0

Сообщение от helter Не уверен, что конкретно J использует хитрые функциональные структуры данных. Тем более, за персистенцию надо платить.

Я не знаю, что использует J. Цена за персистенцию не высока и окупается при использовании рекурсивных функций, даже если не использовать неприятные для массивов функции, типа удаление произвольного элемента. Для тяжёлых операций типа сортировки обычно встроенная (библиотечная) функция копирует вектор в обычный массив, сортирует и копирует обратно в вектор.

0