Version: 7.8.0
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Groups Pages
example19

SALOMEDS methods : FindOrCreateAttribute, Save, SaveAs, Close, Open, FindComponent, FindAttribute, FindObject, LoadWith
Attribute methods : Value, SetValue

#======================================================================
#1. Create Study and add some components to it
#======================================================================
#=================================
# create AttributeReal
#=================================
A = batchmode_geompy.myBuilder.FindOrCreateAttribute(batchmode_geompy.father, "AttributeReal")
if A == None :
raise RuntimeError, "Can't create AttributeReal attribute"
A = A._narrow(SALOMEDS.AttributeReal)
A.SetValue(0.0001)
if A.Value() != 0.0001:
raise RuntimeError, "Error : wrong value of AttributeReal"
# === Geometry ==================================================
import batchmode_geompy
#===================================
# define a box
#===================================
box = batchmode_geompy.MakeBox(0., 0., 0., 100., 200., 300.)
idbox = batchmode_geompy.addToStudy(box,"box")
geom = batchmode_geompy.geom
import SMESH
import batchmode_smesh
ShapeType = batchmode_smesh.ShapeType
# ---- add first face of box in study
subShapeList = batchmode_geompy.SubShapeAll(box,ShapeType["Face"])
face=subShapeList[0]
name = "box_face"
print name
idface= batchmode_geompy.addToStudyInFather(box,face,name)
# ---- add shell from box in study
subShellList= batchmode_geompy.SubShapeAll(box,ShapeType["Shell"])
shell = subShellList[0]
name = "box_shell"
print name
idshell= batchmode_geompy.addToStudyInFather(box,shell,name)
# ---- add first edge of face in study
edgeList = batchmode_geompy.SubShapeAll(face,ShapeType["Edge"])
edge=edgeList[0];
name = "face_edge"
print name
idedge= batchmode_geompy.addToStudyInFather(face,edge,name)
# ---- launch SMESH, init a Mesh with the box
smesh = batchmode_smesh.smesh
# -- Init --
shape = batchmode_geompy.IDToObject(idbox)
mesh=smesh.Init(geom, batchmode_geompy.myStudyId, shape)
orb = batchmode_geompy.orb
idmesh = batchmode_smesh.AddNewMesh( orb.object_to_string(mesh) )
batchmode_smesh.SetName(idmesh, "Meshbox");
batchmode_smesh.SetShape(idbox, idmesh);
# ---- create Hypothesis
print "-------------------------- create Hypothesis ----------------------"
print "-------------------------- LocalLength"
hyp1 = smesh.CreateHypothesis("LocalLength", batchmode_geompy.myStudyId )
hypLen1 = hyp1._narrow(SMESH.SMESH_LocalLength)
hypLen1.SetLength(100)
print hypLen1.GetName()
print hypLen1.GetId()
print hypLen1.GetLength()
idlength = batchmode_smesh.AddNewHypothesis( orb.object_to_string(hypLen1) );
batchmode_smesh.SetName(idlength, "Local_Length_100");
print "-------------------------- NumberOfSegments"
hyp2 = smesh.CreateHypothesis("NumberOfSegments", batchmode_geompy.myStudyId )
hypNbSeg1=hyp2._narrow(SMESH.SMESH_NumberOfSegments)
hypNbSeg1.SetNumberOfSegments(7)
print hypNbSeg1.GetName()
print hypNbSeg1.GetId()
print hypNbSeg1.GetNumberOfSegments()
idseg = batchmode_smesh.AddNewHypothesis( orb.object_to_string(hypNbSeg1) );
batchmode_smesh.SetName(idseg, "NumberOfSegments_7");
print "-------------------------- MaxElementArea"
hyp3 = smesh.CreateHypothesis("MaxElementArea", batchmode_geompy.myStudyId)
hypArea1=hyp3._narrow(SMESH.SMESH_MaxElementArea)
hypArea1.SetMaxElementArea(2500)
print hypArea1.GetName()
print hypArea1.GetId()
print hypArea1.GetMaxElementArea()
idarea1 = batchmode_smesh.AddNewHypothesis( orb.object_to_string(hypArea1) );
batchmode_smesh.SetName(idarea1, "MaxElementArea_2500");
print "-------------------------- MaxElementArea"
hyp3 = smesh.CreateHypothesis("MaxElementArea", batchmode_geompy.myStudyId)
hypArea2 = hyp3._narrow(SMESH.SMESH_MaxElementArea)
hypArea2.SetMaxElementArea(500)
print hypArea2.GetName()
print hypArea2.GetId()
print hypArea2.GetMaxElementArea()
idarea2 = batchmode_smesh.AddNewHypothesis( orb.object_to_string(hypArea2) );
batchmode_smesh.SetName(idarea2, "MaxElementArea_500");
print "-------------------------- Regular_1D"
alg1 = smesh.CreateHypothesis("Regular_1D", batchmode_geompy.myStudyId)
algo1 = alg1._narrow(SMESH.SMESH_Algo)
listHyp = algo1.GetCompatibleHypothesis()
for hyp in listHyp:
print hyp
algoReg=alg1._narrow(SMESH.SMESH_Regular_1D)
print algoReg.GetName()
print algoReg.GetId()
idreg = batchmode_smesh.AddNewAlgorithms( orb.object_to_string(algoReg) );
batchmode_smesh.SetName(idreg, "Regular_1D");
print "-------------------------- MEFISTO_2D"
alg2 = smesh.CreateHypothesis("MEFISTO_2D", batchmode_geompy.myStudyId)
algo2 = alg2._narrow(SMESH.SMESH_Algo)
listHyp=algo2.GetCompatibleHypothesis()
for hyp in listHyp:
print hyp
algoMef=alg2._narrow(SMESH.SMESH_MEFISTO_2D)
print algoMef.GetName()
print algoMef.GetId()
idmef = batchmode_smesh.AddNewAlgorithms( orb.object_to_string(algoMef) );
batchmode_smesh.SetName(idmef, "MEFISTO_2D");
# ---- add hypothesis to edge
print "-------------------------- add hypothesis to edge"
edge = batchmode_geompy.IDToObject(idedge)
submesh = mesh.GetElementsOnShape(edge)
ret = mesh.AddHypothesis(edge,algoReg)
print ret
ret=mesh.AddHypothesis(edge,hypLen1)
print ret
idsm1 = batchmode_smesh.AddSubMeshOnShape( idmesh,
idedge,
orb.object_to_string(submesh),
ShapeType["Edge"] )
batchmode_smesh.SetName(idsm1, "SubMeshEdge")
batchmode_smesh.SetAlgorithms( idsm1, idreg );
batchmode_smesh.SetHypothesis( idsm1, idlength );
print "-------------------------- add hypothesis to face"
face=batchmode_geompy.IDToObject(idface)
submesh=mesh.GetElementsOnShape(face)
ret=mesh.AddHypothesis(face,hypArea2)
print ret
idsm2 = batchmode_smesh.AddSubMeshOnShape( idmesh,
idface,
orb.object_to_string(submesh),
ShapeType["Face"] )
batchmode_smesh.SetName(idsm2, "SubMeshFace")
batchmode_smesh.SetHypothesis( idsm2, idarea2 );
# ---- add hypothesis to box
print "-------------------------- add hypothesis to box"
box=batchmode_geompy.IDToObject(idbox)
submesh=mesh.GetElementsOnShape(box)
ret=mesh.AddHypothesis(box,algoReg)
print ret
ret=mesh.AddHypothesis(box,hypNbSeg1)
print ret
ret=mesh.AddHypothesis(box,algoMef)
print ret
ret=mesh.AddHypothesis(box,hypArea1)
print ret
batchmode_smesh.SetAlgorithms( idmesh, idreg );
batchmode_smesh.SetHypothesis( idmesh, idseg );
batchmode_smesh.SetAlgorithms( idmesh, idmef );
batchmode_smesh.SetHypothesis( idmesh, idarea1 );
# ---- compute box
print "-------------------------- compute box"
ret=smesh.Compute(mesh,box)
print ret
log=mesh.GetLog(0); # no erase trace
for linelog in log:
print linelog
#=======================================================
# add SuperVision component
#=======================================================
from batchmode_SuperV import *
import os
dir= os.getenv("SUPERV_ROOT_DIR")
if dir == None:
raise RuntimeError, "SUPERV_ROOT_DIR is not defined"
xmlfile = dir +"/examples/GraphEssai.xml"
print "--------------\n"+xmlfile+"\n--------------\n"
myGraph = Graph ( xmlfile )
# This DataFlow is "valid" : no loop, correct links between Nodes etc...
print myGraph.IsValid()
# Get Nodes
myGraph.PrintNodes()
Add,Sub,Mul,Div = myGraph.Nodes()
# Load Datas
Addx = Add.Input("x",3.)
Addy = Add.Input("y",4.5)
Subx = Sub.Input("x",1.5)
# Get Output Port
Addz = Add.Port('z')
Subz = Sub.Port('z')
Mulz = Mul.Port('z')
Divz = Div.Port('z')
# This DataFlow is "executable" : all pending Ports are defined with Datas
print myGraph.IsExecutable()
# Starts only execution of that DataFlow and gets control immediatly
print myGraph.Run()
# That DataFlow is running ==> 0 (false)
print myGraph.IsDone()
# Events of execution :
aStatus,aNode,anEvent,aState = myGraph.Event()
while aStatus :
print aNode.Thread(),aNode.SubGraph(),aNode.Name(),anEvent,aState
aStatus,aNode,anEvent,aState = myGraph.Event()
print myGraph.IsDone()
# Wait for Completion (but it is already done after event loop ...)
print "Done : ",myGraph.DoneW()
# Get result
print "Result : ",Divz.ToString()
# Intermediate results :
print "Intermediate Result Add\z : ",Addz.ToString()
print "Intermediate Result Sub\z : ",Subz.ToString()
print "Intermediate Result Mul\z : ",Mulz.ToString()
print " "
#print "Type : print myGraph.IsDone()"
#print " If execution is finished ==> 1 (true)"
res=myGraph.IsDone()
if res != 1:
raise RuntimeError, "myGraph.Run() is not done"
print " "
print "Type : print Divz.ToString()"
print " You will get the result"
Divz.ToString()
print " "
print "Type : myGraph.PrintPorts()"
print " to see input and output values of the graph"
myGraph.PrintPorts()
print " "
print "Type : Add.PrintPorts()"
Add.PrintPorts()
print "Type : Sub.PrintPorts()"
Sub.PrintPorts()
print "Type : Mul.PrintPorts()"
Mul.PrintPorts()
print "Type : Div.PrintPorts()"
print " to see input and output values of nodes"
Div.PrintPorts()
# Export will create newsupervisionexample.xml and the corresponding .py file
tmpdir=os.getenv("TmpDir")
if tmpdir is None:
tmpdir="/tmp"
file = tmpdir + "/newsupervisionexample"
print "--------------\n"+file+"\n--------------\n"
myGraph.Export(file)
ior = batchmode_geompy.orb.object_to_string(myGraph.G)
addStudy(ior)
GraphName = myGraph.Name()
print "Befor save ",
nodes = myGraph.Nodes()
length_bs = len(nodes)
print "ListOfNodes length = ", length_bs
names=[]
for node in nodes:
names.append(node.Name())
print names
#=================================
# save / restore study
#=================================
str= os.getenv("TmpDir")
if str == None:
str = "/tmp"
file = str+"/test.hdf"
#==================================================
#1. SaveAs
#==================================================
print " ------- We will save to", file, "-----------"
batchmode_geompy.myStudyManager.SaveAs(file, batchmode_geompy.myStudy)
batchmode_geompy.myStudyManager.Close(batchmode_geompy.myStudy)
#rename the file and try to reread it again
os.mkdir(str + "/test_dir")
#2.==================================================
print " ------- We rename the file as " + str + "/test_dir/test_new.hdf"
os.rename(file, str+"/test_dir/test_new.hdf")
#==================================================
#3. Open
#==================================================
print " ------- We try to open " + str + "/test_dir/test_new.hdf"
try:
openedStudy=batchmode_geompy.myStudyManager.Open(str+"/test_dir/test_new.hdf")
except Exception:
raise RuntimeError, "Can't open saved study!"
if openedStudy == None:
raise RuntimeError, "Can't open saved study!"
father = openedStudy.FindComponent("GEOM")
if father is None:
raise RuntimeError, "GEOM component is not found! Wrong study is opened."
#=================================
# find AttributeReal
#=================================
res,A=father.FindAttribute("AttributeReal")
if res == 0 or A == None:
raise RuntimeError, "Error: not found AttributeReal"
A = A._narrow(SALOMEDS.AttributeReal)
if A.Value() != 0.0001:
raise RuntimeError, "Error : wrong value of AttributeReal"
#==================================================
# find box
#==================================================
box = openedStudy.FindObject("box")
if box is None :
raise RuntimeError, "box was not found! Wrong study is opened."
edge=openedStudy.FindObject("face_edge")
if edge is None :
raise RuntimeError, "face_edge was not found! Wrong study is opened."
father = openedStudy.FindComponent("MESH")
if father is None:
raise RuntimeError, "MESH component is not found! Wrong study is opened."
hp=openedStudy.FindObject("Meshbox")
if hp is None:
raise RuntimeError, "Meshbox object was not found! Wrong study is opened."
hp = openedStudy.FindObject("Algorithms Definition")
if hp is None:
raise RuntimeError, "Algorithms Definition object was not found! Wrong study is opened."
hp = openedStudy.FindObject("MEFISTO_2D")
if hp is None:
raise RuntimeError, "MEFISTO_2D object was not found! Wrong study is opened."
hp = openedStudy.FindObject("Hypothesis Definition")
if hp is None:
raise RuntimeError, "Algoriths Definition was not found! Wrong study is opened."
#==================================================
# Find Supervisor
#==================================================
father = openedStudy.FindComponent("SUPERV")
if father is None:
raise RuntimeError, "SUPERV component is not found! Wrong study is opened."
SuperV = lcc.FindOrLoadComponent("SuperVisionContainer","Supervision")
Builder = openedStudy.NewBuilder()
Builder.LoadWith(father, SuperV)
import SALOMEDS
aChildIterator = openedStudy.NewChildIterator(father)
#while aChildIterator.More():
anSObject = aChildIterator.Value()
#print "iterate: ", anSObject.GetID()
res, anAttr=anSObject.FindAttribute("AttributeIOR")
if res :
anAttr=anAttr._narrow(SALOMEDS.AttributeIOR)
ior = anAttr.Value()
Graph=SuperV.getGraph(ior)
ListOfNodes=Graph.Nodes()
length_as= len(ListOfNodes)
print "ListOfNodes length = ", length_as
if length_as != length_bs:
raise RuntimeErrror, "defferent length of nodes after study open"
#aChildIterator.Next()
Names = []
for node in ListOfNodes:
Names.append(node.Name())
names.sort()
Names.sort()
if names != Names :
raise RuntimeError, "List of dataflow nodes after save differs from one befor save operation"
#==================================================
#4. Save
#==================================================
batchmode_geompy.myStudyManager.Save(openedStudy)
batchmode_geompy.myStudyManager.Close(openedStudy)
#==================================================
#5. Open
#==================================================
try:
openedStudy = batchmode_geompy.myStudyManager.Open(str+"/test_dir/test_new.hdf")
except Exception:
raise RuntimeError, "Can't open saved study!"
if openedStudy == None:
raise RuntimeError, "Can't open saved study!"
father = openedStudy.FindComponent("GEOM")
if father is None:
raise RuntimeError, "Geom component is not found! Wrong study is opened."
#find AttributeReal
res,A=father.FindAttribute("AttributeReal")
if res == 0 or A == None:
raise RuntimeError, "Error: not found AttributeReal"
A = A._narrow(SALOMEDS.AttributeReal)
if A.Value() != 0.0001:
raise RuntimeError, "Error : wrong value of AttributeReal"
#==================================================
# find box
#==================================================
box = openedStudy.FindObject("box")
if box is None :
raise RuntimeError, "box was not found! Wrong study is opened."
edge=openedStudy.FindObject("face_edge")
if edge is None :
raise RuntimeError, "face_edge was not found! Wrong study is opened."
father = openedStudy.FindComponent("MESH")
if father is None:
raise RuntimeError, "MESH component is not found! Wrong study is opened."
hp=openedStudy.FindObject("Meshbox")
if hp is None:
raise RuntimeError, "Meshbox object was not found! Wrong study is opened."
hp = openedStudy.FindObject("Algorithms Definition")
if hp is None:
raise RuntimeError, "Algorithms Definition object was not found! Wrong study is opened."
hp = openedStudy.FindObject("MEFISTO_2D")
if hp is None:
raise RuntimeError, "MEFISTO_2D object was not found! Wrong study is opened."
hp = openedStudy.FindObject("Hypothesis Definition")
if hp is None:
raise RuntimeError, "Algoriths Definition was not found! Wrong study is opened."
#==================================================
# Find Supervisor
#==================================================
father = openedStudy.FindComponent("SUPERV")
if father is None:
raise RuntimeError, "SUPERV component is not found! Wrong study is opened."
SuperV = lcc.FindOrLoadComponent("SuperVisionContainer","Supervision")
Builder = openedStudy.NewBuilder()
Builder.LoadWith(father, SuperV)
import SALOMEDS
aChildIterator = openedStudy.NewChildIterator(father)
#while aChildIterator.More():
anSObject = aChildIterator.Value()
res, anAttr=anSObject.FindAttribute("AttributeIOR")
if res :
anAttr=anAttr._narrow(SALOMEDS.AttributeIOR)
ior = anAttr.Value()
Graph=SuperV.getGraph(ior)
ListOfNodes=Graph.Nodes()
length_as= len(ListOfNodes)
print "ListOfNodes length = ", length_as
if length_as != length_bs:
raise RuntimeErrror, "defferent length of nodes after study open"
#aChildIterator.Next()
Names = []
for node in ListOfNodes:
Names.append(node.Name())
names.sort()
Names.sort()
if names != Names :
raise RuntimeError, "List of dataflow nodes after save differs from one befor save operation"