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

НовостиФайловые архивы
ПоискАктивные темыТоп лист
ПравилаКто в 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

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

akaGM

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

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

ресурсы этого топика
ссылка на подборку ресурсов, собранных посетителями этого форума
 
то, чем мы решили поделиться
ссылка на страничку программ etc собственного изготовления, которыми любезно делятся наши форумчане


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

Всего записей: 24053 | Зарегистр. 06-12-2002 | Отправлено: 18:11 14-01-2007 | Исправлено: akaGM, 09:47 01-03-2020
KChernov



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

Цитата:
Она по-моему нужна там, где используются старые функции (только не говорите что у вас все уже как минимум на f90), в которых нет вначале  
implicit none
и все константы double precision заданы именно в таком формате.  

Именно это я и скажу:
1. мне не нравится фиксированный формат фортрана (ну не перевариваю я его);
2. попытка использовать вместе фиксированный и свободный форматы у меня часто приводит к глюкам;
3. поэтому если мне приходится дорабатывать чужой код, первым делом я его преобразую в свободный формат (регэкспами) и привожу к модульному виду.
И почему-то после всего этого результат прекрасно работает
 

Цитата:
Конечно, разница несомненно в
Цитата:
Потому что у меня даже массивы так лежат чтобы кэш максимально эффективно использовался. Именно потому я за этот год 8 статей написал, и ещё штук 20 в соавторстве,
 
, но вот и мои примеры.
Пример с жизни - решение ур. Шредингера для гелия в сильном лазерном поле. Одна група решает в лоб- на суперкомпьютере в Англии (кажется несколько десятков тысяч процессоров), другая в Берлине на скромном beowulf кластере с несколькими сотнями процессоров. Методы разные- результаты почти одинаковые.  

Ну так я про то и говорю, что некорректно такие примеры приводить для прямого сравнения эффективности технологий (собственно именно в таком ключе тот пример и был подан).

Всего записей: 2471 | Зарегистр. 20-04-2004 | Отправлено: 17:50 08-12-2009
terminat0r



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

Цитата:
мне не нравится фиксированный формат фортрана (ну не перевариваю я его);  

А кому он нравится, но мы говорим не об этом. Я имею ввиду

Код:
       
      SUBROUTINE QTSMX(F,N,M)
      DIMENSION F(N,M),C(4),F2(4)
      ...
      C(1)=5.D0/16.D0
      C(2)=15.D0/16.D0
      C(3)=-5.D0/16.D0
      C(4)=1.D0/16.D0
      ...
      RETURN
      END

 
где явно не заданы типы массивов или переменных. Такого кода - куча.
 

Цитата:
поэтому если мне приходится дорабатывать чужой код, первым делом я его преобразую в свободный формат (регэкспами) и привожу к модульному виду.  

Ну - это не преобразование в f90. Так наверное все делают.
Но это не освобождает вас от сохранения "D" в  том случае, который я подал выше.
А иначе - точности - гаплык.  
 
Есть конечно контраргумент- компилятор орет благим матом на это все. Но если бы я стал исправлять все что не нравится компилятору сейчас, но нравилось когда-то, то я бы уже поседел.
Я часто использую nag, imsl, hsl  и кучу других. Всех их обьединяет только одно- все сорцы у меня с каких-то 80 годов и пробовать исправлять их и переписывать- просто гиблое дело.

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 18:13 08-12-2009 | Исправлено: terminat0r, 18:19 08-12-2009
FuzzyLogic



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

Цитата:
 Не верю в разницу на порядок (даже от 20 до 40 раз).  

Ну тут я не сильно расстроюсь, это ж не религия чтоб в неё верить или нет
 

Цитата:
Ну так я про то и говорю, что некорректно такие примеры приводить для прямого сравнения эффективности технологий (собственно именно в таком ключе тот пример и был подан).

 
Вообще я пример привёл не для сравнения эффективности технологий, а эффективности подходов, почему я и сказал что люди похоже решают совершенно разные задачи и подход использованный при решении одной задачи может совсем не подходить для другой. Почему такая большая разница в моём случае? (методы решения кстати практически одинаковые) У нас уравнения нормализованы и выписаны хитро ещё на уровне дифуров и двойная точность у меня на 1% переменных и там где у них везде real*16 у меня real*4 => кол-во необходимой памяти и использование кэша => кол-во необходимых компьютеров для обсчёта => разные уровни затыков из за сетевых коммуникаций через MPI.  
 
