Finite Difference Frequency Domain

AI Thread Summary
The discussion focuses on constructing the Finite Difference Frequency Domain (FDFD) method for 3D structures, specifically addressing issues with sparse matrices leading to memory limitations. The user seeks to convert these sparse matrices into a dense format for more efficient problem-solving in Fortran. They inquire about libraries or routines in Fortran similar to MATLAB's sparse() function for this transformation. Suggestions include exploring techniques for handling specific types of sparse matrices, such as tri-diagonal matrices. The conversation emphasizes the need for effective memory management in numerical computations.
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.
 

Thread 'Query on Isentropic relationships (such as PV^gamma = constant)'

Hi all, I have a question. So from the derivation of the Isentropic process relationship PV^gamma = constant, there is a step dW = PdV, which can only be said for quasi-equilibrium (or reversible) processes. As such I believe PV^gamma = constant (and the family of equations) should not be applicable to just adiabatic processes? Ie, it should be applicable only for adiabatic + reversible = isentropic processes? However, I've seen couple of online notes/books, and...
View full post »

Thread 'Dry sump pump issue'

I have an engine that uses a dry sump oiling system. The oil collection pan has three AN fittings to use for scavenging. Two of the fittings are approximately on the same level, the third is about 1/2 to 3/4 inch higher than the other two. The system ran for years with no problem using a three stage pump (one pressure and two scavenge stages). The two scavenge stages were connected at times to any two of the three AN fittings on the tank. Recently I tried an upgrade to a four stage pump...
View full post »

Similar threads

A Editing in Freq domain and applying inverse FFT
Replies
6
Views
2K
The 8 Mode Shapes of a plane finite element
Replies
2
Views
1K
MATLAB Turning Fourier coefficients into an interpolated freq domain function
Replies
8
Views
2K
I Boundary Conditions For Modelling of a Fluid Using Euler's Equations
Replies
16
Views
2K
Finite Difference Solution to Poisson's Equation on Irregular Domain
Replies
2
Views
2K
Unsteady vorticity transport equation: codes available?
Replies
2
Views
2K
A Numerical Hartree Fock with Finite Difference Matrices for Helium
Replies
4
Views
2K
Correct Setup for Finite Difference to Calculate Quantum Tunneling
Replies
2
Views
887
I E-field around charged finite wire - intractable or not?
Replies
10
Views
911
Fortran How Do You Implement MKL's FFT Functions for Multidimensional Data in Fortran?
Replies
2
Views
2K

Hot Threads

Recent Insights

Back
Top