hello,

I've been trying to use the Multigrid PDE solver mgfas, but bumped upon the following problem in the rstrct.f90 function, which begins with :

FUNCTION rstrct(uf)
USE nrtype; USE nrutil, ONLY : assert_eq
IMPLICIT NONE
REAL(DP), DIMENSION(:,:), INTENT(IN) :: uf
REAL(DP), DIMENSION((size(uf,1)+1)/2,(size(uf,1)+1)/2) :: rstrct
nc=(size(uf,1)+1)
rstrct(2:nc-1,2:nc-1)= 0.0

etc.

as soon as the program reaches the last statement ( filling in the matrix rstrct), it Segmentation Faults. What is wrong ?

Thank you for your help,

Nicolas

You can all stop scratching your heads, it was a compiler issue. I installed the intel compiler instead of the gnu one and it worked. Don't ask me why.

Nicolas