> I find that both nf.f90 and capacita.f90 segfault in runtime for any stack 
> size.

On x86_64-apple-darwin10, nf.f90 "works". However if I run it through
valgrind I get

==64815== Memcheck, a memory error detector
==64815== Copyright (C) 2002-2010, and GNU GPL'd, by Julian Seward et al.
==64815== Using Valgrind-3.6.1 and LibVEX; rerun with -h for copyright info
==64815== Command: a.out --max-stackframe=2118496
==64815==
==64815== Warning: set address range perms: large range [0x7ffe6c000000, 
0x7fff5bc01000) (defined)
==64815== Warning: client switching stacks?  SP change: 0x7fff5bffe410 --> 
0x7fff5be0cef0
==64815==          to suppress, use: --max-stackframe=2037024 or greater
==64815== Invalid write of size 8
==64815==    at 0x100003B22: nf2dprecon.1828 (nf.f90:262)
==64815==    by 0x10000C3B6: nfcg_ (nf.f90:279)
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==  Address 0x7fff5be0cee8 is on thread 1's stack
==64815==
==64815== Invalid read of size 8
==64815==    at 0x100001296: trisolve.1839 (nf.f90:238)
==64815==    by 0x10000C3B6: nfcg_ (nf.f90:279)
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==  Address 0x7fff5be0cee8 is on thread 1's stack
==64815==
...
==64815== Invalid write of size 8
==64815==    at 0x100000F38: trisolve.1839 (nf.f90:231)
==64815==    by 0x10000C3B6: nfcg_ (nf.f90:279)
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==    by 0x3FEFFFFFFFFFFFFF: ???
==64815==  Address 0x7fff5bffdde8 is on thread 1's stack
==64815==
==64815== Invalid read of size 8
==64815==    at 0x100000F64: trisolve.1839 (nf.f90:233)
==64815==    by 0xC02396760B69D62B: ???
==64815==    by 0xC02396760B69D62C: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62D: ???
==64815==    by 0xC02396760B69D62B: ???
==64815==    by 0xC02396760B69D62B: ???
==64815==    by 0xC02396760B69D62B: ???
==64815==  Address 0x7fff5bffdde8 is on thread 1's stack
==64815==
...
 Time for setup          1.316
 Time per iteration      3.250
 Total Time            144.295
==64815==
==64815== HEAP SUMMARY:
==64815==     in use at exit: 656 bytes in 13 blocks
==64815==   total heap usage: 411 allocs, 398 frees, 387,314,592 bytes allocated
==64815==
==64815== LEAK SUMMARY:
==64815==    definitely lost: 0 bytes in 0 blocks
==64815==    indirectly lost: 0 bytes in 0 blocks
==64815==      possibly lost: 0 bytes in 0 blocks
==64815==    still reachable: 656 bytes in 13 blocks
==64815==         suppressed: 0 bytes in 0 blocks
==64815== Rerun with --leak-check=full to see details of leaked memory
==64815==
==64815== For counts of detected and suppressed errors, rerun with: -v
==64815== ERROR SUMMARY: 40050 errors from 1000 contexts (suppressed: 0 from 0)

The segfault for capacita.f90 occurs in the subroutine fourir at the line

            write(unit=*, fmt=*) "error in fourier: n=", ntot

AFAICT the problem occurs in the loop

      do m=1,ntot/4-1
        E(m) = exp(m*h)
      end do

If I print ntot, loc(ntot) before it I get

        2048      140734799794712

After the loop loc(ntot) is

              9205357642636066816

and any attemp to print its values yields a segfault.

Dominique

Reply via email to