Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё работало на VC6, в 7ой студии сплошные проблемы с этим фортраном — при импорте функций вылетают конфликты с libc, программа иногда неожиданно просто без сообщений об ошибках прекращает своё выполнение, вобщем бред какой-то... Может кто — нить подскажет где можно найти подобную процеджуру, желательно без необхоодимости включать в проект кучу файлов или огромную навороченную либу

Здравствуйте, Аноним, Вы писали:

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё работало на VC6, в 7ой студии сплошные проблемы с этим фортраном — при импорте функций вылетают конфликты с libc, программа иногда неожиданно просто без сообщений об ошибках прекращает своё выполнение, вобщем бред какой-то... Может кто — нить подскажет где можно найти подобную процеджуру, желательно без необхоодимости включать в проект кучу файлов или огромную навороченную либу

Посмотри в сторону LAPACK (Linear Algebra PACKage)
См http://www.netlib.org/lapack/
Есть порт пакета на С...

Мы как-то на него смотрели, но в итоге решили использовать IMSL.

Здравствуйте, bkat, Вы писали:

B>Здравствуйте, Аноним, Вы писали:


Есть библиотека http://www.gnu.org/software/gsl/gsl.html , единственна засада, нужно править сорцы под VC7


B>Посмотри в сторону LAPACK (Linear Algebra PACKage)
B>См http://www.netlib.org/lapack/
B>Есть порт пакета на С...
B>Мы как-то на него смотрели, но в итоге решили использовать IMSL.


А чем netlib не устроил? LAPACK и BLAS у меня прекрасно lib файлы собрались. Правда, работать с ним еще не пришлось.

Здравствуйте, peterbes, Вы писали:

P>Здравствуйте, bkat, Вы писали:

B>>Здравствуйте, Аноним, Вы писали:


P>Есть библиотека http://www.gnu.org/software/gsl/gsl.html , единственна засада, нужно править сорцы под VC7


B>>Посмотри в сторону LAPACK (Linear Algebra PACKage)
B>>См http://www.netlib.org/lapack/
B>>Есть порт пакета на С...
B>>Мы как-то на него смотрели, но в итоге решили использовать IMSL.


P>А чем netlib не устроил? LAPACK и BLAS у меня прекрасно lib файлы собрались. Правда, работать с ним еще не пришлось.

Нам нужна была одна функция, которой в LAPACK нету, а в IMSL есть.

Здравствуйте, <Аноним>, Вы писали:

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё работало на VC6, в 7ой студии сплошные проблемы с этим фортраном — при импорте функций вылетают конфликты с libc, программа иногда неожиданно просто без сообщений об ошибках прекращает своё выполнение, вобщем бред какой-то... Может кто — нить подскажет где можно найти подобную процеджуру, желательно без необхоодимости включать в проект кучу файлов или огромную навороченную либу


попробуй посмотреть в сторону blitz. это математическая библиотека, вот только не уверен насчет систем уравнений

Когда-то я работал над FEM процессором на C++, так для ускорения разработки использовал SuperLU. http://crd.lbl.gov/~xiaoye/SuperLU/ Библиотека на голом C плюс часть портирована с фортрана (f2c транслятором), поэтому при включении в проект погемороился неслабо. Но в итоге получил сносную либку (VC 6/7). Там есть и версия для распределённых вычислений.

Считает библиотека очень шустро, особенно если проникнуться тем что у них все матрицы столбиками задаются а не как привычно. Но там реально чистый С. C++-ом и не пахло.

Я юзал 2.0, на сайте уже 3.0 лежит, наверное ещё круче

Здравствуйте, Аноним, Вы писали:

А>Когда-то я работал над FEM процессором на C++, так для ускорения разработки использовал SuperLU. http://crd.lbl.gov/~xiaoye/SuperLU/ Библиотека на голом C плюс часть портирована с фортрана (f2c транслятором), поэтому при включении в проект погемороился неслабо. Но в итоге получил сносную либку (VC 6/7). Там есть и версия для распределённых вычислений.

А>Считает библиотека очень шустро, особенно если проникнуться тем что у них все матрицы столбиками задаются а не как привычно. Но там реально чистый С. C++-ом и не пахло.

А>Я юзал 2.0, на сайте уже 3.0 лежит, наверное ещё круче

Тоже для FEM надо Не подкинешь либку ? сюды ...

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++.

Boost.uBLAS не подойдёт?

Вот из boost-овской рассылки:

> I know of no such solvers in or out of uBLAS so I am not of much help.
> However, you may want to consider interfacing uBLAS with libraries which
> are not compatable with uBLAS. You will find such interfaces in the
> BOOST-sandbox CVS. I may be undertaking a similar exercise between uBLAS
> and some sparse system solvers. Also, there exists documentation for
> these bindings, but I cannot find the link

