主题：fortran中 叉乘运算？的编程？
apollo1028
[专家分：0] 发布于 2011-09-19 20:14:00
求 fortran中 叉乘运算？的编程 好方法
回复列表 （共3个回复）
沙发
cgl_lgs [专家分：21040] 发布于 2011-09-19 21:05:00
您是说MatMut么？又或是dgemm？
板凳
yeg001 [专家分：14390] 发布于 2011-09-19 21:47:00
楼主可能是讲向量叉乘. 内部函数记得是没的, 要自己写. 按照数学计算方法写就是了. 要重载到运算符那可能就麻烦一点.
3 楼
adda [专家分：1520] 发布于 2011-09-20 12:54:00
function cross_product( vect1, vect2 ) result( vectout )
implicit none
real, dimension(3), intent(in) :: vect1, vect2
real, dimension(3) :: vectout
vectout(1) = vect1(2) * vect2(3) - vect1(3) * vect2(2)
vectout(2) = vect1(3) * vect2(1) - vect1(1) * vect2(3)
vectout(3) = vect1(1) * vect2(2) - vect1(2) * vect2(1)
end function cross_product
我来回复
您所在位置：社区首页 — Fortran讨论区 — fortran中 叉乘运算？的编程？