Вернемся к разным задачам и игре в верю-не-верю. Так вот, если бы сетка была раз эдак в 30 меньше и помещалась на один физ. компьютер (в плане памяти) то мой код работал бы может раза в 3 быстрее (именнно этот случай вы и подразумеваете когда сравниваете эффективность кодов). На данной сетке он работает быстрее на порядок, а вот если сетку увеличить ещё раз в десять (это мы подберёмся к критической точке для моего кода и используемых нами конкретных ресурсов итд), то я думаю что разница в счёте будет уже на два порядка после чего эта разница видимо останется приблизительно неизменной. Это к тому что "универсальность" оно конечно хорошо но вот реалии они зачастую отличаются от теории.

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 18:51 08-12-2009 | Исправлено: FuzzyLogic, 18:54 08-12-2009
KChernov



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

Цитата:
А кому он нравится, но мы говорим не об этом

Все подобные вещи я гребу туда же и так же и преобразую.
 

Цитата:
А иначе - точности - гаплык.

Ав чем проблема после убирания D прописать точность на уровне компилятора?
У меня так пока работало.
 
FuzzyLogic

Цитата:
Так вот, если бы сетка была раз эдак в 30 меньше и помещалась на один физ. компьютер (в плане памяти) то мой код работал бы может раза в 3 быстрее (именнно этот случай вы и подразумеваете когда сравниваете эффективность кодов).  

И какая у вас и у них сетка?

Всего записей: 2471 | Зарегистр. 20-04-2004 | Отправлено: 19:05 08-12-2009
FuzzyLogic



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

Цитата:
 И какая у вас и у них сетка?

одна и та же, просто у меня порог, после которого кривая скорости счета в зависимости от размера сетки начинает отклоняться от линейной, гораздо дальше чем у них за счёт более эффективного использования памяти, да и сама кривая падает не столь круто.
 

Цитата:
Ав чем проблема после убирания D прописать точность на уровне компилятора?
У меня так пока работало.  

Именно в том что в большинстве задач эта точность нужна далеко не везде а только для нескольких процентов переменных. А память она не резиновая.

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:06 08-12-2009 | Исправлено: FuzzyLogic, 19:11 08-12-2009
MegoChelovek



Junior Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Здрасте, у меня такая как бы проблема.......на Фортране я программирую буквально пару дней, нам по учёбе задали задачу : В заданном одномерном массиве найти :наибольший отрицательный и наименьший элементы.
Я написал код, тока вот компилятор мне ошибку выдает....пытался гуглить, нету такого, ну или я плохо искал.....если кто знает в чём ошибка напишите.....буду благодарен )
 
Тут тока первая половина кода, вторую половину не стал писать ибо все равно ошибку выдает )
 
 
integer j, i, k
print*, ' vvedite razmer massiva '
read*, i
integer a (i)
print*, ' vvodite elementi massiva '
 
do j=1, i
read*, a(j)
end do
 
! находим первый отрицательный элемент массива
do j=1, i
if a(j)<0 then
k= a(j)
!end if
end do
 
! ищет набольший отрицательный элемент
 
do j=1, i-1
if ( a(j)<0 .and. a(j)>k ) then
k=a(j)
end if
end do
print*, ' набольший отрицательный элемент '
print*, k
 
end
 

Всего записей: 33 | Зарегистр. 08-12-2009 | Отправлено: 19:12 08-12-2009
terminat0r



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

Цитата:
Ав чем проблема после убирания D прописать точность на уровне компилятора?  

Потому-что компиляторы бывают очень разные и  полагаться в этом вопросе на компилятор-  неправильно.
А в стандарте было когда-то записано  

Цитата:
To cause the constant assigned to X to be treated as DOUBLE PRECISION, you must write "X = 0.0025D0", where the "D" exponent denotes a DOUBLE PRECISION constant.

значит так будет на ЛЮБОМ компиляторе.

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:13 08-12-2009
FuzzyLogic



Advanced Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
MegoChelovek
У вас перед end if стоит восклицательный знак, что превращает его в комментарий.
Вы б хоть написали какую ошибку выдаёт и каким компилятором пользуетесь. И ещё, когда постите код кладите его в тэг code или more плиз.
А ну и ещё скобки поставьте ...  if a(j)<0 then     должно быть  if (a(j)<0) then

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:18 08-12-2009 | Исправлено: FuzzyLogic, 19:19 08-12-2009
terminat0r



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

Код:
!end if  

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:18 08-12-2009
Vskazka

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

Цитата:
integer j, i, k
print*, ' vvedite razmer massiva '
read*, i
integer a (i)

 
надо

Цитата:
integer j, i, k
allocatable, integer:: a(:)
print*, ' vvedite razmer massiva '
read*, i
allocate (a (i))

 

