Mã Octree ở Fortran

Tôi mới làm quen với máy tính khoa học. Tôi đang tìm kiếm một triển khai Fortran (tốt nhất là f90) của một Octree.

Vấn đề của tôi yêu cầu một Octree phân chia miền của tôi cho đến khi không có nhiều hơn một số hạt N (hoặc nguồn mà tôi biết giá trị mật độ có thể được cắm trong một phương trình tích phân để giải phương trình như phương trình Laplace) trong bất kỳ hộp nào.

Tất cả những gì tôi có thể tìm thấy là việc thực hiện C ++ này . Tôi muốn biết nếu đã có thư viện Fortran hiện có mà tôi có thể sử dụng.

Ngoài ra, các đề xuất giấy tốt về triển khai Octree trong Fortran có thể được sử dụng để áp dụng các phương pháp tích phân nhanh trên hộp tính toán trở nên dễ dàng, sẽ được đánh giá cao!

Tôi không có kinh nghiệm với Octrees, nhưng bất cứ khi nào có thư viện C ++ đẹp mà tôi muốn sử dụng trong Fortran, tôi chỉ cần viết trình điều khiển C đơn giản --- điển hình là một vài hàm C thực hiện chính xác những gì tôi cần. Sau đó, tôi gọi họ từ Fortran bằng mô-đun iso_c_binding. Điều này có lợi thế lớn là bạn sử dụng lại một thư viện được thử nghiệm tốt với cộng đồng xung quanh nó.

Để biết thông tin về cách sử dụng iso_c_binding, tôi đã viết một số thông tin ở đây:

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

thông tin cụ thể này là để gọi Fortran từ C. Tôi có thể thấy rằng tôi nên cập nhật nó. Đây là mã của tôi nơi tôi gọi thư viện C ++ từ Fortran, liên kết đầu tiên là trình điều khiển C đơn giản, liên kết thứ hai là giao diện Fortran và sau đó tôi chỉ gọi nó như bất kỳ chương trình con Fortran nào khác:

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

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

Bạn dường như rất sẵn sàng sử dụng Fortran. Octrees, khi được triển khai hiệu quả, là các cấu trúc dữ liệu khá phức tạp và, do đó, phù hợp hơn với các ngôn ngữ lập trình có nhiều hỗ trợ hơn cho việc này, chẳng hạn như C / C ++. Có một số triển khai chất lượng rất cao trong C / C ++ mà bạn có thể sử dụng.