__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.')
