- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

In the file PROGRAMA07021011PIVOTGAUSS.f90 I'm usinggetrfand getrs.In the file PROGRAMA07022011SPARSE.f90 I'm resolved by sparse matrices.The sparse file, once run

generates five files val.dat, col.dat, rowindex.dat, ld.dat and sp.txt.SP.txt file is the matrix Ax = b, A system.ld.txt is the vector b.While val.dat are non-zero values of A,

rowIndex col.dat your column and row index.I have found these three on hand, and appear well.

The solution is given in the file solutionsparse.dat and are huge numbers.The solution of the file PROGRAMA07021011PIVOTGAUSS.f90 (apparently correct) iscalled solutiongauss.dat.

I tried to solve the system with matlab, inverting the matrix A and multiplying by solutiongauss.dat by match.

Does anyone know why I get these results?

Compile files with this line:

ifort-openmp-O3-o $ MKLINCLUDE PROGRAMA07022011SPARSE-I-L $ MKLPATH PROGRAMA07022011SPARSE.f90-lmkl_lapack95-Wl, - start-group $ MKLPATH / libmkl_intel.a $ MKLPATH / libmkl_intel_thread.a $ MKLPATH / libmkl_core.a-Wl,- end-group-lpthread-liomp5

and with:

ifort-openmp-O3-o $ MKLINCLUDE PROGRAMA07022011PIVOTGAUSS-I-L $ MKLPATH PROGRAMA07022011PIVOTGAUSS.f90-lmkl_lapack95-Wl, - start-group $ MKLPATH / libmkl_intel.a $ MKLPATH / libmkl_intel_thread.a $ MKLPATH / libmkl_core.a-Wl,- end-group-lpthread-liomp5

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Link Copied

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

*Hi,Your matrix is almost singular:Singular values:max= 52.7334033722656 min= 3.009615796461843E-014Condition number = 1.752163961734249E+015Theoretically, solving such matrices in double precision cannot produce stable results.Results,Konstantin*

Moreover, in the archives of triangulation would have a bug, right?

Thanks for your help

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

*PROGRAMA09022011.f90(654): error #6633: The type of the actual argument differs from the type of the dummy argument.*

if (dss_solve_real( hDSS, MKL_DSS_FORWARD_SOLVE, P, 1, X) &

-------------------------------------------------------^

PROGRAMA09022011.f90(661): error #6633: The type of the actual argument differs from the type of the dummy argument.

if (dss_solve_real( hDSS, MKL_DSS_FORWARD_SOLVE, P, 1, X) &

-------------------------------------------------------^

PROGRAMA09022011.f90(661): error #6633: The type of the actual argument differs from the type of the dummy argument.

if (dss_solve_real( hDSS, MKL_DSS_BACKWARD_SOLVE, P, 1, X ) &

--------------------------------------------------^

PROGRAMA09022011.f90(661): error #6633: The type of the actual argument differs from the type of the dummy argument.

if (dss_solve_real( hDSS, MKL_DSS_BACKWARD_SOLVE, P, 1, X ) &

--------------------------------------------------------^

PROGRAMA09022011.f90(668): error #6633: The type of the actual argument differs from the type of the dummy argument.

if (dss_solve_real( hDSS, MKL_DSS_DEFAULTS, P, 1, X ) &

--------------------------------------------^

PROGRAMA09022011.f90(668): error #6633: The type of the actual argument differs from the type of the dummy argument.

if (dss_solve_real( hDSS, MKL_DSS_DEFAULTS, P, 1, X ) &

--------------------------------------------------^

compilation aborted for PROGRAMA09022011.f90 (code 1)

**Any idea how to calculate my system sparse matrix?**

Thanks, for your help.

Regards

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

The program will have to solve very large sparse systems, from a finite difference scheme. For this case, I can resolution with getrfand getrs, but when I grow up I'll have to calculate the system with sparse.

In theory, Is it normal that the system of good results using getfrf (and getrs) and bad results using sparse?

On Other hand, How matrix ill-conditioned test them? For viewing the finite difference schemes should have a good ill-conditioned?

Sorry for the inconvenience

Thanks and Regards!

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

How I can get dss values, as I do with the file SP.txt (write (52.52)?

Thanks, and I'm sorry for the inconvenience

Regards

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Anyway, thank you very much Konstantine

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page