Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
ARTEMIS
[v7p2r0] Unable to run test case chwac1
Flat Threaded

TOPIC: [v7p2r0] Unable to run test case chwac1

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25295

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
OK - I can confirm that chwac1 works both in scalar mode and in parallel (including the use of MUMPS).

Let me know if you still have troubles.

Sébastien.
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25304

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
I am afraid I am still having the same issue with the trunk. Both scalar and parallel. Other examples work fine.
Here is my configuration file:

File Attachment:

File Name: systel_gfmsmpi.cfg
File Size: 5 KB

Any help would be appreciated.
Regards,
Costas
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25305

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
Just in case - please make sure you remove the local executable (when you switch from config to another).

In fact you can just have the parallel version on your system and use it with 1 or more processor.

Also, remember to change SOLVER = 8 into 9 for the parallel.

Séabstien.
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25312

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Hello Sebastien,

I usually check and remove older executables as part of the troubleshooting process. Regarding the solver, ARTEMIS throws an error if I mix it up.
However, I just noticed that ARTEMIS has now a DEBBUGER option available. So the latest run gives:

File Attachment:

File Name: output_2017-02-20.txt
File Size: 15 KB

Which points to LECWAC1 for errors.
I have added a WRITE statement (WRITE(LU,*) 'LECWAC: START') at the beginning of LECWAC1 (which btw has 2 "USE DECLARATIONS_SPECIAL" statements) as such:
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      USE DECLARATIONS_TELEMAC
      USE DECLARATIONS_ARTEMIS
      USE INTERFACE_HERMES
      USE DECLARATIONS_SPECIAL
      USE BIEF, ONLY : FIND_VARIABLE
!
      USE DECLARATIONS_SPECIAL
      IMPLICIT NONE
!
      CHARACTER(LEN=16),ALLOCATABLE :: VARNAME(:),VARUNIT(:)
      INTEGER      NBI1,NBCL
      INTEGER      IPLAN,IFF,ID,KK,J,IDEC,JD
      DOUBLE PRECISION :: COORD(MAXDIR*MAXFRE)
      REAL         CL1(MAXDIR,MAXFRE,1)
      REAL         XTWC(MAXDIR,MAXFRE),YTWC(MAXDIR,MAXFRE)
      REAL         STOCKD(MAXDIR)
      REAL         PI,DEGRAD,EPS
      DOUBLE PRECISION :: READTWC(MAXDIR*MAXFRE)
      INTEGER :: IERR
      CHARACTER(LEN=8) :: FFORMAT
!
!-----------------------------------------------------------------------
!
      WRITE(LU,*) 'LECWAC: START'
!> SEB @ HRW: ALGORITHMIC DIFFERENTIATION
      PI = 4.D0 * ATAN( 1.D0 )
      DEGRAD = PI / 180.D0
!      PARAMETER(PI = 3.1415926535897932384626433D0 , DEGRAD=PI/180.D0)
!< SEB @ HRW
!
!-----------------------------------------------------------------------
However, it does not print. Is it possible to have an error in the declarations (that is not announced by the compiler)?

Regards,
Costas
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25313

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
... I am not sure about the compilation -- I know intel lets you get away with mistakes.
Try taking the USE DECLARATIONS_SPECIAL out (one of them) -- or add that file to your PRINCI for further debugging.

Let me know.
Sébastien
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25314

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
I already included several core subroutines in my princi file in order to add debbuging messages until I noticed they were already there... Removing either or both USE DECLARATIONS_SPECIAL did not make any difference.

What compiler do you use for the trunk validation?

Costas
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25315

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
We use gfortran 6.2.

Yo ucan try with the following options:
gfortran -c -g -cpp -fbounds-check -Wall -fbacktrace -finit-real=nan

For windows you can download a good install from here:
www.equation.com
The administrator has disabled public write access.
The following user(s) said Thank You: cyamin

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25316

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Gotcha! The compilation options you posted worked for the scalar debug configuration.
Here is a typical parallel configuration I use, which apparently is a bit outdated.
cmd_obj:    x86_64-w64-mingw32-gfortran -cpp -march=nehalem -O3 -fopenmp -c -fno-range-check -ffixed-line-length-132 -fconvert=big-endian -frecord-marker=4 -DHAVE_MPI -DHAVE_MUMPS <mods> <incs> <f95name>
cmd_exe:    x86_64-w64-mingw32-gfortran -static -march=nehalem -O3 -fopenmp -fconvert=big-endian -frecord-marker=4 -v -o <exename> <objs> -Xlinker --start-group <libs>
What would be the up-to-date commands?
Costas
The administrator has disabled public write access.

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25317

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
My configuration on windows is as follows (I am still not able to compile MUMPS on windows).
cmd_obj:    gfortran -c -Og -cpp -fbounds-check -ffpe-trap=zero,overflow,underflow -Wall -fbacktrace -finit-real=nan -fconvert=little-endian -DHAVE_MPI -frecord-marker=4 <mods> <incs> <f95name>
cmd_exe:    gfortran -fconvert=little-endian -frecord-marker=4 -v -lm -lz -o <exename> <objs> <libs>

On Linux (CentOS7) I have:
cmd_obj:    gfortran -c -g -cpp -fconvert=big-endian -DHAVE_MPI -DHAVE_MUMPS -frecord-marker=4 <mods> <incs> <f95name>
cmd_exe:    mpif90 -fconvert=big-endian -frecord-marker=4 -v -lm -lz -o <exename> <objs> <libs>

hope this helps.

Sébastien.
The administrator has disabled public write access.
The following user(s) said Thank You: huyquangtran

[v7p2r0] Unable to run test case chwac1 7 years 9 months ago #25318

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
... remove -fopenmp to see.
The administrator has disabled public write access.
Flat Threaded
Forum
System Modules
ARTEMIS
[v7p2r0] Unable to run test case chwac1

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.