Перейти из форума на сайт.

НовостиФайловые архивы
ПоискАктивные темыТоп лист
ПравилаКто в on-line?
Вход Забыли пароль? Первый раз на этом сайте? Регистрация
Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)

Модерирует : ShIvADeSt

 Версия для печати • ПодписатьсяДобавить в закладки
Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330

Открыть новую тему     Написать ответ в эту тему

akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Обсуждаются все вопросы, связанные с программированием на ФОРТРАН, как общего так и конкретного характера.
Постарайтесь дать как можно больше информации о возникшей проблеме -- это в конце концов в ваших же интересах чтобы вам помогли...

прежде чем просить помощи в задании
платное решение задач

ресурсы этого топика
ссылка на подборку ресурсов, собранных посетителями этого форума
 


если вам вдруг не отвечают или ответ вас не устраивает
и вообще полезно прочитать всем спрашивающим
 
просьба к пишущим и отвечающим все большие листинги оформлять тегом more
и отключать графические смайлики при размещении фортран-кода

Всего записей: 24961 | Зарегистр. 06-12-2002 | Отправлено: 18:11 14-01-2007 | Исправлено: akaGM, 15:26 15-05-2024
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
На libgen залили интересную книгу по сабжу. Перезалил
Gregory Ruetsch, Massimiliano Fatica-CUDA Fortran for Scientists and Engineers_ Best Practices for Efficient CUDA Fortran Programming-Morgan Kaufmann (2013).pdf  
http://rghost.net/50463392

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:01 26-11-2013
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
 
молодец, спасибо,
тока в шапку кинь...

Всего записей: 24961 | Зарегистр. 06-12-2002 | Отправлено: 20:45 26-11-2013
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
закинул.
 
Ув ALL
есть у меня системка ОДУ - ничего сложного - десяток немного длинных уравнений, но очень жесткая .  Математике с BDF это конечно как семечки, а вот у меня опыта дальше RK как бы и нет - все всегда решалось в лоб - тем что под рукой было.   Если у кого есть experience и отзывы по ODEPACK, Sundials, ...? для жестких уравнений, - т.е. по LSODA, CVODE, поделитесь. Без якобиана мне не обойтись, кажется, и вот гложут меня сомнения - делать его мне аналитически или будет достаточно и на коленке - численно. А если аналитически то что делать, если система поменяется? Прикручивать авто-дифференцирование?
Попробовал с FGSL - была у меня заготовка для системы ОДУ  - стало ясно точно - RK4 не подходит. Остальные как бы stiff solvers оттуда - просто висят и что-то там считают пол часа (кажется шаг по времени уменьшают) - не подходит. Брать NAG? - не очень и хочется - сорцы возможно надо будет и отдать потом.
 
Кстати в c++сников есть отличная (но слишком ООП-шная) odeint-v2 которая кажется даже CUDA умеет. Вот только тянуть все эти boostы не хочется.
 
Вот такие мысли вслух...

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 15:28 29-11-2013
mechmath

Junior Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
"есть у меня системка ОДУ - ничего сложного - десяток немного длинных уравнений, но очень жесткая .  Математике с BDF это конечно как семечки, а вот у меня опыта дальше RK как бы и нет - все всегда решалось в лоб - тем что под рукой было.   Если у кого есть experience и отзывы по ODEPACK, Sundials, ..."
 
Как показала практика решения жестких систем ОДУ лучше всего подходит решения методом Розенброка. В частности 4 стадийный метод.  
 
Кроме того, есть спектральные методы решение жестких систем ОДУ.  
S. M. Cox and P. C. Matthews. Exponential Time Differencing for Stiff Systems. J.Comput. Phys., pp.430-455, 2002.
Последний метод очень эффективен, только если дифф. уравнения можно представить в виде комбинации линейной и нелинейной частей. Причем вся жесткость должна сидеть нелинейной части.
 
Есть пакет программ (свободный) для решения на Фортране и Матлабе. Если интересно, могу прислать.

