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

[Игорь]

Для конструирования статистических алгоритмов может оказаться полезным пакет программ "Математические и инженерные компоненты 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

[p2004r]

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

Похоже это вечная беда Попытался взять в свое время тексты программ для решения интегральных уравнений из Верланя-Сизикова (вариант фортран кода), так да --- тест проходит, а например автоподбор невязки как класс не функционирует. Причем взятый из другого источника листинг в электронном виде в таком же состоянии. После того как пришлось по межбибу выписывать и ждать сборник стандартных подпрограмм был несколько разочарован. Ошибкой случайной такой эффект не объяснишь.