__version__ = "3.0" __email__ = "*****@*****.**" feProblem = xc.FEProblem() preprocessor = feProblem.getPreprocessor # Defining nodes points = preprocessor.getCad.getPoints pt1 = points.newPntIDPos3d(1, geom.Pos3d(1.0, 0.0, 0.0)) pt2 = points.newPntIDPos3d(2, geom.Pos3d(0.0, 0.0, 0.0)) eje = geom.Recta3d(pt1.getPos, pt2.getPos) ang = math.pi / 2.0 trfs = preprocessor.getCad.getGeometricTransformations rot1 = trfs.newTransformation("rotation") rot1.setAxisAng(eje, ang) pt3 = points.newPntIDPos3d(3, rot1.getTrfPos(geom.Pos3d(0.0, 0.0, 1.0))) pt4 = points.newPntIDPos3d(4, geom.Pos3d(0.0, 1.0, 0.0)) dist34 = mi.distance_2points(preprocessor, 3, 4) # print "dist(3,4)= ",dist34 import os from miscUtils import LogMessages as lmsg fname = os.path.basename(__file__) if (dist34 < 1e-12): print "test ", fname, ": ok." else: lmsg.error(fname + ' ERROR.')
__email__= "*****@*****.**" feProblem= xc.FEProblem() preprocessor= feProblem.getPreprocessor # Defining nodes points= preprocessor.getMultiBlockTopology.getPoints pt1= points.newPntIDPos3d(1,geom.Pos3d(0,0,0)) pt2= points.newPntIDPos3d(2,geom.Pos3d(0,0,1)) pt3= points.newPntIDPos3d(3,geom.Pos3d(1,0,0)) axis= geom.Line3d(pt1.getPos,pt2.getPos) ang= math.pi/2.0 trfs= preprocessor.getMultiBlockTopology.getGeometricTransformations rot1= trfs.newTransformation("rotation") rot1.setAxisAng(axis,-ang) pt4= points.newPntIDPos3d(4,rot1.getTrfPos(geom.Pos3d(0.0, 1.0, 0.0))) dist34= mi.distance_2points(preprocessor, 3,4) #print "dist(3,4)= ",dist34 import os from miscUtils import LogMessages as lmsg fname= os.path.basename(__file__) if (dist34<1e-12): print "test ",fname,": ok." else: lmsg.error(fname+' ERROR.')
pts = preprocessor.getMultiBlockTopology.getPoints ptA = pts.newPntIDPos3d(1, geom.Pos3d(0.0, 0.0, 3.0)) ptB = pts.newPntIDPos3d(2, geom.Pos3d(0.0, 0.0, 4.0)) trfs = preprocessor.getMultiBlockTopology.getGeometricTransformations transl = trfs.newTransformation("translation") v = ptB.getPos - ptA.getPos transl.setVector(v) ptTrfA = pts.newPntIDPos3d(10, transl.getTrfPos(ptA.getPos)) ptTrfB = pts.newPntIDPos3d(11, transl.getTrfPos(ptB.getPos)) ''' for p in pts: pos= p.getPos print "id: ",p.getNombre," x= ",pos.x,," y= ",pos.y,," z= ",pos.z ''' dist110 = mi.distance_2points(preprocessor, 1, 10) # print "dist(1,10)= ",dist110 dist211 = mi.distance_2points(preprocessor, 2, 11) # print "dist(2,11)= ",dist211 import os from misc_utils import log_messages as lmsg fname = os.path.basename(__file__) if (abs(dist110 - 1.0) < 1e-12) & (abs(dist211 - 1.0) < 1e-12): print "test ", fname, ": ok." else: lmsg.error(fname + ' ERROR.')
trfs= preprocessor.getMultiBlockTopology.getGeometricTransformations transl= trfs.newTransformation("translation") v= ptB.getPos-ptA.getPos transl.setVector(v) ptTrfA= pts.newPntIDPos3d(10,transl.getTrfPos(ptA.getPos)) ptTrfB= pts.newPntIDPos3d(11,transl.getTrfPos(ptB.getPos)) ''' for p in pts: pos= p.getPos print "id: ",p.getNombre," x= ",pos.x,," y= ",pos.y,," z= ",pos.z ''' dist110= mi.distance_2points(preprocessor, 1,10) # print "dist(1,10)= ",dist110 dist211= mi.distance_2points(preprocessor, 2,11) # print "dist(2,11)= ",dist211 import os from miscUtils import LogMessages as lmsg fname= os.path.basename(__file__) if (abs(dist110-1.0)<1e-12) & (abs(dist211-1.0)<1e-12): print "test ",fname,": ok." else: lmsg.error(fname+' ERROR.')