====== Présentation ======

  * [[http://hsl.rl.ac.uk/archive/hslarchive.html|Site HSL]]
  * soumis à [[http://www.hsl.rl.ac.uk/hsl2007/academic_licence.html|licence]], donc impossible de l'installer sur un espace commun...
  * [[http://mumps.enseeiht.fr/doc_n7/ma41R.pdf|ma41ad]]: solveur direct (direct sparse solver for unsymmetric matrices)
  * [[http://www.hsl.rl.ac.uk/specs/mi26.pdf|mi26ad]]: solveur itératif (BiCGStab)

====== sur nemo ======

  * ./projects/tests/hsl

fichiers fournis par J.D:
<code>
nemo-henry% sccs diffs Makefile

------- Makefile -------
4,5c4,5
< FORTRAN  = gfortran
< OPTS     = -O5
---
> FORTRAN  = f90
> OPTS     = -O5 -m64

nemo-henry% module list
Currently Loaded Modulefiles:
  1) ss12u1   2) hpc82
nemo-henry% make |& tee MAKE.LOG
f90 -O5 -m64 -c depf77.for
f90 -O5 -m64 -c hsl_zd11_double.f90
f90 -O5 -m64 -c hsl_ma48_double.f90
f90 -O5 -m64 -c mi26ad.for
f90 -O5 -m64 -c ma41ad.for
f90 -O5 -m64 -c mi24ad.for
f90 -O5 -m64 -c mi11ad.for
ar cr libhsl.a depf77.o hsl_zd11_double.o hsl_ma48_double.o mi26ad.o ma41ad.o mi24ad.o mi11ad.o
ranlib libhsl.a

Copier  libhsl.a  dans /usr/local/lib/, les .mod dans /usr/local/lib/mods/
et utiliser '-fintrinsic-modules-path/usr/local/lib/mods' comme option de gfortran
pour compiler les programmes utilisant cette librairie

</code>
sur nemo, on fera:
  * Copier  libhsl.a dans /local/apps/lib/amd64/lib/
  * les .mod dans /local/apps/lib/amd64/lib/mods/
  * utiliser '-M/local/apps/lib/amd64lib/mods' ou MODDIR environment variable [[http://docs.sun.com/app/docs/doc/820-7600/aevog?a=view|4.9.1 Searching for Modules]]

puis:
<code>
nemo-henry% make -f Makefile.hmm
f90 -O5 -M ../droniou/hsl2007 -m64 -c mod_interface_hsl.f90
f90 -O5 -M ../droniou/hsl2007 -m64 -c benchencours.f90
f90 -O5 -M ../droniou/hsl2007 -m64 benchencours.o mod_interface_hsl.o -R../droniou/hsl2007 -L../droniou/hsl2007 -lhsl -o BENCHENCOURS
nemo-henry% file BENCHENCOURS
BENCHENCOURS:   ELF 64-bit LSB executable AMD64 Version 1 [SSE2 SSE FXSR CMOV FPU], dynamically linked, not stripped
</code>

====== Un cas test ======

<code>
============================
imesh =            5
============================
 entre  lecture maillage imesh =            5
 typemail=           -1
Mesh  ../MESHES/hexa_32x32x32_locraf.msh                                             
 nom du maillage grid_5_good                                                                    
 taille voisins,faces/volumes       370944
 taille sommets/volume       491520
 taille sommets/face       758784
 coefficient de planarite des faces   1.00000000000000000     volume total   1.00000000000000000     ok lecture maillage
 ok lecturemail
 ok calcul des ask :  nhyb         8448  nfbar       181248  nlin        69888
 ok calculask
 ok calcul numeq
 lcvn =           17
 nptmat avant      3379987
 nptmat apres      1680292
 ok calculptmat
 ok calculcoef
 ok calculsl
Operating system error: Cannot allocate memory
Out of memory 
</code>