Welcome, Guest
Username: Password: Remember me

TOPIC: Mesh allocation error with large mesh, problem with mxptel.f ?

Mesh allocation error with large mesh, problem with mxptel.f ? 10 years 1 month ago #14427

  • pilou1253
  • pilou1253's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 584
  • Thank you received: 106
Hi all!

I am having trouble with a large mesh of ca 1 080 000 elements. The program stops abruptly in ALMESH. I made a traceback debugging by adding some WRITE in the subroutine, and it appears that the problem lies in the call to MXPTEL:
CALL MXPTEL(MXPTVS,MXELVS,IKLES,IELM,
     &              NPOIN,NELEM,NDP,IPOBO,.TRUE.)

I made a test by adding a WRITE before the first statements (after declarations) in mxptel.f which is not printed out. The problem should therefore be in one of the argument, but here stop my skills...

Any idea of what is happenning? This model has already been ran on another machine.

I have Win7 64bit, 24GO free disk space, 16GO RAM.
I use version 6.3.2.

I also made a looser mesh which worked normally.

Any help would be greatly appreciated!

Regards
PL
The administrator has disabled public write access.

Mesh allocation error with large mesh, problem with mxptel.f ? 10 years 1 month ago #14432

  • jmhervouet
  • jmhervouet's Avatar
Hello,

No idea! I you send me the case I can try to run it on various compilers and see if I can reproduce the problem.

Regards,

Jean-Michel Hervouet
The administrator has disabled public write access.

Mesh allocation error with large mesh, problem with mxptel.f ? 10 years 1 month ago #14440

  • pilou1253
  • pilou1253's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 584
  • Thank you received: 106
Hi!

Here comes the case:
dl.dropboxusercontent.com/u/36016083/PL.zip

Would be interesting to see if you can reproduce the problem.
I have gfortran and same model could be ran on machines also using gfortran.

Thanks in advance!

Regards
PL
The administrator has disabled public write access.

Mesh allocation error with large mesh, problem with mxptel.f ? 10 years 1 month ago #14454

  • jmhervouet
  • jmhervouet's Avatar
Hello,

Sorry I tried and it works here, with Nag compiler and with Intel compiler crammed with various checkings, so it is difficult to say what happens, but obviously it is a memory allocation error due to the mesh size. If others can try...

Regards,

JMH
The administrator has disabled public write access.
Moderators: pham

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