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

@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, жаль нет кода для решения дучп

Новые блоги и статьи Все статьи Все блоги /
Работа с объемным DOM в javascript Htext 04.04.2025 Сегодня прочитал статью тут о расходах памяти в JS, ее утечках и т. п. И вот что вспомнил из своей недавней практики. Может, кому пригодится. Хотя, в той статье об этом тоже есть. Дело в том, что я. . .	Оптимизация производительности Node.js с помощью кластеризации run.dev 04.04.2025 Масштабирование приложений для обработки тысяч и миллионов запросов — обыденная задача для многих команд. Node. js, благодаря своей асинхронной событийно-ориентированной архитектуре, стал популярной. . .	Управление зависимостями в Python с Poetry py-thonny 04.04.2025 Стандартный инструмент для установки пакетов в Python - pip - прекрасно справляется с базовыми сценариями: установил пакет командой pip install и используешь его. Но что произойдёт, когда разные. . .	Мониторинг с Prometheus в PHP Jason-Webb 04.04.2025 Prometheus выделяется среди других систем мониторинга своим подходом к сбору и хранению метрик. В отличие от New Relic, который использует агентный подход и отправляет данные во внешнее хранилище,. . .	Пакет Context в Golang: Управление потоками и ресурсами golander 04.04.2025 Работа с горутинами в Go часто напоминает управление непослушными детьми - они разбегаются кто куда, делают что хотят и не всегда завершаются вовремя. К счастью, в Go 1. 7 появился пакет context,. . .
Контейнеризация React приложений с Docker Reangularity 03.04.2025 Контейнеризация позволяет упаковать приложение со всеми его зависимостями в автономный контейнер, который можно запустить на любой платформе с установленным Docker. Это существенно упрощает процессы. . .	Свой попап в SwiftUI mobDevWorks 03.04.2025 SwiftUI, как декларативный фреймворк от Apple, предоставляет множество инструментов для создания пользовательских интерфейсов. В нашем распоряжении есть такие API как alerts, popovers, action sheets. . .	Антипаттерны микросервисной архитектуры ArchitectMsa 03.04.2025 Хорошо спроектированная микросервисная система может выдержать испытание временем, оставаясь гибкой, масштабируемой и устойчивой к большинству проблем. Такая архитектура обладает высоким уровнем. . .	std::mutex в C++: Советы и примеры использования bytestream 03.04.2025 std::mutex - это механизм взаимного исключения, который гарантирует, что критический участок кода выполняется только одним потоком в каждый момент времени. Это простое, но могущественное средство. . .	Не удержался от оценки концепции двигателя Стирлинга. Hrethgir 03.04.2025 Сколько не пытался - она выдавала правильные схемы, причём случайно рисовала горячие области в середине, холодные по краям, трубки с краёв в низ и магнит в соединяющей, но при этой выдавала описание. . .

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

@eugrita 5 / 6 / 4 Регистрация: 18.11.2009 Сообщений: 661
	22.09.2010, 20:02 [ТС]
	спасибо конечно за сочувствие. Вообще-то мне в ближайшее время нужно от фортрана применение в численных методах задач механики (конечно-разностные, прогонка, конечных элементов, собственные частоты) т.е в 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
	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
	Уважаемый eugrita. Куда отправлять пакет программ? Я его отыскал. Приношу извинения за задержку и несколько неправильное название пакета, он называется SSP. Кстати, тексты IMSL тоже отыскал. 0

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

@AlexKazancev 7 / 7 / 2 Регистрация: 09.09.2012 Сообщений: 73
	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. А тут две библиотеки на русском с примерами. 0

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