Всего записей: 382 | Зарегистр. 24-11-2003 | Отправлено: 19:22 08-12-2009 | Исправлено: Vskazka, 19:22 08-12-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Vskazka
FuzzyLogic
После работы все бросились задачки решать

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:24 08-12-2009
MegoChelovek



Junior Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
proga3.f90(13): error FOR3852: syntax error detected between IF and A
Error executing fl32.exe.
 
вот ошибка...
 
а насчет ! end if , там вроде бы и не надо его прописывать, это я просто закоментил его
 
Ок спасибо, щас попробую исправить

Всего записей: 33 | Зарегистр. 08-12-2009 | Отправлено: 19:25 08-12-2009 | Исправлено: MegoChelovek, 19:26 08-12-2009
FuzzyLogic



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

Цитата:
 После работы все бросились задачки решать  

Хорошо вам, после работы ... а у меня 11.30 утра, перерыв на кофе/подумать, будем считать это за "подумать"

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:29 08-12-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
FuzzyLogic
Ага, зато вы за морями-океанами
В принципе у меня тоже рабочий день еще не закончился, но все не так трагично как у вас.
 

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:42 08-12-2009 | Исправлено: terminat0r, 19:43 08-12-2009
MegoChelovek



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

Код:
 
 integer j, i, k  
 !allocatable integer: a(  
 print*, ' vvedite razmer massiva '  
 read*, i  
 integer a (i)
 print*, ' vvodite elementi massiva '
 
 do j=1, i
  read*, a(j)
 end do
 
 ! находим первый отрицательный элемент массива
 do j=1, i
  if a(j)<0  then
   k= a(j)
  !end if
 end do
 
 ! ищет набольший отрицательный элемент
 
 do j=1, i-1
  if ( a(j)<0 .and. a(j)>k ) then
   k=a(j)
  end if
 end do
 print*, ' набольший отрицательный элемент '
 print*, k
 end
 

 
я тут маленько по другому массив объявил, ибо как вы сказали у нас в методичке такого нету
и всё равно копилятор выдает ошибку в 14 строчке ( if a(j)<0  then ) , ошибка: error FOR3852: syntax error detected between IF and A
 
что означает эта ошибка, и почему она вылазит....так то вроде всё правильно наверно.....
 

Всего записей: 33 | Зарегистр. 08-12-2009 | Отправлено: 19:52 08-12-2009
FuzzyLogic



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

Цитата:
и всё равно копилятор выдает ошибку в 14 строчке ( if a(j)<0  then )

Вы посмотрите где я вам написал скобки поставить и где вы поставили ...
 
Добавлено:

Цитата:
я тут маленько по другому массив объявил, ибо как вы сказали у нас в методичке такого нету

Так как вы объявили делать нельзя, можно так:

Код:
 
integer   N
const (N=1000)
integer a(N)
 

и засуньте это сразу после integer i,j,k  "объявления" надо писать до самого кода

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:54 08-12-2009 | Исправлено: FuzzyLogic, 19:57 08-12-2009
Vskazka

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

Цитата:
ибо как вы сказали у нас в методичке такого нету

Ну ежели нет, то прикинте, сколько по-максимум у Вас может быт членов массива а (например 1000) и тогда
Цитата:
integer j, i, k
 integer a(1000)
print*, ' vvedite razmer massiva '
read*, i  

В фортране объявления всего на свете должно быть до первых исполняемых команд. Не плюсы чай.
 

Всего записей: 382 | Зарегистр. 24-11-2003 | Отправлено: 19:59 08-12-2009
MegoChelovek



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

теперь другие ошибки вылазят в 5 строчке (  integer a (i)  ), ошибки :
1. error FOR2588: Type declaration statement cannot occur here
2. error FOR2942: symbol I is not permitted in a specification expression
 
как их править, подскажите ?

Всего записей: 33 | Зарегистр. 08-12-2009 | Отправлено: 20:02 08-12-2009
Vskazka

Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
MegoChelovek
Только что Вам написали два варианта начала программы... Смотрите

Всего записей: 382 | Зарегистр. 24-11-2003 | Отправлено: 20:04 08-12-2009 | Исправлено: Vskazka, 20:04 08-12-2009
FuzzyLogic



Advanced Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
MegoChelovek
Вы как-то читаете что вам пишут через строчку ...

Цитата:
 1. error FOR2588: Type declaration statement cannot occur here  

Это когда объявление переменных стоит после исполняемого кода

Цитата:
 2. error FOR2942: symbol I is not permitted in a specification expression  

нельзя объявлять массивы переменного размера, читайте выше

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 20:06 08-12-2009
Открыть новую тему     Написать ответ в эту тему

Страницы: 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

Компьютерный форум 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