sympy

[Solved] Do the projection (with Jacobian) and marginalisation (inversion of matrix and remove a row/column and reinversion) commute?

by

To do np.dot last dimension of first matrix must be the same as first dimension of second one. They are not, so you are getting ValueError, that shapes are not aligned. Everything seems to be fine as you printed, but then you forgot about lines: j_temp = np.copy(J_2_SYM) # Add row/col into J_2_SYM j_temp = … Read more

[Solved] triangular pyramid, trigonal pyramid offset-3d-Plane-Intersection By vbasolver

by

I check from sympy import * var(‘m’) def myUnitVector(myPoint3D): myL=myPoint3D.distance((0, 0)) return Point3D(myPoint3D.x/myL,myPoint3D.y/myL,myPoint3D.z/myL) def myHtoP(myHairetu): return Point3D(myHairetu[0],myHairetu[1],myHairetu[2]) def my3PLaneIntersection(PTO,PTA,PTB,PTC,RA,RB,RC): vA =myUnitVector(myHtoP(Plane(PTO, PTB, PTC).normal_vector)) PLA=Plane(PTO + RA * vA, normal_vector=vA) vB =myUnitVector(myHtoP(Plane(PTO, PTC, PTA).normal_vector)) PLB=Plane(PTO + RB * vB, normal_vector=vB) vC =myUnitVector(myHtoP(Plane(PTO, PTA, PTB).normal_vector)) PLC=Plane(PTO + RC * vC, normal_vector=vC) return PLC.intersection(PLB.intersection(PLA)[0]) def myProjection(PTO,PTA,PTC,PTOO): v = … Read more