如何解决如何将以下 Mathematica 代码转换为 sympy?
我正在尝试将一些 Mathematica 代码转换为 python,但我不知道如何处理以下函数。有没有办法使用 sympy 或其他 python 模块? TPUV、AM1UV和GUV是由其他伴随函数计算的向量。
perpUnitVector[]:=Module[
{sol,a3,b3,c3,a4,b4,c4},sol=Quiet[Solve[
{a3,c3}.TPUV==0&&
{a3,c3}.AM1UV==0&&
{a3,c3}.{a3,c3}==1&&
Sign[Cross[-GUV,{a3,c3}].Cross[{1,0},-GUV]]>0,c3},Reals]];
{a4,c4}=Flatten[{a3,c3}/.sol[[1]],1]
]
测试用例: 当 TPU = {1,3,4},AM1UV = {4,5,6},GUV = {1,1,1} 时,Mathematica 给出 a4 = -0.408248,b4 = 0.816497,c4 = -0.408248。
我的同情尝试在这里但似乎不起作用
import numpy as np
import sympy as sym
GUV = [1.0,1.0,1.0];
AM1UV = [4.0,5.0,6.0];
TPUV = [1.0,2.0,3.0];
a3,c4 = sym.symbols('a3 b3 c3 a4 b4 c4')
sol = sym.solveset(np.dot([a3,c3],TPUV),np.dot([a3,AM1UV),[a3,c3])-1,np.sign(np.dot(np.cross(-GUV,c3]),np.cross([1,0],-GUV)))>0,S.Reals)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。