Матричная и линейная алгебра, ПО в помощь исследователям Опции

[Игорь]

Для конструирования статистических алгоритмов может оказаться полезным пакет программ "Математические и инженерные компоненты ME.com". Как и ПО AtteStat, данный пакет бесплатен и доступен для загрузки с сайта AtteStat. Он также работает в Excel всех версий. В его отладке также принимали участие заинтересованные специалисты, благодаря чему точность матричных расчетов неплоха. Доступны полные исходные тексты.

Приведем пример. Дана плохо обусловленная матрица. Найти обратную. Перемножить с исходной. Должна получиться единичная матрица. Первый расчет сделан ME.com. Второй расчет сделан с помощью известной программы Mathematica. Можно сравнить.

Сообщение отредактировал Игорь - 15.08.2010 - 14:42

Прикрепленные файлы

 Mathematica.zip ( 4 килобайт ) Кол-во скачиваний: 399

 

[Игорь]

LAPACK основан на BLAS ( Basic Linear Algebra Subprograms) говорить что это стандарт дефакто "спорно" можно только шутя.

Вполне серьезно. Цитирую: "BLAS (англ. Basic Linear Algebra Subprograms - базовые подпрограммы линейной алгебры) - стандарт де-факто интерфейса программирования приложений..." http://ru.wikipedia.org/wiki/Basic_Linear_...bra_Subprograms "Де-факто" именно и означает, что применяют все, но стандартом это не является. Кстати, иногда стандартами объявляют IMSL, NAG и прочее. И это тоже неверно.

Про "математические функции" Excek еще забавнее Первый вопрос --- какой версии этого самого Excel (он меняется (с моей точки зрения) чаще чем курс партии) и где референсная реализация? Пока нет кода программы никаким стандартом такой "черный ящик" вообще являться не может по определению.

Так что надеюсь что Вы шутите.

Правильно надеетесь. Т.к. о математических способностях Excel можно только шутя... Ну зато масштабы... Кстати, и Excel часто объявляют стандартной программой, хотя не стандартизован не только сам Excel, но даже формат его файлов.

PS ну а библиотека ESSL всего навсего дополнение к BLAS для процессоров поверпц от IBM

Не могу согласиться с формулировкой, хотя понимаю, что Вы имеете в виду - а именно: BLAS входит в ESSL. А то, пользуясь такой логикой, можно объявить автомобиль дополнением к его колесам или, ладно, согласен на большее, к его мотору.

Предлагаю точнее формулировать свою точку зрения.

Есть просьба. Не имею информации, как соотносятся SSP и ESSL. Если имеете сведения по конкретно данному вопросу, прошу поделиться.

Сообщение отредактировал Игорь - 23.12.2011 - 09:04

[p2004r]

Есть просьба. Не имею информации, как соотносятся SSP и ESSL. Если имеете сведения по конкретно данному вопросу, прошу поделиться.

вряд ли окажу конструктивную помощь, я думаю что оригинальная библиотека представляет скорее археологический интерес

http://www.ebyte.it/library/codesnippets/I...tem360_SSP.html

конечно для старого софта ищут многие

http://permalink.gmane.org/gmane.comp.emul...ules390.vm/2713
http://forum.academ.org/lofiversion/index.php?t209547.html

[Игорь]

вряд ли окажу конструктивную помощь, я думаю что оригинальная библиотека представляет скорее археологический интерес

http://www.ebyte.it/library/codesnippets/I...tem360_SSP.html

Спасибо, это хорошая ссылка. У меня есть по SSP русский двухтомник (качественная ксерокопия). Но нигде нет определенной информации, является ли ESSL обновленной или переписанной SSP или, напротив, ESSL - совершенно новый продукт. Что ESSL включает в себя известные библиотеки BLAS и LAPACK - такая информация есть. С другой стороны, у IBM есть руководство по переходу с LAPACK на ESSL. Документального подтверждения ни одной точки зрения не обнаружено.

От SSP в свое время отказался по причине того, что большинство функций в имеющемся пакете (также и введенные вручную из упомянутого выше двухтомника) оказались нерабочими, хотя описания алгоритмов (включая формулы) - верные. Полагаю, порча программ была сделана намеренно и очень неквалифицированно, а уверения переводчиков о якобы имевшей место проверке спорны. Поэтому все необходимые для наших задач функции были написаны нами заново, позднее переписаны на Си. Добавление новых алгоритмов как своих, так и переведенных на Си программ Брандта и еще некоторых авторов послужило тому, что сейчас представляет математическая и статистическая библиотека, лежащая в основе AtteStat и ME.com.

Перевод отдельных программ на Visual Basic for Application для Excel, предпринятый нами в порядке эксперимента, успеха не имел ввиду крайней медлительности VBA и совершенной его непригодности для решения каких-либо вычислительных задач. Использование же Excel в качестве front-end для программ на Си оказалось удобным. С уходом из Windows это решение, однако, потеряло актуальность.

Сообщение отредактировал Игорь - 23.12.2011 - 12:27

[100$]

Перевод отдельных программ на Visual Basic for Application для Excel, предпринятый нами в порядке эксперимента, успеха не имел ввиду крайней медлительности VBA и совершенной его непригодности для решения каких-либо вычислительных задач. Использование же Excel в качестве front-end для программ на Си оказалось удобным. С уходом из Windows это решение, однако, потеряло актуальность.

Не желая навязывать собеседникам никому не нужные дискуссии, замечу, что мой опыт общения с VBA, похоже, чуть более успешный: выборку из 36 000 случайных чисел (больше мой Excel не выдал без объяснения причин) реализованным на VBA методом сортировки вставками отсортировал за незаметное для глаза время. Ту же выборку за мгновение проверил на нормальность опять же самодельными критериями Колмогорова и Шапиро-Уилка. Программа Statistica 6.0 работать с таким массивом наотрез отказалась.

И еще. Есть такой эконометричсекий пакет EViews (7.0) Я по нему тестирую такие серьезные в вычислительном отношении процедуры, как тесты временных рядов на стационарность, состоятельное оценивание матриц вариации-ковариации регрессоров методом Уайта и Ньюэя-Уэста, тест Йохансена на коинтеграцию временных рядов. Пока что результаты совпадают до 9-10-го знака после запятой. Не шучу.

Ну, а такие пустяки, как критерий Фридмана, etc., уж не обессудьте, SPSS'ом проверял. SVD разложения матриц - MathLab'ом. Тоже без жалоб.

Сообщение отредактировал 100$ - 23.12.2011 - 14:46