Cassiopee Forum
CFD python modules
S'il vous plaît connectez vous ou Enregistrement'.
Date et heure en ce moment : mardi, 12 décembre 2017, 07:45

Cassiopee Forum :: Discussions, suggestions, bug report :: Post :: Bonjour J ai un probleme avec la fonction  ::
Eric Garnier
Newbie
Image


Messages: 2
Bonjour J ai un probleme avec la fonction (jeudi, 18 juillet 2013, 15:21)  
Bonjour

J ai un probleme avec la fonction extractpoint ou avec SplitSize


import Converter.PyTree as CV
import XTree as XT
import Post.PyTree as PS
import Converter.Internal as Internal
import Transform.PyTree as TF
#Je relis un maillage en sortie d ICEM
A = CV.convertFile2PyTree(Mesh)
A = CV.addVars(A, 'centers:cellN')
A = CV.initVars(A, 'centers:cellN', 1.)
listPoints = [(220., 0., 70.)]
PS.extractPoint(A, listPoints, order=2, tol=1.e-6)
[[1.0000000000000002]]
#Jusque la tout va bien il a trouvé que les coordonnées du point etaient dans mon maillage et il m a interpolé cellN
# Ensuite , je decoupe en bloc de 250000 points
A = TF.splitSize(A, 250000)
PS.extractPoint(A, listPoints, order=2, tol=1.e-6)
Warning: extractPoint: cannot interpolate point 220.0000006 0.0000006 70.0000006.
[[0.0]]
# La fonction extractPoint ne me retrouve plus le point a l interieur du domaine.

Cela ne me semble pas tres normal ?


eric

Ip enregistré Statut: déconnecté Profil | Site Web 
Ordre des réponses: Première réponse en dernier :: Première réponse en premier
StephaniePeron
Full Member
Image


Messages: 159
RE: Bonjour J ai un probleme avec la fonction (vendredi, 19 juillet 2013, 08:48)  
Effectivement.
Peux-tu me dire où trouver ton fichier "Mesh" ?
Quelle version de Cassiopee utilises-tu ?
Merci
a+

Ip enregistré Statut: déconnecté Profil | Site Web 
StephaniePeron
Full Member
Image


Messages: 159
RE: Bonjour J ai un probleme avec la fonction (vendredi, 19 juillet 2013, 09:56)  
En fait, le problème vient de ton champ en centres (cellN) :
les zones donneuses (A) sont converties de noeuds en centres.
Au niveau des raccords, il y a donc des trous.
Et le point reste tel quel, car la localisation centres/noeuds n'a pas de sens pour un pt. Il tombe donc dans le trou.
Avec un champ en noeuds, ça marche.
Exemple : A = CV.initVars(A,'{f}={CoordinateX}**2')
Si tu veux interpoler des champs en centres, il faut donner un maillage avec cellules fictives. Par ex, tu fais un
X.connectMatch(), suivi de C.addGhostCells(). Tu peux enlever les ghost cells avec C.rmGhostCells()
a+


Ip enregistré Statut: déconnecté Profil | Site Web 
Nouvelle discussion - Discussion barrée, pas de réponse possible

Software PBLang 4.65 © 2002-2003 by Martin Senftleben
Image