FuzzyLogic Advanced Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору KChernov Ну потому видимо, что гвозди тоже можно микроскопом забивать, и в лаборатории где нет ничего кроме микроскопов этот метод даже будет оптимальным ... Всего записей: 1920 | Зарегистр. 27-07-2002 | Отправлено: 19:50 20-12-2010

aleksiu Newbie Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Andrew10 Цитата: Если будут еще ссылки, тоже принимаются, буду выбирать наиболее подходящий вариант.   Не так давно писал подобную вещь для себя ссылка на библиотеку и пример libFortran.zip при желании могу дать сам код библиотеки   Всего записей: 8 | Зарегистр. 10-12-2006 | Отправлено: 15:33 21-12-2010

Andrew10 Advanced Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору akaGM Цитата: а чисто фортран -- это, я думаю, только для тебя как разработчика или для "конечно-элементных" пользователей тоже Как раз возникла необходимость передать программу людям, которые могут только задавать входные данные и считать разные варианты, в том числе и разные свойства среды, заданные с помощью формул. Они не в состоянии каждый раз перекомпилировать программу при расчете нового варианта.     Конечно необязательно, чтобы то, что мне хочется, было реализовано обязательно на Фортране, но расчетный модуль написан на нем и тогда требуется передать в основной модуль значение функции по заданным координатам точки расчетной области и конечному числу параметров, от которых функция зависит.   Хотелось бы еще, чтобы была реализована работа с комплексными переменными, а в тех ссылках, которые здесь уже появились, насколько я увидел, доступна только действительная арифметика.   Конечно это уже "капризы" , но все же.   aleksiu Большое спасибо! Я пока посмотрю на все доступные варианты и, при необходимости, попрошу Ваш код.     Всего записей: 780 | Зарегистр. 26-02-2005 | Отправлено: 15:58 21-12-2010

Bear87 Newbie Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору   Вечер добрый всем    Подскажите пожалуйста как выйти из сложившейся ситуации.      Нужно проинтегрировать большую и страшную функцию f(x,y,z,r,...) по переменной х от a=0 до b=1. Функция имеет особенности в этих точках.        Нашел соответствующую процедуру в библ IMSL - QDAGS(f,a,b,res,ErrAba,ErrRel). Запустил для f(x)=dlog(x)/dsqrt(x). Все работает отлично.      Но интегрировать функцию f(x,y,z,r,...) она отказывается. Я так понимаю, функция нужна функция от одной переменной.      Пробовал:    -указать ей переменную интегрирования (call QDAGS(f(x),a,b,res,ErrAba,ErrRel) )- не работает;    -создать функцию g(x)=f(x,y,z,r,...) и проинтегрировать ее, при этом остальные переменные в основной программе задал. Интегрирует но, но при этом ост переменные, хоть убейте, но передавать не хочет и принимает их как нули. В итоге значение интеграла не верное.      Прошу помощи Заранее благодарю Всего записей: 31 | Зарегистр. 11-01-2010 | Отправлено: 18:54 21-12-2010

plazmon Junior Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Код:         FUNCTION g(x)       IMPLICIT NONE       Real g, f, x, y, z, r         y = ....       z = .....       r = ..... ...............       g = f(x, y, z, r)       RETURN       END FUNCTION   И в подпрограмму QDAGS передаешь g Всего записей: 107 | Зарегистр. 29-05-2007 | Отправлено: 19:32 21-12-2010

Bear87 Newbie Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Цитата:   y = ....       z = .....       r = .....      Дело в том что так я задавать не могу - вся эта подпрограмма по подсчету интеграла должна быть в цикле, предположим, по y на 3000 значений.      Может быть есть какой-нибудь другой способ передать в эту функцию значение величин y,z,r... кроме как в скобках в заголовке функции? Тогда все, по идее, должно решиться...    Может можно каким то образом объявить некоторые переменные, так, чтобы они были доступны для всей программы,тогда их и не нужно будет передавать? Всего записей: 31 | Зарегистр. 11-01-2010 | Отправлено: 19:55 21-12-2010 | Исправлено: Bear87, 20:01 21-12-2010

terminat0r Silver Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Bear87 Цитата: Может можно каким то образом объявить некоторые переменные, так, чтобы они были доступны для всей программы,тогда их и не нужно будет передавать?   Ура! Мы в очередной раз наблюдаем процесс изобретения колеса! Конечно же есть такие способы и о них здесь уже десяток раз писали. И выше один уже неявно упомянут. ключевые слова common block, module. И это описывается во всех учебниках языка фортран   Добавлено: Andrew10 Цитата: Хотелось бы еще, чтобы была реализована работа с комплексными переменными Как допилите до состояния модуля/библиотеки, не забудьте поделится. В хозяйстве это будет штука полезная. aleksiu Цитата: Не так давно писал подобную вещь для себя   Если возможно, то тоже бы не отказался. Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 20:38 21-12-2010 | Исправлено: terminat0r, 20:44 21-12-2010