bindigs docs are also in the sandbox, in libs/numeric/bindings.
Introduction: http://tinyurl.com/66wb3
Traits: http://tinyurl.com/4vy6g
ATLAS bindings: http://tinyurl.com/5ty8f
LAPACK bindings: http://tinyurl.com/443zv
UMFPACK bindings: http://tinyurl.com/4qfol

Здравствуйте, Tonal-, Вы писали:

А>>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++.

T>Boost.uBLAS не подойдёт?

T>Вот из boost-овской рассылки:
T>

>> I know of no such solvers in or out of uBLAS so I am not of much help.
>> However, you may want to consider interfacing uBLAS with libraries which
>> are not compatable with uBLAS. You will find such interfaces in the
>> BOOST-sandbox CVS. I may be undertaking a similar exercise between uBLAS
>> and some sparse system solvers. Also, there exists documentation for
>> these bindings, but I cannot find the link

T>bindigs docs are also in the sandbox, in libs/numeric/bindings.
T>Introduction: http://tinyurl.com/66wb3
T>Traits: http://tinyurl.com/4vy6g
T>ATLAS bindings: http://tinyurl.com/5ty8f
T>LAPACK bindings: http://tinyurl.com/443zv
T>UMFPACK bindings: http://tinyurl.com/4qfol

Ну не знаю я boost'а Надо что нить что легко заюзать и просто вклеить в проект ... Пока в 6-ой сутлии сделал длл-ку, с одной функцией, в которую статически влинкован вызов функции из IMSL фортрановской.. Так вроде нормально работает, хотя это ясное дело времянка, заглушка что бы дальше двигаться

Здравствуйте, Аноним, Вы писали:

не забивайте себе мозги хитроумными boost-ми ( ей, право, гемора будет больше ). Если нужен лапак и блас по полной программе, берите http://www.netlib.org/lapack/ (хотя критика продукта уже здесь звучала), если нужно что либо попроще берите sgl (ссылки выше). Если поискать что круче и сильней то идите на mathtool или mathnet ( название забыл, но если есть интерес, то вспомню и дам точную ссылку )

Здравствуйте, Аноним, Вы писали:

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё

так сделайте COM в Compaq Fortran и вызывайте из любого языка

Здравствуйте, Lepsik, Вы писали:

L>Здравствуйте, Аноним, Вы писали:

А>>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё

L>так сделайте COM в Compaq Fortran и вызывайте из любого языка

Ну спасибо за совет, у меня этого фортрана вообще нет слава богу, тока старые либы от него остались

Здравствуйте, Аноним, Вы писали:

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё работало на VC6, в 7ой студии сплошные проблемы с этим фортраном — при импорте функций вылетают конфликты с libc, программа иногда неожиданно просто без сообщений об ошибках прекращает своё выполнение, вобщем бред какой-то... Может кто — нить подскажет где можно найти подобную процеджуру, желательно без необхоодимости включать в проект кучу файлов или огромную навороченную либу


MTL — Matrix Template Library. Доки и примеры в ИНет есть. Одна заковыка — сделана под STL. То есть, надо ещё и разобраться. Удачи

Здравствуйте, Аноним, Вы писали:

А>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений.
http://www.mathtools.net/

Здравствуйте, Аноним, Вы писали:

А>Здравствуйте, Аноним, Вы писали:

А>>Нужно найти срочно процедурку решения систем линейных алгебраических уравнений. Желательно также процедурку для решения разреженных СЛАУ на С++. Раньше использовал IMSL фортрановский, импортировала одну процедурку в С и всё работало на VC6, в 7ой студии сплошные проблемы с этим фортраном — при импорте функций вылетают конфликты с libc, программа иногда неожиданно просто без сообщений об ошибках прекращает своё выполнение, вобщем бред какой-то... Может кто — нить подскажет где можно найти подобную процеджуру, желательно без необхоодимости включать в проект кучу файлов или огромную навороченную либу


А>MTL — Matrix Template Library. Доки и примеры в ИНет есть. Одна заковыка — сделана под STL. То есть, надо ещё и разобраться. Удачи

http://www.osl.iu.edu/research/mtl/intro.php3
Системы она решать вроде не умеет...

А>MTL — Matrix Template Library. Доки и примеры в ИНет есть. Одна заковыка — сделана под STL. То есть, надо ещё и разобраться. Удачи

В свое время когда писал дисер с использованием FEM пришлось возится с разряженными матрицами, использовал сначала MTL, но потом перешел на uBLAS по скоростям выигрыш на моей проге был значителен. Да и mtl не развивается уже более 2-х лет.в uBLAS есть обычный LU солвер
смотреть эти функции
ublas::lu_factorize
ublas::lu_substitute
сам по себе uBLAS удобен для всевозможного преобразования матриц...а уже когда приходит время считать, то солвер можно прилепить через bindings. Удачи.