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

НовостиФайловые архивы
ПоискАктивные темыТоп лист
ПравилаКто в 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
и отключать графические смайлики при размещении фортран-кода

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

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

Код:
 
...
  OPEN(N,FILE=STRING(1:ISTE), status='old')
 

 
и вообще непонятно... N читается с консоли, причём тут тогда входной файл?
 
 
-----------
гы
попробуй ещё
  OPEN(N,FILE=STRING(1:ISTE), status='old', share='denynone')
или просто shared (есть ли такое дело в gfortran, не знаю)
может прокатит...

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 18:31 21-08-2007 | Исправлено: akaGM, 19:37 21-08-2007
FuzzyLogic



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

Код:
 
program Test
  open(unit=10, file='test.in', status='old')
  open(unit=11, file='test.in', status='old')
  print *, 'Success'
end program Test
 

... в общем не работает ни в каких комбинациях
Про такие вещи как share/shared gfortran похоже не знает и ведь прав собака, по стардарту то :\

Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:57 21-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
раз уж речь зашла...
 
чем так хорош gfortran, ну в сравнении с тем же g77/g95?
 
просто как-то все эти самодеятельные порты... начинает студент какой-нибудь, потом бросает...

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 21:11 21-08-2007
terminat0r



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

Цитата:
раз уж речь зашла...  
чем так хорош gfortran, ну в сравнении с тем же g77/g95?

g77 больше не поддерживается.  
g95- это ответвеление несогласных с gnu gfortran. Производит неоптимизированый код (это конечно мое неофициальное сравнение на моей программе) и уступает gfortran который идет в ноздрю с интелл по производительности (я конечно не говорю о ОпенМП -там ситация иная немного).
Кроме того мне так и не удалось откомпиллировать g95 с последним gnu компилятором а пользоваться i386 компилятором на х86_64 машине нехочется только сейчас заметил что на странице наконец появились билды для других платформ- почти год ждать пришлось однако!  
gfortran хорош для меня например тем что это gnu, что он идет в стандартном линуксе вместе с системой.
А у нас в институте нелинукс стоит только один -возле сканера

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 22:06 21-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
спас, прояснил...
только всё же наши линуксоиды сидят на интел-фортране, он вроде тоже как свободный для Л
а с другой стороны -- куда им деваться, если всё идёт из-под виндов от него же...

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 12:35 22-08-2007
terminat0r



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

Цитата:
только всё же наши линуксоиды сидят на интел-фортране, он вроде тоже как свободный для Л

не у всех машин процессоры -интеловские
Время от времени приходится считать на самых разных кластерах. При этом очень важно что программа компилируется почти любым фортран компилятором. По сути у нас f77 программа с небольшими f90 усовершенствованиями- (динамические массивы etc). Дальше не идем умышленно чтобы была совместимость с разными компиляторами.
gfortran среди них наиболее универсальный пока. Конечно он отстает в введении новых инструкций языка, но это не критично. Параллелизация у нас ручная- одна точка на графике- один узел на кластере , поэтому и OpenMP пока не критично. Конечно со временем и это поджимать будет, но на горизонте пока ничего нету.  
Кстати, вот линк на тестирование производительности
http://gcc.gnu.org/wiki/GFortranResults
http://www.polyhedron.co.uk/55696351213626

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 16:27 22-08-2007 | Исправлено: terminat0r, 16:44 22-08-2007
akaGM

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

Цитата:
При этом очень важно что программа компилируется почти любым фортран компилятором

более чем согласен...
у меня (мои) программы ещё с VAX VMS (preF77 VAX-extension) и ЕС (F-IV) бегают...
 

Цитата:
одна точка на графике- один узел на кластере

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


тьфу ты, блин... только что допёр...
 

Цитата:
и вообще непонятно... N читается с консоли, причём тут тогда входной файл?

 
terminat0r, FuzzyLogic
а я думаю, что это вы никто мне не отвечаете...

Код:
 
      OPEN(5,FILE='input.in')
      CALL FOPEN
...
      READ(5,'(I3,A)') N,STRING  
 

 
FuzzyLogic

Код:
 
 program Test
  open(unit=10, file='test.in', status='old')
  open(unit=11, file='test.in', status='old')
  print *, 'Success'
end program Test
 
 
поставь юнит побольше >100...
как нас учили,  до 10 -- внутренние нужды...

Цитата:
...с консоли...

зря стандартные потоки перенаправлять...
может в этом дело?
недаром в файле числа дескрипторов с 20+ начинаются...

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 16:53 22-08-2007 | Исправлено: akaGM, 18:16 22-08-2007
terminat0r



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

Цитата:
что ж вы такое считаете, если не секрет?

рассеивание частиц в tddft

Цитата:
поставь юнит побольше >100...  
как нас учили,  до 10 -- внутренние нужды...

да нет, все равно не пашет,  действительно, gfortran- образец стандарта (правда недоделаный )

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 18:45 22-08-2007 | Исправлено: terminat0r, 18:47 22-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
вообще, опять же, как нас учили, данные определяют поток кода (ООП подождёт)...
значит, что-то с вашими входными данными не того...
 

