Вопросы с тегом «fortran»

В этом примере программы я делаю одно и то же (по крайней мере, мне так кажется) двумя разными способами. Я запускаю это на своем компьютере с Linux и наблюдаю за использованием памяти с помощью top. Используя gfortran, я обнаружил, что в первом случае (между «1» и «2») используемая память составляет …

13 performance  fortran 

Вывод здесь: Насколько лучше на самом деле компиляторы Фортрана? в том, что gfortran и gcc так же быстры для простого кода. Поэтому я хотел попробовать что-то более сложное. Я взял пример спектральной нормы. Сначала я рассчитываю 2D матрицу A (:, :), а затем вычисляю норму. (Это решение не допускается в …

13 fortran  c 

Тематические исследования в моей докторской диссертации требуют, чтобы у меня были расширения интервалов подпрограмм Фортрана в CHEMKIN-II (извинения за ссылку; это лучшее, что я мог найти для пакета, который больше не распространяется Национальной лабораторией Сандиа). CHEMKIN-II - химический пакет, используемый в химии горения; поскольку это стандартный пакет, используемый для сжигания, …

13 fortran  c  automatic-differentiation  interval-arithmetic 

Я просто заправлен в Фортран 95 для моделирования квантовой механики. Честно говоря, я был избалован Октавой, поэтому я воспринял возведение в матрицу как должное. Учитывая (малую, ) косоэрмитову матрицу размера n × n , каков наиболее эффективный способ использования LAPACK для решения этой проблемы? Я не использую оболочку LAPACK95, просто …

11 fortran  lapack 

Вот типичное использование указателей на функции в C. Я хотел бы сделать нечто подобное в Fortran. У меня есть некоторые идеи, но я хотел бы знать, есть ли какой-то канонический способ сделать это. Указатели и контексты функций, переданные пользователем, сохраняются, а затем вызываются позже. typedef PetscErrorCode (*TSIFunction)(TS,PetscReal,Vec,Vec,Vec,void*); PetscErrorCode TSSetIFunction(TS ts,Vec …

11 petsc  fortran 

Я новичок в научных вычислениях. Я ищу Fortran (предпочтительно f90) реализацию Octree. Моя проблема требует Octree, которая делит мою область, пока в любом блоке не будет больше, чем несколько N частиц (или источников, где я знаю значение плотности, которое можно включить в метод интегрального уравнения для решения уравнения, подобного уравнению …

10 fortran  integral-equations 

Недавно я задал вопрос в том же духе для косоэрмитовых матриц. Вдохновленный успехом этого вопроса, и после того, как пару часов ударился головой о стену, я смотрю на экспоненциальную матрицу вещественных асимметричных матриц. Путь к поиску собственных значений и собственных векторов кажется довольно запутанным, и я боюсь, что я заблудился. …

10 linear-algebra  fortran  lapack 

Недавно я столкнулся со странной проблемой с FORTRAN95. Я инициализировал переменные X и Y следующим образом: X=1.0 Y=0.1 Позже я складываю их вместе и печатаю результат: 1.10000000149012 После изучения переменных кажется, что 0,1 не представляется в двойной точности с полной точностью. Есть ли способ избежать этого?

10 fortran  floating-point 

Мне было интересно, как лучше всего научиться программировать на Фортране? Я искал книги, примеры, которые показались мне подходящими, были «Фортран 95/2003 Объяснено» и «Фортран для ученых и инженеров» . Я не знаю, что вы, ребята, думаете об этом? Еще одна вещь, которая меня интересует, это какой компилятор мне следует использовать? …

9 fortran 

Как видно из ответа Джеда Брауна « Затраты на поиск по сравнению с вычислениями» , использование векторизованных и не векторизованных операций с плавающей запятой приводит к гораздо более быстрому коду. Многие современные компиляторы утверждают, что они могут выполнять автоматическую векторизацию. Как мне определить, какие части моего кода успешно векторизованы?

9 fortran  c  compiling  vectorization