Finite Difference Frequency Domain

Click For Summary
SUMMARY

The discussion focuses on the implementation of the Finite Difference Frequency Domain (FDFD) method for 3D structures, specifically addressing the challenge of handling sparse matrices in Fortran. The user seeks a library or routine to convert sparse matrices into a dense format to resolve memory issues. They reference MATLAB's sparse() routine as a potential solution and inquire about equivalent functionality in Fortran. Techniques for managing tri-diagonal matrices are also mentioned as a possible approach.

PREREQUISITES
  • Understanding of Finite Difference Frequency Domain (FDFD) methods
  • Familiarity with matrix operations and sparse matrix representation
  • Proficiency in Fortran programming
  • Knowledge of numerical methods for solving linear systems
NEXT STEPS
  • Research Fortran libraries for sparse matrix manipulation, such as PETSc or Eigen
  • Explore techniques for converting sparse matrices to dense matrices in numerical computing
  • Learn about tri-diagonal matrix algorithms and their applications
  • Investigate memory management strategies in Fortran for large data structures
USEFUL FOR

Researchers and developers working on numerical simulations, particularly in computational physics or engineering, who require efficient matrix handling in Fortran.

la4361
Messages
1
Reaction score
0
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!
 
Engineering news on Phys.org
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.
 

Similar threads

Graduate Editing in Freq domain and applying inverse FFT
  • · Replies 6 ·
Replies
6
Views
2K
The 8 Mode Shapes of a plane finite element
  • · Replies 2 ·
Replies
2
Views
2K
MATLAB Turning Fourier coefficients into an interpolated freq domain function
  • · Replies 8 ·
Replies
8
Views
3K
Undergrad Boundary Conditions For Modelling of a Fluid Using Euler's Equations
  • · Replies 16 ·
Replies
16
Views
3K
Graduate Finite Difference Solution to Poisson's Equation on Irregular Domain
  • · Replies 2 ·
Replies
2
Views
2K
Graduate Numerical Hartree Fock with Finite Difference Matrices for Helium
  • · Replies 4 ·
Replies
4
Views
3K
Correct Setup for Finite Difference to Calculate Quantum Tunneling
  • · Replies 2 ·
Replies
2
Views
1K
Undergrad E-field around charged finite wire - intractable or not?
  • · Replies 10 ·
Replies
10
Views
2K
Question about an LTI System and its Frequency Response
  • · Replies 2 ·
Replies
2
Views
2K
Fortran How Do You Implement MKL's FFT Functions for Multidimensional Data in Fortran?
  • · Replies 2 ·
Replies
2
Views
2K