Код Октре у місті Фортран

Я новачок у наукових обчисленнях. Я шукаю Фортран (бажано f90) впровадження указу.

Моя проблема потребує Октре, який розділяє мій домен, доки в будь-якому полі не буде більше ніж декілька N частинок (або джерел, де я знаю значення щільності, яке можна включити в метод інтегрального рівняння для вирішення рівняння, як рівняння Лапласа).

Все, що я міг знайти, це реалізація C ++ . Хочу знати, чи вже існують бібліотеки Fortran, якими я можу користуватися.

Також будуть оцінені хороші рекомендації з паперу щодо впровадження Octree у Fortran, які можна використовувати так, що швидке застосування методів інтегрального рівняння на обчислювальному полі стане просто!

Я не маю досвіду роботи з Octrees, але всякий раз, коли є якась приємна бібліотека C ++, яку я хочу використовувати у Fortran, я просто пишу простий драйвер C --- як правило, кілька функцій C, які роблять саме те, що мені потрібно. Потім я телефоную їм з Fortran за допомогою модуля iso_c_binding. Це має велику перевагу в тому, що ви повторно використовуєте добре перевірену бібліотеку із спільнотою навколо неї.

Щоб дізнатися, як використовувати iso_c_binding, я записав тут трохи інформації:

http://fortran90.org/src/best-practices.html#interfacing-with-c

ця конкретна інформація призначена для виклику Fortran від C. Я можу побачити, що мені слід оновити її. Ось мій код, де я викликаю бібліотеку C ++ з Fortran, перше посилання - простий драйвер C, друге посилання - інтерфейс Fortran, а потім я просто називаю його, як і будь-яку іншу підпрограму Fortran:

https://github.com/certik/hfsolver/blob/master/src/qc_libint.cpp

https://github.com/certik/hfsolver/blob/master/src/qc.f90

Ви, здається, дуже налаштовані на використання Fortran. При ефективній реалізації Octrees є досить складною структурою даних і, як така, краще підходить для мов програмування, які мають більшу підтримку для цього, наприклад, C / C ++. Існує ряд дуже якісних реалізацій в C / C ++, які ви могли використовувати.