Hi Adam,<br><br>Could you let me know what version of the gfortran compiler you are using? On our system here, we aren't seeing the genmap errors. Also, if they are pointing to undeclared arrays, we would probably want to go ahead and declare those arrays too.<br>
<br>I am still adjusting the 2d_eig tests in Analysis.py. Thanks for the input and soon that test and the MOAB tests should be fully functional. As Aleks mentioned, we are having an issue with the MHD case and versions after v780.<br>
<br>Would you send me your scripts off list (<a href="mailto:heisey@mcs.anl.gov">heisey@mcs.anl.gov</a>) with the different path structure? If it makes it easier to test in other environments, I think we would be interested in adding your changes to the repo.<br>
<br>Thanks so much!<br>Katie<br><br><div class="gmail_quote">On Fri, Aug 24, 2012 at 10:29 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Adam,<br>
<br>
Thanks for the detailed email -- I'll discuss it with developers.<br>
<br>
But let me comment on certain issues you have raised.<br>
<br>
We are close to finalizing MOAB tests so they will be added shortly.<br>
<br>
There are some issues with semg MHD cases after revision 780. Additive Shwatrz (p43=1) works fine, and it is good to know that FDM works too.<br>
<br>
Thanks.<br>
Aleks<br>
<div><div><br>
<br>
<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
To: "nekton" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>><br>
Sent: Friday, August 24, 2012 8:58:09 AM<br>
Subject: [Nek5000-users] testing suite<br>
<br>
Hi Neks<br>
<br>
I have a few questions about test suite in nekton (in nek5_svn/tests/).<br>
Nek5000 is a part of CRESTA project and I need simple method for testing<br>
new versions of the code on different machines. Tests in nek5_svn/tests/<br>
are perfect, but I've got some problems:<br>
<br>
1) genmap work perfect with pgi compiler, but with gnu compilers returns<br>
segmentation fault for conj_ht test. It's due to not initialized key<br>
array in subroutine spec_bis_conn. It would be good to add<br>
-finit-local-zero flag for gfortran in maketools<br>
<br>
2) tests in 2d_eigtest work fine for pgi, but fail for gfortran, since<br>
the difference in the lambda error is bigger than 10^-6. It happens for<br>
N=3 and 4 only; eig1.err:<br>
GNU:<br>
2 19 0 5.7831860E+00 2.1065324E-02 1.6329577E-01 lambda<br>
3 34 0 5.7831860E+00 6.8588903E-04 1.3936160E-01 lambda<br>
4 61 0 5.7831860E+00 6.4127192E-07 2.6444060E-02 lambda<br>
5 84 0 5.7831860E+00 4.7561843E-08 3.3011994E-02 lambda<br>
6 121 0 5.7831860E+00 3.2211656E-10 1.1637760E-02 lambda<br>
7 178 0 5.7831860E+00 2.8581582E-12 1.6504072E-03 lambda<br>
8 243 0 5.7831860E+00 5.4534155E-13 4.8188495E-04 lambda<br>
9 314 0 5.7831860E+00 8.6153307E-14 2.9147672E-05 lambda<br>
10 400 0 5.7831860E+00 3.2018832E-12 6.8264261E-06 lambda<br>
11 478 0 5.7831860E+00 7.3994144E-12 4.4899419E-07 lambda<br>
12 614 0 5.7831860E+00 6.7972294E-12 4.0358872E-08 lambda<br>
13 715 0 5.7831860E+00 2.7814195E-11 2.0542709E-09 lambda<br>
14 829 0 5.7831860E+00 1.8316904E-11 1.1108244E-10 lambda<br>
15 883 0 5.7831860E+00 1.1938006E-11 3.0416639E-12 lambda<br>
16 900 0 5.7831860E+00 1.6848745E-12 5.3768924E-12 lambda<br>
PGI:<br>
2 19 1 5.7831860E+00 2.1065324E-02 1.6329577E-01 lambda<br>
3 37 1 5.7831860E+00 6.8588903E-04 1.4437432E-01 lambda<br>
4 61 1 5.7831860E+00 6.4127192E-07 2.8477333E-02 lambda<br>
5 89 1 5.7831860E+00 4.7561947E-08 3.3011994E-02 lambda<br>
6 119 1 5.7831860E+00 3.2096015E-10 1.1637760E-02 lambda<br>
7 178 1 5.7831860E+00 2.7320368E-12 1.6504072E-03 lambda<br>
8 235 1 5.7831860E+00 1.0373924E-12 4.8188495E-04 lambda<br>
9 313 1 5.7831860E+00 5.7545080E-12 2.9147672E-05 lambda<br>
10 400 1 5.7831860E+00 1.0382806E-12 6.8264260E-06 lambda<br>
11 477 1 5.7831860E+00 2.2977176E-12 4.4899373E-07 lambda<br>
12 614 1 5.7831860E+00 2.7853275E-12 4.0358850E-08 lambda<br>
13 713 1 5.7831860E+00 1.4867219E-11 2.0542157E-09 lambda<br>
14 826 1 5.7831860E+00 1.4868107E-11 1.1256349E-10 lambda<br>
15 884 1 5.7831860E+00 1.1343815E-11 3.8686936E-12 lambda<br>
16 900 1 5.7831860E+00 5.1991300E-11 1.1914564E-11 lambda<br>
<br>
Have you got the same result? Should the tolerance in Analysis.py be<br>
bigger in this case? Moreover this test was failing for me with<br>
<br>
[Example 2d_eig/INT2: Serial-iter/err]...I couldn't find all the<br>
requested value in the log file...<br>
[Example 2d_eig/INT2: Serial-iter/err]... 2 , 3 , 4 , 10 ,<br>
were not found...<br>
Example 2d_eig/INT2: Serial-iter/err : F<br>
<br>
due to single list<br>
value = [[' 2 ',19,2,6],<br>
[' 2 ',1.6329577E-01,1e-06,2],<br>
[' 3 ',37,2,6],<br>
[' 3 ',1.4437432E-01,1e-06,2],<br>
[' 4 ',61,2,6],<br>
[' 4 ',2.8477333E-02,1e-06,2],<br>
[' 10 ',400,2,6],<br>
[' 10 ',6.8264260E-06,1e-06,2]]<br>
<br>
in Analysis.py. To get rid of this I had to divide it to 2 lists<br>
value = [[' 2 ',19,2,6],<br>
[' 3 ',37,2,6],<br>
[' 4 ',61,2,6],<br>
[' 10 ',400,2,6]]<br>
<br>
value = [[' 2 ',1.6329577E-01,1e-06,2],<br>
[' 3 ',1.4437432E-01,1e-06,2],<br>
[' 4 ',2.8477333E-02,1e-06,2],<br>
[' 10 ',6.8264260E-06,1e-06,2]]<br>
<br>
3) mhd-gpf (case without fast diagonalization method lfdm=0) fails both<br>
for pgi and gnu compilers. It seems pressure does not converge giving<br>
nan for gfortran. Cases with lfdm=1 give no errors.<br>
<br>
4) I found most of moab tests commented in Analysis.py. Does it mean<br>
moab is not fully implemented yet?<br>
<br>
5) I found testing scripts tuned for single environment (./ in the PATH,<br>
moab and matlab locations hidden in the script), so I've modified them<br>
slightly to make them easier to run for others. Is it possible to add my<br>
version of scripts to the repository?<br>
<br>
Regards<br>
<br>
Adam<br>
<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br>