Hi,
We are working with Jérôme on the Hmat features for kriging. However, we
noticed that the binaries that come from conda-forge are not compiled with
HMAT. This is revealed by the following script :
import openturns as ot
covModel = ot.SquaredExponential(1)
myRegularGrid = ot.RegularGrid(0.0, 0.1, 100)
covarianceMatrix = covModel.discretize(myRegularGrid)
ot.ResourceMap_SetAsBool("HMatrix-ForceSequential", True)
ot.ResourceMap_SetAsUnsignedInteger("HMatrix-MaxLeafSize", 10)
hmatrixParameters = ot.HMatrixParameters()
hmatrixParameters.setAssemblyEpsilon(1.e-6)
hmatrixParameters.setRecompressionEpsilon(1.e-6)
nuggetFactor = 1.e-5
covarianceMatrix =
covModel.discretizeHMatrix(myRegularGrid,nuggetFactor,hmatrixParameters)
which is partially adapted from :
https://github.com/openturns/openturns/blob/master/lib/test/t_HMatrix_LU.cxx
The script produces
RuntimeError: NotYetImplementedException : In
CovarianceModelImplementation::discretizeHMatrix, OpenTURNS had been compiled
without HMat support
This leads to 2 questions :
· We used the following statement inside a kriging script and noticed
that nothing changed :
ResourceMap.Set("KrigingAlgorithm-LinearAlgebra", "HMAT")
This is bizarre : we might expect that the lack of HMAT support would generate
an exception or a warning. Instead, the script worked fine, even though we
actually did not use hmat.
· Is there a way of getting OpenTURNS Anaconda binaries with HMAT
enabled ? The workaround is to compile OpenTURNS, but this is inconvenient for
many users.
Best regards,
Michaël
Michaël BAUDIN
Ingénieur - Chercheur
EDF - R&D
Département PRISME
6, quai Watier
78401 CHATOU
[email protected]<mailto:[email protected]>
Tél. : 01 30 87 81 82
Fax : 01 30 87 82 13
Ce message et toutes les pièces jointes (ci-après le 'Message') sont établis à
l'intention exclusive des destinataires et les informations qui y figurent sont
strictement confidentielles. Toute utilisation de ce Message non conforme à sa
destination, toute diffusion ou toute publication totale ou partielle, est
interdite sauf autorisation expresse.
Si vous n'êtes pas le destinataire de ce Message, il vous est interdit de le
copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou partie. Si
vous avez reçu ce Message par erreur, merci de le supprimer de votre système,
ainsi que toutes ses copies, et de n'en garder aucune trace sur quelque support
que ce soit. Nous vous remercions également d'en avertir immédiatement
l'expéditeur par retour du message.
Il est impossible de garantir que les communications par messagerie
électronique arrivent en temps utile, sont sécurisées ou dénuées de toute
erreur ou virus.
____________________________________________________
This message and any attachments (the 'Message') are intended solely for the
addressees. The information contained in this Message is confidential. Any use
of information contained in this Message not in accord with its purpose, any
dissemination or disclosure, either whole or partial, is prohibited except
formal approval.
If you are not the addressee, you may not copy, forward, disclose or use any
part of it. If you have received this message in error, please delete it and
all copies from your system and notify the sender immediately by return message.
E-mail communication cannot be guaranteed to be timely secure, error or
virus-free.
_______________________________________________
OpenTURNS users mailing list
[email protected]
http://openturns.org/mailman/listinfo/users
