[CIG-SHORT] Problem running user compiled version of Relax

[Sylvain Barbot]

Hi Oliver,

I've never had this issue before. The code is optimized for gcc and
gfortran, but I compiled it successfully with ifort and the kml
library on my system. If you can, perhaps you can run it without the
queue system, or you can try with a different compiler. We have
successful builds with gcc-48 gcc-49 and gcc-5. ifort and mkl only
provide a factor of two performance over gcc with fftw.

Cheers,
Sylvain


On Thu, Jun 18, 2015 at 6:38 AM, Oliver Boyd <olboyd at usgs.gov> wrote:
> Hello,
> I compiled a version of relax using source obtained from the git repository
> (6/16/15) with the intel compilers and mkl library in order to run jobs on a
> cluster. Everything seemed to go fine during compilation. I wasn't sure how
> to install it so I simply copied relax to the parent directory of build in
> which resides setup.sh (not sure that this is a problem). I then ran
> setup.sh and ran an example. The example runs fine with the binary version
> available from the CIG website (v1.0.7).
>
> The configure options I use are
> ./waf configure --proj-dir=/usr/bin --gmt-dir=/usr/bin
> --gmt-incdir=/usr/include --gmt-libdir=/usr/lib64
> --mkl-libdir=/opt/intel/mkl/lib/intel64/ --mkl-incdir=/opt/intel/mkl/include
> --check-c-compiler=icc --check-fortran-compiler=ifort
> --check-cxx-compiler=icpc
>
> When I submit the job to the queueing system, I get
> #
> ----------------------------------------------------------------------------
> # RELAX: nonlinear postseismic relaxation with Fourier-domain Green's
> function
> #     * Intel MKL implementation of the FFT
> #     * parallel OpenMP implementation with 001/002 threads
> #     * export to longitude/latitude text format cancelled
> (--no-proj-output)
> #     * export to VTK format cancelled                     (--no-vtk-output)
> #
> ----------------------------------------------------------------------------
> # grid dimension (sx1,sx2,sx3)
>   512  512  512
> # sampling (dx1,dx2,dx3), smoothing (beta, nyquist)
>   1.00E+0  1.00E+0  1.00E+0  2.00E-1  2.00E+0
> # origin position (x0,y0) and rotation
>   0.00E+0  0.00E+0  0.00E+0
> # observation depth (displacement and stress)
>   0.00E+0  5.00E+0
> # output directory
>  ./NewMadridThrust_Afterslip_GMMh01 (time report:
> ./NewMadridThrust_Afterslip_GMMh01/time.txt)
> # lambda, mu, gamma (gamma = (1 - nu) rho g / mu)
>   3.50E+01  3.50E+01  5.68E-04
> # time interval, (positive time step) or (negative skip, scaling)
>   2 –5 0.5
> # output every   1.00E+1 time unit
> # number of observation planes
>     0
> # number of observation points
>     0
> # number of stress observation segments
>     0
> # number of prestress interfaces
>     0
> # number of linear viscous interfaces
>     0
> # number of nonlinear viscous interfaces
>     0
> # number of fault creep interfaces
>     1
> #
> ----------------------------------------------------------------------------
> # n    depth   gamma0 (a-b)sig friction cohesion
> #
> ----------------------------------------------------------------------------
> 001  0.00E+0  1.00E+0  1.20E+0  0.00E+0  0.00E+0
> # number of afterslip planes
>     1
> #
> ----------------------------------------------------------------------------
> # n        x1       x2       x3   length    width strike    dip   rake
> #
> ----------------------------------------------------------------------------
> 001   9.65E+1 -3.57E+1  1.50E+1  2.00E+2  1.61E+2  173.2   31.9   90.0
> # number of inter-seismic strike-slip segments
>     0
> # number of inter-seismic tensile segments
>     0
> # number of events
>     1
> # number of coseismic strike-slip segments
>     2
> #
> ----------------------------------------------------------------------------
> # n     slip       xs       ys       zs  length   width strike   dip   rake
> #
> ----------------------------------------------------------------------------
> 001  5.00E+0  1.89E+1 -1.05E+1  5.10E+0 4.00E+1 1.87E+1  173.2  31.9   90.0
> 002  6.63E+0  1.91E+1 -8.94E+0  0.00E+0 4.00E+1 5.33E+0  173.2  73.0   90.0
> #
> ----------------------------------------------------------------------------
> # number of coseismic tensile segments
>     0
> # number of coseismic dilatation point sources
>     0
> # number of surface loads
>     0
> # max sampling size (hor.,vert.): 1.60E+1 2.13E+0
> #
> ----------------------------------------------------------------------------
> forrtl: severe (174): SIGSEGV, segmentation fault occurred
> Image              PC                Routine            Line        Source
> relax              00000000005AF571  Unknown               Unknown  Unknown
> relax              00000000005ADCC7  Unknown               Unknown  Unknown
> relax              000000000055E794  Unknown               Unknown  Unknown
> relax              000000000055E5A6  Unknown               Unknown  Unknown
> relax              000000000050D184  Unknown               Unknown  Unknown
> relax              000000000051132D  Unknown               Unknown  Unknown
> libpthread.so.0    0000003F77C0F710  Unknown               Unknown  Unknown
> relax              0000000000454AB7  Unknown               Unknown  Unknown
> relax              00000000004143E1  Unknown               Unknown  Unknown
> relax              0000000000410C6E  Unknown               Unknown  Unknown
> libc.so.6          0000003F7701ED5D  Unknown               Unknown  Unknown
> relax              0000000000410B79  Unknown               Unknown  Unknown
> srun: error: compute87: task 0: Exited with exit code 174
>
> Please let me know what additional information you need.
>
> Thank you for your help,
> Oliver
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short