Finite Difference Frequency Domain

[la4361]

Hello everybody!
I am trying to construct the FDFD method for 3D structures. I have already constructed the general formulation and specifically I have set the complete matrix form. Due to the fact that the matrices are too sparse, and my system is out of memory, I am trying to set the problem in a dense mode.
To be more specific I want to transform the sparse matrices into dense form and to solve the problem efficiently. I have developed the problem in fortran, so I would like to ask you if anybody knows any library, which I can use in order to achieve this transformation. I know for example that in MATLAB using the sparse() routine I can achieve this goal. Is there any similar routine in fortran?

Thanks in advance!

 

[FredGarvin]

I don't know about the FORTRAN end of things, but depending on what your sparse matrix looks like, I know there are techniques for dealing with them. The first, and only real one that comes to mind, is in the case of a tri-diagonal or an almost tri-diagonal matrix.

 

[oswaler]

Hello there! It's great to hear that you are working on constructing the FDFD method for 3D structures. This is a widely used and powerful technique in computational electromagnetics. Regarding your question about transforming sparse matrices into dense form, there are a few options you can explore in Fortran. One option is to use the LAPACK library, which contains routines for dense matrix operations. Another option is to use the Intel Math Kernel Library (MKL), which also has efficient routines for dense matrix operations. Both of these libraries have Fortran interfaces and can be easily integrated into your code. I hope this helps and good luck with your research!