Всего записей: 50 | Зарегистр. 21-03-2013 | Отправлено: 20:58 29-11-2013
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
mechmath
Извините, отвлекся на революцию...
Спасибо за инфо! Очень познавательно.

Цитата:
Есть пакет программ (свободный) для решения на Фортране и Матлабе. Если интересно, могу прислать.  

Да, конечно!  terminat0r at yandex.ru

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 16:10 02-12-2013
BagaBaga

Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Хм, а можно и мне на notprinciple [злая собака] на yahoo.com
?

Всего записей: 463 | Зарегистр. 14-11-2005 | Отправлено: 20:49 02-12-2013
nadgol

Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Я конечно не специалист, но насколько я знаю, LSODE, ODEPACK, и т.д. это все одни и теже методы BDF, только слегка измененным  интерфейсом. А вообще-то все начиналось с GEAR. В Nag  и IMSL те же методы.
 
Советую взять с Netlib что либо по-новее (я в основном использую ODEPACK) и посмотреть численный Якобиан, на моих задачах считает нормально.
 
Успехов

Всего записей: 453 | Зарегистр. 28-09-2002 | Отправлено: 12:21 03-12-2013
BagaBaga

Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Я тут подумал, но всё же - вдруг тутошние гуру могут подсказать "правильные" фортрановские либы для решения (системы) PDE (да еще на сфере)?

Всего записей: 463 | Зарегистр. 14-11-2005 | Отправлено: 08:29 04-12-2013
mechmath

Junior Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
BagaBaga

Цитата:
Я тут подумал, но всё же - вдруг тутошние гуру могут подсказать "правильные" фортрановские либы для решения (системы) PDE (да еще на сфере)?

 
Что значит на сфере? В сферической системе координат?

Всего записей: 50 | Зарегистр. 21-03-2013 | Отправлено: 11:48 04-12-2013
BagaBaga

Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Тут сложнее. Есть моделька (унаследованная из 80х). В ней используются три системы координат - географическая сферическая, геомагнитная сферическая и геомагнитная дипольная. Система координат - широта (в действительности - коширота), долгота, высота (ну и геомагнитная широта, геомагнитная долгота, высота над поверхностью). Так что да, система координат - сферическая. Периодически в ряде расчётов возникают численные проблемы в районе полюсов (в районе географического полюса и в районе геомагнитного полюса). Есть желание поэкспериментировать, но переезжать на сетку Аракавы пока как-то не готов (если менять сетку, придётся пол-модели переписывать).

Всего записей: 463 | Зарегистр. 14-11-2005 | Отправлено: 19:27 04-12-2013
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
BagaBaga
Дальше 2D в PDEшках на фортране я не ходил, но в вашем случае посмотрел бы, как решает вашу задачу тот же Comsol Multiphysics например и сравнил с вашей реализацией. У меня достаточно хорошее впечатление от этого пакета - успешно решал на нем задачу диффузии для образца со сложной геометрией и структурой.
 
Наверное вы в курсе, но все же на всякий случай  
http://www2.cisl.ucar.edu/resources/libraries/classic
 spherepack, fishpack90, mudpack итп
А как универсальный сольвер нахваливали недавно на конференции http://www.mcs.anl.gov/petsc
Ах да и в принципе на полюсах кажется у всех проблемы, если это вас утешит.

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 14:30 05-12-2013 | Исправлено: terminat0r, 14:42 05-12-2013
BagaBaga

Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r,
спасибо.
 
Как решают проблему у полюсов ...  большинство используют решатели на основе различных модификаций метода контрольного объема. У полюсов чаще всего (в реальных моделях) пользуют искусственную вязкость, либо просто не считают выше 75-80 градусов. Хотя есть, правда, и спектральные модельки. Вроде как на С-сетках Аракавы таких проблем нет, но и моделей верхней атмосферы, её использующих, мне не известно Хотя вроде как "набирает обороты" в моделях океана и атмосферы.

