The solution of the system is
x( 1 ) = -4.186020128680939E-002
x( 2 ) = -3.413124159279098E-003
x( 3 ) = 0.117250376805018
x( 4 ) = -0.112639579923180
x( 5 ) = 2.417224446137141E-002
x( 6 ) = -0.107633340356223
x( 7 ) = 0.198719673273585
x( 8 ) = 0.190382963551205
Is this the right solution?
I also tried to call Pardiso myself in the file test.f90.
There I get the following output:
The array ia:
1 2 3 4 5 6
The array a:
(1.000000,0.0000000E+00) (1.000000,0.0000000E+00) (1.000000,0.0000000E+00)
The array ja:
1 2 3 4 5
The array x:
If I got the ja, ia and a arrays right, I don't understad the output.
I'm using Visual Studio 2008and the MKL 10.3 under Windows XP Professional SP2 .
OK pardiso_sym_f90.f90 works correct (of course). I made a mistake in testing it.
When I changed the type complexin test.f90 into complex*16 it gave me the right output. But how can I use single precision? In my application I only need single precision. Why is it not working with the type complex resp. where is my mistake?