+ j*a^RIHtHtHtHtHtHtHtHt H t ^RI H t ^RI Ht^RIHt^RIHt^RIHtR.t!RR]]4t!R R ]4tR tR #) ) SsympifyexpandsqrtAddzerosacosImmutableMatrixsimplify)trigsimp) Printable) filldedent) EvalfMixin) prec_to_dpsVectorca]tRt^toRtRtRtRt]R4t Rt Rt Rt Rt R tR tR tR tR tRtRtR(RltRtRt] t]tRtRt] P]n]tRt]P]nRt]P]nR(RltR)Rlt Rt!Rt"Rt#Rt$Rt%Rt&Rt'R t(R!t)R"t*R#t+R$t,R%t-R&t.Vt/R'#)*raThe class used to define vectors. It along with ReferenceFrame are the building blocks of describing a classical mechanics system in PyDy and sympy.physics.vector. Attributes ========== simp : Boolean Let certain methods use trigsimp on their outputs Fc .VnV^8Xd.p\V\4'dTpMO/pVFGpV^,V9d%W#^,;;,V^,, uu&K5V^,W#^,&KI VP4F5wrEV\ .RO48wgKVPP WT34K7 R#)a_This is the constructor for the Vector class. You should not be calling this, it should only be used by other functions. You should be treating Vectors like you would with if you were doing the math by hand, and getting the first 3 from the standard basis vectors from a ReferenceFrame. The only exception is to create a zero vector: zv = Vector(0) N)rrr)args isinstancedictitemsMatrixappend)selfinlistdinpkvs&& \/home/luca/screen-checking/.venv/lib/python3.14/site-packages/sympy/physics/vector/vector.py__init__Vector.__init__s Q;F fd # #AAq6Q;!fIQ'I #AA!fI  GGIDAF9%%   !(c \#)zReturns the class Vector. )rrs&r func Vector.func<s  r#c>\\VP44#N)hashtuplerr%s&r __hash__Vector.__hash__AsE$))$%%r#c vV^8XdV#\V4p\VPVP,4#)zThe add operator for Vector. ) _check_vectorrrrothers&&r __add__Vector.__add__Ds0 A:Ke$dii%**,--r#c h^RIHpHp\W4'dcV!V4p\ ^4pVP F=pWE^,V^,,V^,P V4,, pK? V#\V4p\PpVP FhpVP FUpWh^,PV^,PV^,4,V^,,^,, pKW Kj \P'd\VRR7#V#)aDot product of two vectors. Returns a scalar, the dot product of the two Vectors Parameters ========== other : Vector The Vector which we are dotting with Examples ======== >>> from sympy.physics.vector import ReferenceFrame, dot >>> from sympy import symbols >>> q1 = symbols('q1') >>> N = ReferenceFrame('N') >>> dot(N.x, N.x) 1 >>> dot(N.x, N.y) 0 >>> A = N.orientnew('A', 'Axis', [q1, N.x]) >>> dot(N.y, A.y) cos(q1) Dyadic _check_dyadicT) recursive)sympy.physics.vector.dyadicr6r7rrrdotr/rZeroTdcmsimpr ) rr1r6r7olroutv1v2s && r r: Vector.dotKs8 F e $ $!%(EBZZdQqTkQqTXXd^44 Ie$ff))BjjARUYYr!u%56"Q%@!DD! ;;;C40 0Jr#c NVP\PV, 4#)z6This uses mul and inputs self and 1 divided by other. )__mul__rOner0s&&r __truediv__Vector.__truediv__xs||AEEEM**r#c V^8Xd \^4p\V4pTP.8XdTP.8XdR#TP.8XgTP.8XdR#TP^,^,pTF+p\ Y, P T44^8wgK*R# R# \dR#i;i)a<Tests for equality. It is very import to note that this is only as good as the SymPy equality test; False does not always mean they are not equivalent Vectors. If other is 0, and self is empty, returns True. If other is 0 and self is not empty, returns False. If none of the above, only accepts other as a Vector. FT)rr/ TypeErrorrrr:)rr1framers&& r __eq__ Vector.__eq__|s A:1IE !%(E IIO%**"2ii2o5::#3 ! QAt|((+,1  s B44 CCc \VP4p\V4p\\ V44F*pWV,^,,W#,^,3W#&K, \ V4#)aLMultiplies the Vector by a sympifyable expression. Parameters ========== other : Sympifyable The scalar to multiply this Vector with Examples ======== >>> from sympy.physics.vector import ReferenceFrame >>> from sympy import Symbol >>> N = ReferenceFrame('N') >>> b = Symbol('b') >>> V = 10 * b * N.x >>> print(V) 10*b*N.x )listrrrangelenr)rr1newlistis&& r rEVector.__mul__sT,tyy/s7|$A!*Q-/A?GJ%gr#cVR,#r%s&r __neg__Vector.__neg__s byr#c  ^RIHp\V4pV!^4pVPEF(pVPEFpW2!V^,^,V^,^,,V^,PV^,P3.4, pW2!V^,^,V^,^,,V^,PV^,P 3.4, pW2!V^,^,V^,^,,V^,PV^,P 3.4, pW2!V^,^,V^,^,,V^,P V^,P3.4, pW2!V^,^,V^,^,,V^,P V^,P 3.4, pW2!V^,^,V^,^,,V^,P V^,P 3.4, pW2!V^,^,V^,^,,V^,P V^,P3.4, pW2!V^,^,V^,^,,V^,P V^,P 3.4, pW2!V^,^,V^,^,,V^,P V^,P 3.4, pEK EK+ V#)aHOuter product between two Vectors. A rank increasing operation, which returns a Dyadic from two Vectors Parameters ========== other : Vector The Vector to take the outer product with Examples ======== >>> from sympy.physics.vector import ReferenceFrame, outer >>> N = ReferenceFrame('N') >>> outer(N.x, N.x) (N.x|N.x) )r6)r9r6r/rxyz)rr1r6r?rrBs&& r outer Vector.outersm* 7e$ AYAjjfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEEfqtAwAq11Q4662a577CDEE! r#c VPp\V4^8Xd \^4#.pVEFRpR EFGpV^,V,^8Xd3VPRV^,PV,,4KKV^,V,R 8Xd3VPRV^,PV,,4KV^,V,^8wgKVP V^,V,4p\ V^,V,\4'd RV,pV^,R8Xd VR,pRpMRpVPWv,V^,PV,,4EKJ EKU RPV4pVPR4'd VR,pV#VPR4'd VR,pV#) zLatex Printing method.  +  - (%s)-rWNNNN rrWrX) rrQstrr latex_vecs_printrrjoin startswith) rprinterarr?rjarg_str str_startoutstrs && r _latex Vector._latexskYY r7a<q6M AQ47a<IIeadooa&889qT!W]IIeadooa&889qT!W\&nnQqT!W5G!!A$q'3//"(7"2qzS(")"+$) $) IIi1AaDOOA4FFG%(   U # #BZF   s # #BZF r#c  ,aa^RIHo.pVV3RlpVPFwrE\^4FpWF,^8XdKWF,^8Xd+VP S!VP V,44KIWF,R8Xd8VP S!R4S!VP V,4,4KVP V!WF,VP V,44K K V'dSP !V!pV#S!R4pV#)zPretty Printing method. ) prettyFormc<SPV4pSPV4pVP'dS!VP4!pSPW#.RR7#)rk) delimiter)rpis_Addparens _print_seq)abpapbr|rss&& r juxtapose!Vector._pretty..juxtaposesP"B"Bxxx-%%rh#%> >r#z-10rX) sympy.printing.pretty.stringpictr|rrPr pretty_vecsr2) rrstermsrMNrS pretty_resultr|s &f @r _prettyVector._prettys? ?IIDA1X419TQYLLAMM!,TRZLLD!1Jq}}Q?O4P!PQLL14q1A!BC &..6M'sOMr#c"RV,V,#rVrYr0s&&r __rsub__Vector.__rsub__!sT U""r#c V'd\VP4^8Xd\VP4pM\VP4^8XdVP^4#VPUu/uFqD^,V^,bK pp\ VP 4RR7p.pVFpVP WW,V34K .pVEFYpR EFNp V^,V ,^8Xd3VP RV^,PV ,,4KKV^,V ,R 8Xd3VP RV^,PV ,,4KV^,V ,^8wgKVPV^,V ,4p \V^,V ,\4'd RV ,p V ^,R8Xd V R,p Rp MRp VP W,R,V^,PV ,,4EKQ EK\ R PV4p V PR4'd V R ,p V #V PR 4'd V R,p V #uupi)zPrinting method. cVP#r))index)r]s&r "Vector._sympystr..,s!''r#)keyrcrdrerfrg*rhrirkrlrX) rQrrOrpsortedkeysrstr_vecsrrrqrr) rrsorderrtrrrrr?rurvrwrxs &&& r _sympystrVector._sympystr$sDII!+diiB ^q >>!$ $%)YY/Y1qtYA/!&&((9:DB 163-( AQ47a<IIeadmmA&667qT!W]IIeadmmA&667qT!W\&nnQqT!W5G!!A$q'3//"(7"2qzS(")"+$) $) IIi1C7!A$--:JJK%(   U # #BZF   s # #BZF ?0s2Jc 2VPVR,4#)zThe subtraction operator. rX)r2r0s&&r __sub__Vector.__sub__Ls||EBJ''r#c \^RIHpHp\W4'dsV!V4pV!^4p\ VP 4FGwrVWF^,VP V^,4PV^,4,, pKI V#\V4pVP .8Xd \^4#Rp.pVP p V FpV^,Pp V^,Pp V^,Pp WV .VPV 4VPV 4VPV 4.\V.4PV 4\V.4PV 4\V.4PV 4..p W!V 4P , pK \V4#)aThe cross product operator for two Vectors. Returns a Vector, expressed in the same ReferenceFrames as self. Parameters ========== other : Vector The Vector which we are crossing with Examples ======== >>> from sympy import symbols >>> from sympy.physics.vector import ReferenceFrame, cross >>> q1 = symbols('q1') >>> N = ReferenceFrame('N') >>> cross(N.x, N.y) N.z >>> A = ReferenceFrame('A') >>> A.orient_axis(N, q1, N.x) >>> cross(A.x, N.y) N.z >>> cross(N.y, A.x) - sin(q1)*A.y - cos(q1)*A.z r5cnV^,^,V^,^,V^,^,,V^,^,V^,^,,, ,V^,^,V^,^,V^,^,,V^,^,V^,^,,, ,,V^,^,V^,^,V^,^,,V^,^,V^,^,,, ,,#)zThis is needed as a little method for to find the determinant of a list in python; needs to work for a 3x3 list. SymPy's Matrix will not take in Vector, so need a custom function. You should not be calling this. rY)mats&r _detVector.cross.._detxsF1IQSVAY!6QSVAY9N!NO!fQi3q6!9s1vay#83q6!9F1I<$!$Qc!fQi#a&).CF1IAq )/*"++ ,r#)r9r6r7r enumeratercrossr`r/rr]r^r_r:)rr1r6r7r?rSrroutlistrttempxtempytempztempms&& r r Vector.crossPs`: F e $ $!%(EB!%**-dtzz!A$/66qt<==.Ie$ :: !9  , ZZAaDFFEaDFFEaDFFEU+xx%Iqc{u-vqc{u/Eqc{u-/0E tE{'' 'Ggr#c Z/pVPFp\V.4W^,&K V#)as The constituents of this vector in different reference frames, as per its definition. Returns a dict mapping each ReferenceFrame to the corresponding constituent Vector. Examples ======== >>> from sympy.physics.vector import ReferenceFrame >>> R1 = ReferenceFrame('R1') >>> R2 = ReferenceFrame('R2') >>> v = R1.x + R2.x >>> v.separate() == {R1: R1.x, R2: R2.x} True )rr)r componentsr]s& r separateVector.separates/( A%qc{Jt r#c$VPV4#r))r:r0s&&r __and__Vector.__and__sxxr#c$VPV4#r))rr0s&&r __xor__Vector.__xor__zz%  r#c$VPV4#r))r`r0s&&r __or__ Vector.__or__rr#c F^RIHpV!V4\V4p.pVPFpV^,pV^,pW8XdWWP V4V3., pK7V'd0VP V4P V4\ ^^48XdWWP V4V3., pK\V.4PV4p V P^,^,P V4p V\W3.4P, pK \V4#)aReturns the partial derivative of the vector with respect to a variable in the provided reference frame. Parameters ========== var : Symbol What the partial derivative is taken with respect to. frame : ReferenceFrame The reference frame that the partial derivative is taken in. var_in_dcm : boolean If true, the differentiation algorithm assumes that the variable may be present in any of the direction cosine matrices that relate the frame to the frames of any component of the vector. But if it is known that the variable is not present in the direction cosine matrices, false can be set to skip full reexpression in the desired frame. Examples ======== >>> from sympy import Symbol >>> from sympy.physics.vector import dynamicsymbols, ReferenceFrame >>> from sympy.physics.vector import init_vprinting >>> init_vprinting(pretty_print=False) >>> t = Symbol('t') >>> q1 = dynamicsymbols('q1') >>> N = ReferenceFrame('N') >>> A = N.orientnew('A', 'Axis', [q1, N.y]) >>> A.x.diff(t, N) - sin(q1)*q1'*N.x - cos(q1)*q1'*N.z >>> A.x.diff(t, N).express(A).simplify() - q1'*A.z >>> B = ReferenceFrame('B') >>> u1, u2 = dynamicsymbols('u1, u2') >>> v = u1 * A.x + u2 * B.y >>> v.diff(u2, N, var_in_dcm=False) B.y ) _check_frame) sympy.physics.vector.framerrrdiffr=rrexpress) rvarrK var_in_dcmrrvector_componentmeasure_numbercomponent_framereexp_vec_compderivs &&&& r r Vector.diffsR <Ucl $ -a0N.q1O'//4e<=="eii&@&E&Ec&J&+Aqk'2 3 3C 8/JKKF%+-=,>%?%G%G%NN*//215::3?Efun%56;;;F!* f~r#c "^RIHpV!WVR7#)a Returns a Vector equivalent to this one, expressed in otherframe. Uses the global express method. Parameters ========== otherframe : ReferenceFrame The frame for this Vector to be described in variables : boolean If True, the coordinate symbols(if present) in this Vector are re-expressed in terms otherframe Examples ======== >>> from sympy.physics.vector import ReferenceFrame, dynamicsymbols >>> from sympy.physics.vector import init_vprinting >>> init_vprinting(pretty_print=False) >>> q1 = dynamicsymbols('q1') >>> N = ReferenceFrame('N') >>> A = N.orientnew('A', 'Axis', [q1, N.y]) >>> A.x.express(N) cos(q1)*N.x - sin(q1)*N.z )r) variables)sympy.physics.vectorr)r otherframerrs&&& r rVector.expresss8 1t9==r#c ~\VUu.uFq PV4NK up4P^^4#uupi)aReturns the matrix form of the vector with respect to the given frame. Parameters ---------- reference_frame : ReferenceFrame The reference frame that the rows of the matrix correspond to. Returns ------- matrix : ImmutableMatrix, shape(3,1) The matrix that gives the 1D vector. Examples ======== >>> from sympy import symbols >>> from sympy.physics.vector import ReferenceFrame >>> a, b, c = symbols('a, b, c') >>> N = ReferenceFrame('N') >>> vector = a * N.x + b * N.y + c * N.z >>> vector.to_matrix(N) Matrix([ [a], [b], [c]]) >>> beta = symbols('beta') >>> A = N.orientnew('A', 'Axis', (beta, N.x)) >>> vector.to_matrix(A) Matrix([ [ a], [ b*cos(beta) + c*sin(beta)], [-b*sin(beta) + c*cos(beta)]]) )rr:reshape)rreference_frameunit_vecs&& r to_matrixVector.to_matrixsDJ&(&/7xx)&())0A 7(s:c a/pVPF(pV^,PV3Rl4W#^,&K* \V4#)z(Calls .doit() on each term in the Vectorc(<VP!R/SB#)NrY)doit)r]hintss&r rVector.doit..Hsqvvr#)r applyfuncr)rrrrs&l r r Vector.doitDs9 Adnn%>?AdGayr#c ^RIHpV!W4#)z Returns a Vector which is the time derivative of the self Vector, taken in frame otherframe. Calls the global time_derivative method Parameters ========== otherframe : ReferenceFrame The frame to calculate the time derivative in )time_derivative)rr)rrrs&& r dt Vector.dtKs 9t00r#c x/pVPFp\V^,4W^,&K \V4#)zReturns a simplified Vector.)rr r)rrrs& r r Vector.simplify\s2 AqtnAdGayr#c /pVPF%pV^,P!V/VBW4^,&K' \V4#)zSubstitution on the Vector. Examples ======== >>> from sympy.physics.vector import ReferenceFrame >>> from sympy import Symbol >>> N = ReferenceFrame('N') >>> s = Symbol('s') >>> a = N.x * s >>> a.subs({s: 2}) 2*N.x )rsubsr)rrkwargsrrs&*, r r Vector.subscs> Adii00AdGayr#c 6\VPV44#)zReturns the magnitude (Euclidean norm) of self. Warnings ======== Python ignores the leading negative sign so that might give wrong results. ``-A.x.magnitude()`` would be treated as ``-(A.x.magnitude())``, instead of ``(-A.x).magnitude()``. )rr:r%s&r magnitudeVector.magnitudexsDHHTN##r#c d\VP.,4VP4, #)z9Returns a Vector of magnitude 1, codirectional with self.)rrrr%s&r normalizeVector.normalizes!dii"n%(888r#c \V4'g \R4h/pVPF$pV^,PV4W#^,&K& \ V4#)z/Apply a function to each component of a vector.z`f` must be callable.)callablerJrrr)rfrrs&& r rVector.applyfuncsL{{34 4 AdnnQ'AdGayr#c zVP4pVP4p\VPV44pV#)a; Returns the smallest angle between Vector 'vec' and self. Parameter ========= vec : Vector The Vector between which angle is needed. Examples ======== >>> from sympy.physics.vector import ReferenceFrame >>> A = ReferenceFrame("A") >>> v1 = A.x >>> v2 = A.y >>> v1.angle_between(v2) pi/2 >>> v3 = A.x + A.y + A.z >>> v1.angle_between(v3) acos(sqrt(3)/3) Warnings ======== Python ignores the leading negative sign so that might give wrong results. ``-A.x.angle_between()`` would be treated as ``-(A.x.angle_between())``, instead of ``(-A.x).angle_between()``. )rr r:)rvecvec1vec2angles&& r angle_betweenVector.angle_betweens3B~~}}TXXd^$ r#c 8VPV4P#)anReturns the free symbols in the measure numbers of the vector expressed in the given reference frame. Parameters ========== reference_frame : ReferenceFrame The frame with respect to which the free symbols of the given vector is to be determined. Returns ======= set of Symbol set of symbols present in the measure numbers of ``reference_frame``. )r free_symbols)rrs&&r rVector.free_symbolss$~~o.;;;r#c ^RIHpV!WR7#)aReturns the free dynamic symbols (functions of time ``t``) in the measure numbers of the vector expressed in the given reference frame. Parameters ========== reference_frame : ReferenceFrame The frame with respect to which the free dynamic symbols of the given vector is to be determined. Returns ======= set Set of functions of time ``t``, e.g. ``Function('f')(me.dynamicsymbols._t)``. )find_dynamicsymbols)r)!sympy.physics.mechanics.functionsr)rrrs&& r free_dynamicsymbolsVector.free_dynamicsymbolss& J"4IIr#cVP'gV#.p\V4pVPF(wrEVPVPVR7V.4K* \ V4#))n)rrrevalfr)rprecnew_argsdpsrrKs&& r _eval_evalfVector._eval_evalfsVyyyK$))JC OOSYYY-u5 6$hr#c .pVPF(wr4VPV4pVPW4.4K* \V4#)aReplace occurrences of objects within the measure numbers of the vector. Parameters ========== rule : dict-like Expresses a replacement rule. Returns ======= Vector Result of the replacement. Examples ======== >>> from sympy import symbols, pi >>> from sympy.physics.vector import ReferenceFrame >>> A = ReferenceFrame('A') >>> x, y, z = symbols('x y z') >>> ((1 + x*y) * A.x).xreplace({x: pi}) (pi*y + 1)*A.x >>> ((1 + x*y) * A.x).xreplace({x: pi, y: 2}) (1 + 2*pi)*A.x Replacements occur only if an entire node in the expression tree is matched: >>> ((x*y + z) * A.x).xreplace({x*y: pi}) (z + pi)*A.x >>> ((x*y*z) * A.x).xreplace({x*y: pi}) x*y*z*A.x )rxreplacerr)rrulerrrKs&& r r Vector.xreplacesCL))JC,,t$C OOSL )$hr#)rN)T)F)0__name__ __module__ __qualname____firstlineno____doc__r> is_numberr!propertyr&r,r2r:rGrLrErZr`ryrrrrr__radd____rmul__rr__rand__rrrrrrrr rrrrrrrrr __static_attributes____classdictcell__) __classdict__s@r rrs  DI):&.+Z+:8&P D>#&P(@DHH2kkGOH!mmGO!]]FN@D>>&7P1"* $9$L<(J. * * r#c2aa]tRtRtoV3RltRtVtV;t#)VectorTypeErroric z<\R\V4: RV: R\V4: R24p\SV` V4R#)zExpected an instance of z, but received object 'z' of .N)rtypesuperr!)rr1wantmsg __class__s&&& r r!VectorTypeError.__init__s0*.t*eT%[JK r#rY)r rrrr!rr __classcell__)r"rs@@r rrsr#rcH\V\4'g \R4hV#)zA Vector must be supplied)rrrJ)r1s&r r/r/#s eV $ $344 Lr#N)sympyrrrrrrr r rr sympy.simplify.trigsimpr sympy.printing.defaultsr sympy.utilities.miscrsympy.core.evalfrmpmath.libmp.libmpfr__all__rrJrr/rYr#r r-sN888,-+'+ *J  Y J  Zir#