[CuriousBoy]

Hello Everybody,

I am solving a set of linear equations arising from the FVM discretisation of the 1D Poisson Equation which is part of a larger problem. In my case, the pressure is dependent on the density, which in turn is dependent on the temperature.

The matrix is not symmteric due to the boundary conditions imposed on the temperature field. As it also turns out, the matrix is not diagonal domnant but has the property ap(i) = ae(i) + aw(i), where ap is the main diagonal element, ae is the superdiagonal element and aw is the subdiagonal element in the i-th row. As a result, TDMA, Gauss-Seidel, Jacobi, SOR are not working.

My entire program is in Fortran 90. I am searching for a Fortran 90 BiCGSTAB code to solve the same problem. I tried with the commonly available F77 codes but apparently, they are not interfacing well with F90.

I tried to obtain SMLIB, but the website appears to be out of order.

I will be very greatful or your help in this matter. Thank you in advance.

[ankeshgoyal]

i also trying to find the sparskit.f90 code but i am unable to find it .thats why i am working on f77 code of sparskit given by yousef saad. i need to parallelize this code. as i tried my best to parallelize this code by parallelizing the bcgstab solver (subroutine).but it is not getting parallelized can some one help me to parallelize this code with open mp directives . thanks in advance