Всего записей: 463 | Зарегистр. 14-11-2005 | Отправлено: 19:34 05-12-2013
nandr

Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Тесть обратился за помощью.
Он использует DIGITAL Visual Fortran Professional v6.0. Под XP он работает глючно, а под Win7 вроде бы вообще не ставится.
Какую более современную версию Фортрана ему посоветовать на замену, чтобы все его старые программы в ней работали без больших переделок, и чтобы при этом ему не очень сложно было перейти на нее (учитывая, что ему 84 года)? Его программы посвящены статистической обработке метеорологической информации.

Всего записей: 258 | Зарегистр. 01-11-2007 | Отправлено: 21:45 06-12-2013
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
nandr
 
если он не использовал QuickWin и какие-то сторонные библиотеки, то подойдёт свободный gfortran... ну или любую версию ворованного Интел Фортрана...
 
да, а как он компилирует, командной строкой или сидел в среде?
кстати, наверное не Дигитал, а DEC или Compaq, судя по номеру версии...

Всего записей: 24961 | Зарегистр. 06-12-2002 | Отправлено: 01:59 07-12-2013
nandr

Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:
если он не использовал QuickWin и какие-то сторонные библиотеки
Нет, сторонние библиотеки он не использует.

Цитата:
то подойдёт свободный gfortran...  
Я так понял, что это чисто компилятор, без Visual-среды?

Цитата:
да, а как он компилирует, командной строкой или сидел в среде?
Да, он в среде сидит.

Цитата:
кстати, наверное не Дигитал, а DEC или Compaq, судя по номеру версии...
Так DEC - это и есть Digital Equipment Corp.
Я скопировал обозначение из ридми.

Всего записей: 258 | Зарегистр. 01-11-2007 | Отправлено: 13:38 07-12-2013 | Исправлено: nandr, 13:58 07-12-2013
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
nandr
 
вообще-то компаковский 6.5 под вынь7 ставится, как -- не знаю, спроси здесь:
http://forum.ru-board.com/topic.cgi?forum=33&topic=4836&start=0
в крайнем случае всегда есть VM,
ну а за ворованнным интелом&студией тебе в варезник...

Всего записей: 24961 | Зарегистр. 06-12-2002 | Отправлено: 14:02 07-12-2013
kerpendel

Newbie
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Здравствуйте. Ребят подскажите пожалуйста, как случайные числа расставить в порядке возрастания. Нужно для того чтобы посчитать площадь методом трапеций. У меня сгенерировано 1000 точек на промежутке [-180:170]  
Спасибо=)

Всего записей: 2 | Зарегистр. 02-12-2013 | Отправлено: 00:07 08-12-2013
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
kerpendel
 
на Интеловском Фортране и подобных есть встроенные п/п sortqq() и qsort()

Всего записей: 24961 | Зарегистр. 06-12-2002 | Отправлено: 00:30 08-12-2013
nandr

Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:
ну а за ворованнным интелом&студией тебе в варезник...
Скачал с Интеловского сайта триальную версию Intel Visual Fortran Studio XE 2013 SP1. Установил ее. В "Программах" появилась папка "Intel Parallel Studio XE 2013". Но ни внутри ее, ни за ее пределами я не вижу заветной строчки для запуска Фортрана. Как его вызывать?

Всего записей: 258 | Зарегистр. 01-11-2007 | Отправлено: 17:59 08-12-2013
kerpendel

Newbie
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
AkaGM, Спасибо)

Всего записей: 2 | Зарегистр. 02-12-2013 | Отправлено: 00:53 09-12-2013
Открыть новую тему     Написать ответ в эту тему

Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330

Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)


Реклама на форуме Ru.Board.

Powered by Ikonboard "v2.1.7b" © 2000 Ikonboard.com
Modified by Ru.B0ard
© Ru.B0ard 2000-2024

BitCoin: 1NGG1chHtUvrtEqjeerQCKDMUi6S6CG4iC

Рейтинг.ru