Решение систем линейных уравнений

@eugrita · Регистрация: 18.11.2009

Author24 — интернет-сервис помощи студентам

у меня в этом казалось простом вопросе возникло 2 проблемы
1)с выбором подключения к IMSL в разных оболочках
как мне известно, в старом фортране 77 или еще более была программа Linb
здесь она недоступна .Единственн вариант - через IMSL
в CVF (Cjnpac Fisual Fortran 6.5) она есть но не смог подключиться и нет хелпа в скачке
пришлось остановиться на MS Fortran Power Station (более древнем)
(подключился оператором USE MSIMSL)
там проблема была в том, что все программы решения СЛАУ IMSL ждут в качестве входных параметров
массивы той размерности решаемой системы. Поскольку динамических массивов нет
фокусы типа

Fortran

1
2
3

                   DIMENSION A(100,100),C(100),X(100)
                   Read *,N
                   CALL LSARG (N,A,N,c,1, X)

не проходят -LSARG выдает ошибку типа матрица сингулярна
Ладно и из этого выкрутился введением оператора parameter
parameter (N=3)

Fortran

1 2	DIMENSION A(N,N),C(N),X(N) CALL LSARG (N,A,N,c,1, X)

правда приходится перетранслировать при каждом изменении размерности
Вопросы
1)можно ли передать в программу 2-мерный массив постоянной размерности но скажем N единиц по каждому измерению , например с N:100 по строкам и столбцам ?
2)как все таки работать с IMSL в CVF
3) Фортран имеет единственный плюс (по сравнению с С++) именно благодаря мощной поддержке математических библиотек. Недостатки - в сложности построения графика.
Правильно я выбрал версию фортрана. Может в более поздних MS Intell Fortran под Visual Studio 2003 и др удобнее работать с графиком ?
Я скачал версии CNL для C Visual Studio 2003 и IMS для Visual Studio 2007
может такая технология лучше (зная C) правдп придется ставить тяжеловесную студию.

@Игорь С · 22.09.2010, 12:22

Вообще то пользоваться библиотекой IMSL - неудобно. Во-первых та проблема, которой Вы упомянули. Во-вторых эти подпрограммы в случае ошибки (у меня такое было) сообщают о ней и приводят к останову приложения. Поэтому я в своих разработках использую старые программы в исходных текстах (известен пакет СSS таких программ, описанный еще в советской книжке). Например, для решения СЛАУ общего вида можно использовать программу GELG (Исходник у меня где-то есть). Там есть программы и для ленточных матриц. Еще много всяких подпрограмм было в пакете Белорусского института математики, но у меня нет исходников. Был еще американский пакет NAG, но пользоваться им менее удобно, т.к. там много неавтономных процедур, т.е. одна программа ссылается на другую и.т.д. хотя этот пакет очень обширный.

@eugrita · 22.09.2010, 20:02 **[ТС]**

спасибо конечно за сочувствие. Вообще-то мне в ближайшее время нужно от фортрана
применение в численных методах задач механики (конечно-разностные, прогонка, конечных элементов, собственные частоты) т.е в 1 очередь
1) программа рещения СЛАУ общего вида
2)Обращение 3-диагональных матриц (или решение СЛАУ)
3)собственные частоты и формы (не обязательно все достаточно с мин собств ч-той)
Ну а все таки по-поводу моиих вопросов?
1)Оставаться ли мне на MS Power Station или лучше на Compaq Visual Fortr
2)по поводу передачи матриц в подпрограмму
3)Готов пользоваться любым доступным для скачки или с исходниками математич пакетом
Возможно ли качнуть Css?

@Игорь С · 23.09.2010, 11:08

1. GELG - решает СЛАУ с системой общего вида, GELB - с ленточной несимметричной матрицей, GELS - с симметричной, MINV - оращение матриц общего вида, MCHB - два действия с ленточной симметричной положительно определенной матрицей: приведение к треугольному виду, и обратный ход (Ведь Вам для МКЕ, кстати, я сам этим занимаюсь).
2. Обращение 3-диагонольных матриц в CSS по-моему нет, а решение таких СЛАУ - есть.
3. Собственные частоты - это, кажется, определение собственных чисел?

1. Пока Вам не "тесно" в MS FPS то оставайтесь в нем, переход к Compaq VF прост.
Собственно говоря - если Вы "свободно" используете элементы Фортрана 95, то, конечно, переход необходим.
2. В CSS все матрицы в подпрограммах допускают переменные размеры.
3. CSS попробую добыть (где-то затерял).

@Игорь С · 28.09.2010, 17:58

Уважаемый eugrita. Куда отправлять пакет программ? Я его отыскал.
Приношу извинения за задержку и несколько неправильное название пакета, он называется SSP.
Кстати, тексты IMSL тоже отыскал.

@Игорь С · 05.10.2010, 19:04

Пакет

@eugrita · 06.10.2010, 22:17 **[ТС]**