pir0texnik2 Junior Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Bear87, я тебя не знаю? вот, сделай модуль Код:   MODULE moi_chudo_perenennie   IMPLICIT NONE   real*8 a,b,c,d,e Integer i,j,klm COMMON a,b,c,d,e,f,i,j,klm   END MODULE moi_chudo_perenennie   в функции напиши f Код: USE moi_chudo_perenennie и тоже в вызывающей процедуре, там же инить все переменные модуля. Если я тебя знаю, то интегрируй ты своего Фенхеля, эээ, т.е.... Ханкеля QUADPACK-ом, ну его тот ИМСЛ, ради одного интеграла тянуть. Живет он тут http://people.scs.fsu.edu/~burkardt/f_src/quadpack/quadpack.f90. Там интересны ф-ции QAG и QAGI и может QAGP. Почти все обозначения те же, что и в ИМСЛ. Если не знаю, то не надо интегрировать ф-цию Ханкеля вообще.   Комплексных авто-интеграторов я че-то в интах не видел, приходится компл. ф-цию РУКАМИ разделять на реальную и мнимую части. Да, неудобно с одной стороны, но с другой имно это быстрее, чем если бы это считалось автоматом... Кстати, интересно можно ли распараллелить этот квадпак... руки до этого не доходят все. Всего записей: 173 | Зарегистр. 27-02-2008 | Отправлено: 22:16 21-12-2010

terminat0r Silver Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору pir0texnik2 COMMON по-моему лишний, хотя это дело вкуса, а вот save кажется рекомендуют ставить даже в модуле с комонблоками   module moi_chudo_peremennie implicit none save real*8 a,b,c,d,e integer i,j,klm end module moi_chudo_peremennie Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 22:58 21-12-2010 | Исправлено: terminat0r, 22:59 21-12-2010

pir0texnik2 Junior Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору terminat0r /intel.forprodocs/intel.fprocompilerdocs/main_for/lref_for/source_files/rfsave.htm Цитата: Description In Intel® Fortran, certain variables are given the SAVE attribute, or not, by default:   The following variables are saved by default:   COMMON variables Module variables .....   Всего записей: 173 | Зарегистр. 27-02-2008 | Отправлено: 23:54 21-12-2010

terminat0r Silver Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору pir0texnik2 Не Intel-ом единым жив человек в даном случае вы правы но если бы там было COMMON /BNAME/ a,b,c,d,e,f,i,j,klm   то стандарт говорит, что с одной стороны Цитата: Execution of a RETURN or END statement might cause data objects in a named common block to become undefined unless the common block has the SAVE attribute, but never causes data objects in blank common to become undefined   а с другой Цитата:  A variable, common block, or procedure pointer declared in the scoping unit of a main program, module, or submodule implicitly has the SAVE attribute, which may be confirmed by explicit specification.   По-этому лучше  всегда уточнить компилятору, если есть такая возможность. Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 03:20 22-12-2010 | Исправлено: terminat0r, 03:21 22-12-2010

pir0texnik2 Junior Member Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Но ведь он таки безымянный - значит уже save. Но в принципе да, с сейвом было бы наверно спокойней. PS А раз уж такой разговор пошел, значит ли то , что save такой безопасный, что все время при заходе в процедуру делается его копия, а по выходу из нее его копируют обратно? Просто если предположить, что данных много, а процедура вызывается часто тоже много не скажется ли безопасность на скорости? Всего записей: 173 | Зарегистр. 27-02-2008 | Отправлено: 07:33 22-12-2010

aleksiu Newbie Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Andrew10   Забыл написать какие функции реализованы asin,asinh,sin,sinh,acos,acosh,cos,cosh,atan,atanh,tan,tanh,exp,lg,ln,sqrt,lg2,abs,sgn,int,ceil,floor,round,rnd,deg,rad   Добавлено: По просьбам трудящихся реализована комлексных числел libFortran.zip ввод комплексного числа 1+2i или 1+2j функции работающие с комплексными числами sin,cos,tan,exp,lg,ln,sqrt,lg2,abs Всего записей: 8 | Зарегистр. 10-12-2006 | Отправлено: 09:11 22-12-2010

Страницы: 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 (ФОРТРАН)

Переход по форумам  Закладки   » Компьютеры Цифровое изображение Программы Microsoft Windows UNIX Другие OS Hardware В помощь сисадмину Прикладное программирование Игры Форумные игры   » Интернет Web-программирование В помощь вебмастеру Графика Хостинг Зацените-ка!   » Тематические Ikonboard v.2 Ikonboard v.3 Invision Board Другие форумы PHP-Nuke и другие порталы Мобила   » Общие Спорт Флейм Темы по региональным встречам Музыка и Кино Юмор Литература и Живопись   » Ru.Board Новости Помощь по форуму   » Андеграунд eBookz Варезник Раздачи и акции Мобильный варезник Андеграунд Игры и фильмы   » Специальные Тестирование