|June 14, 2011, 11:50||
URGENT: Please Help! Fortran 90 code for bicgstab algorithm
Join Date: Jun 2011
Location: Bordeaux, France
Posts: 1Rep Power: 0
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.
Last edited by CuriousBoy; June 16, 2011 at 08:18.
|Thread||Thread Starter||Forum||Replies||Last Post|
|fortran code needed||HaKu||Main CFD Forum||2||February 15, 2010 13:02|
|How to parallelize this fortran code by openmp||ronac||Main CFD Forum||0||January 18, 2008 22:36|
|3D fortran Navier-Stokes code||emrah||Main CFD Forum||4||December 15, 2007 10:53|
|Help on fortran code for FD 4th order||Miguel||Main CFD Forum||0||November 14, 2005 11:58|
|Design Integration with CFD?||John C. Chien||Main CFD Forum||19||May 17, 2001 15:56|