спасибо за пакет. Позже попробую

@AlexKazancev · 23.01.2016, 13:49

Подниму тему, инфу нашёл на сайте http://fortran-90.narod.ru/slu.html#LA10 (кто-то давал ссылку в другой теме, продублирую)

Добавлено через 38 минут
Ну и там же Библиотека Численного Анализа НИВЦ МГУ http://num-anal.srcc.msu.ru/lib_na/libnal.htm
З.Ы. Не сочтите за флуд, просто устал листать описания функций LINPACK и LAPACK. А тут две библиотеки на русском с примерами.

@bobah16 · 25.01.2016, 14:34

AlexKazancev, жаль нет кода для решения дучп

@Игорь С 4 / 4 / 0 Регистрация: 01.03.2010 Сообщений: 181
	22.09.2010, 12:22	2
	Вообще то пользоваться библиотекой IMSL - неудобно. Во-первых та проблема, которой Вы упомянули. Во-вторых эти подпрограммы в случае ошибки (у меня такое было) сообщают о ней и приводят к останову приложения. Поэтому я в своих разработках использую старые программы в исходных текстах (известен пакет СSS таких программ, описанный еще в советской книжке). Например, для решения СЛАУ общего вида можно использовать программу GELG (Исходник у меня где-то есть). Там есть программы и для ленточных матриц. Еще много всяких подпрограмм было в пакете Белорусского института математики, но у меня нет исходников. Был еще американский пакет NAG, но пользоваться им менее удобно, т.к. там много неавтономных процедур, т.е. одна программа ссылается на другую и.т.д. хотя этот пакет очень обширный. 1

@eugrita 5 / 6 / 4 Регистрация: 18.11.2009 Сообщений: 661
	22.09.2010, 20:02 [ТС]	3
	спасибо конечно за сочувствие. Вообще-то мне в ближайшее время нужно от фортрана применение в численных методах задач механики (конечно-разностные, прогонка, конечных элементов, собственные частоты) т.е в 1 очередь 1) программа рещения СЛАУ общего вида 2)Обращение 3-диагональных матриц (или решение СЛАУ) 3)собственные частоты и формы (не обязательно все достаточно с мин собств ч-той) Ну а все таки по-поводу моиих вопросов? 1)Оставаться ли мне на MS Power Station или лучше на Compaq Visual Fortr 2)по поводу передачи матриц в подпрограмму 3)Готов пользоваться любым доступным для скачки или с исходниками математич пакетом Возможно ли качнуть Css? 0

@Игорь С 4 / 4 / 0 Регистрация: 01.03.2010 Сообщений: 181
	23.09.2010, 11:08	4
	1. GELG - решает СЛАУ с системой общего вида, GELB - с ленточной несимметричной матрицей, GELS - с симметричной, MINV - оращение матриц общего вида, MCHB - два действия с ленточной симметричной положительно определенной матрицей: приведение к треугольному виду, и обратный ход (Ведь Вам для МКЕ, кстати, я сам этим занимаюсь). 2. Обращение 3-диагонольных матриц в CSS по-моему нет, а решение таких СЛАУ - есть. 3. Собственные частоты - это, кажется, определение собственных чисел? 1. Пока Вам не "тесно" в MS FPS то оставайтесь в нем, переход к Compaq VF прост. Собственно говоря - если Вы "свободно" используете элементы Фортрана 95, то, конечно, переход необходим. 2. В CSS все матрицы в подпрограммах допускают переменные размеры. 3. CSS попробую добыть (где-то затерял). 0

@Игорь С 4 / 4 / 0 Регистрация: 01.03.2010 Сообщений: 181
	28.09.2010, 17:58	5
	Уважаемый eugrita. Куда отправлять пакет программ? Я его отыскал. Приношу извинения за задержку и несколько неправильное название пакета, он называется SSP. Кстати, тексты IMSL тоже отыскал. 0

@eugrita 5 / 6 / 4 Регистрация: 18.11.2009 Сообщений: 661
	06.10.2010, 22:17 [ТС]	7
	спасибо за пакет. Позже попробую 0

@AlexKazancev 7 / 7 / 2 Регистрация: 09.09.2012 Сообщений: 73
	23.01.2016, 13:49	8
	Подниму тему, инфу нашёл на сайте http://fortran-90.narod.ru/slu.html#LA10 (кто-то давал ссылку в другой теме, продублирую) Добавлено через 38 минут Ну и там же Библиотека Численного Анализа НИВЦ МГУ http://num-anal.srcc.msu.ru/lib_na/libnal.htm З.Ы. Не сочтите за флуд, просто устал листать описания функций LINPACK и LAPACK. А тут две библиотеки на русском с примерами. 0

@bobah16 373 / 343 / 42 Регистрация: 14.07.2015 Сообщений: 2,890
	25.01.2016, 14:34	9
	AlexKazancev, жаль нет кода для решения дучп 0