pguerron

03-10-2009, 01:26 PM

Dear all,

I'm solving a system of equations using an eigenvalue decomposition technique. I use the ISML routines for eigenvctor (dgvcrg) and matrix inversion (dlincg). The result is almost correct except for one row in the solution matrix.

The code is enclosed, see main.f90. The solution starts in line number 174. The actual solution is stored in the matrices hxaux (line 308) and gxaux (line 310). This matrices differ with their counterparts from a solution in either matlab or mathematica in row #19. The rest of the solution is correct. I've written the code many times and get back to the same incorrect answer. Moreover, the logic behind the algorithm in fortran and matlab is the same.

Any suggestions what is wrong?

Thanks,

Pablo

I'm solving a system of equations using an eigenvalue decomposition technique. I use the ISML routines for eigenvctor (dgvcrg) and matrix inversion (dlincg). The result is almost correct except for one row in the solution matrix.

The code is enclosed, see main.f90. The solution starts in line number 174. The actual solution is stored in the matrices hxaux (line 308) and gxaux (line 310). This matrices differ with their counterparts from a solution in either matlab or mathematica in row #19. The rest of the solution is correct. I've written the code many times and get back to the same incorrect answer. Moreover, the logic behind the algorithm in fortran and matlab is the same.

Any suggestions what is wrong?

Thanks,

Pablo