Код:
 
22 ../../det/sppHeE0025.wp OLD
23 ../../det/sppHeE0025.wp OLD
 

 
смысл?
 

Цитата:
tddft

эээ...
термодинамическое что? фурье-дискрет-трансформ? гы

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 19:00 22-08-2007 | Исправлено: akaGM, 19:16 22-08-2007
terminat0r



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

Цитата:
смысл?

универсальность
 
 

Цитата:
эээ...  
термодинамическое что? фурье-дискрет-трансформ? гы

tddft

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:06 22-08-2007
akaGM

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

Цитата:
универсальность

я имел ввиду смысл такого формата, вернее, одного и того же файла на открытие...
 

Цитата:
tddft

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

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 19:12 22-08-2007
terminat0r



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

Цитата:
я имел ввиду смысл такого формата, вернее, одного и того же файла на открытие...

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

Цитата:
мы Шредингера напишем, упростим, пренебрежём

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

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:21 22-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
ну вот для этих "иногда" и напишите не меняя структуры файлов:
 
if (filename1 .equ. filename2) then
  call exclusive_sub()
else
  call ordinary_sub()

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 19:30 22-08-2007 | Исправлено: akaGM, 20:46 22-08-2007
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
это и пробую сделать сейчас. Вы вывели меня все-таки на путь истинный
Спасибо!

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 19:33 22-08-2007 | Исправлено: terminat0r, 19:38 22-08-2007
akaGM

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

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 20:46 22-08-2007
djdtyfhu



Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
В Фортране формирую массив строк:
   CHARACTER*4 string(4,250)
 
Как надо определить тип для такого массива в Delphi, чтобы можно потом было с помощью DLL вызвать процедуру Фортрана и на выходе обратно в Delphi в массиве иметь данные сформированные в вызываемой внешней процедуре.
 
Пытался сделать так:

Код:
type
  string_type=array[1..16,1..250] of char;
type
  pstring_type=^string_type;

 
Не помогает
 
 
P.S.  
Подробнее...

Всего записей: 422 | Зарегистр. 28-12-2004 | Отправлено: 14:05 23-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
djdtyfhu
не очень понимаю зачем тебе фортран нужен?
делай всё в Дельфи...

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 14:31 23-08-2007
djdtyfhu



Full Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
 
Смысл в том, что есть бинарные файлы сделанные программой, написанной на Фортране, и набор функций для чтения этих файлов. Вот и мучаюсь. Все победил: и числа, и массивы как целочисленные, так и с вещественными данными, даже отдельные строки вроде могу передавать, а с массивом строк пока не выходит.

Всего записей: 422 | Зарегистр. 28-12-2004 | Отправлено: 14:50 23-08-2007
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
djdtyfhu
а вот со строками может и не получится...
если будешь упорствовать, то совет:
переходи на "нуль-терминатед стрингс" , на С-строки то есть...
 

Цитата:
набор функций для чтения этих файлов.


Код:
 
 subroutine CALL_READ_STRING(pa_string,N_File,Name_REC)
   pointer(pa_string,string)
   OPEN(8,FILE=N_File,ACCESS='DIRECT',status='old',RECL=4000)
   CALL FREADC(8,Name_REC,string)
end subroutine CALL_READ_STRING
 

если это -- одна из таких функций, то в чём проблема?
назначение этой функции -- прочитать из файла блок данных, единственный вызов
FREADC() (Лахей что ли?)
 
ну так и пиши:

Код:
 
type
 TMyType =  array[1..16, 1..250] of char; //здесь определяй тот тип, который нужен
// TMyType =  array[1..4, 1..250] of widechar; //или что там тебе надо?
const
 RECLEN = 4000; // вообще лучше sizeof(TMyType)
var
 f: file of TMyType;
 mybuf: TMyType;
begin
  assign(f, filename);
  reset(f);//, RECLEN
  read(f, mybuf); //читаем _одну_ запись размера RECLEN
  closefile(f);
end.
 

 
и всё, переписана одна фор-функция CALL_READ_STRING()

Всего записей: 24055 | Зарегистр. 06-12-2002 | Отправлено: 15:52 23-08-2007 | Исправлено: akaGM, 15:55 23-08-2007
djdtyfhu



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

Цитата:
CALL FREADC(8,Name_REC,string)  

 
Это идет вызов как раз одной из тех функций, которые и обеспечивают сервис чтения бинарного файла данных. По вызываемому параметру Nаme_REC находится место в файле, с которого надо считывать данные в string. Хотя в итоге все равно все сводится к обычному  

Код:
Read(num_file,rec_num) string

 
Спасибо за наводку насчет самому написать процедуру чтения, попробую.

Всего записей: 422 | Зарегистр. 28-12-2004 | Отправлено: 16:11 23-08-2007
Открыть новую тему     Написать ответ в эту тему

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