natural Fortran order
This commit is contained in:
parent
72ba4645cb
commit
5267aff1a2
|
@ -737,13 +737,13 @@ pure function eu2om(eu) result(om)
|
||||||
s = sin(eu)
|
s = sin(eu)
|
||||||
|
|
||||||
om(1,1) = c(1)*c(3)-s(1)*s(3)*c(2)
|
om(1,1) = c(1)*c(3)-s(1)*s(3)*c(2)
|
||||||
om(1,2) = s(1)*c(3)+c(1)*s(3)*c(2)
|
|
||||||
om(1,3) = s(3)*s(2)
|
|
||||||
om(2,1) = -c(1)*s(3)-s(1)*c(3)*c(2)
|
om(2,1) = -c(1)*s(3)-s(1)*c(3)*c(2)
|
||||||
om(2,2) = -s(1)*s(3)+c(1)*c(3)*c(2)
|
|
||||||
om(2,3) = c(3)*s(2)
|
|
||||||
om(3,1) = s(1)*s(2)
|
om(3,1) = s(1)*s(2)
|
||||||
|
om(1,2) = s(1)*c(3)+c(1)*s(3)*c(2)
|
||||||
|
om(2,2) = -s(1)*s(3)+c(1)*c(3)*c(2)
|
||||||
om(3,2) = -c(1)*s(2)
|
om(3,2) = -c(1)*s(2)
|
||||||
|
om(1,3) = s(3)*s(2)
|
||||||
|
om(2,3) = c(3)*s(2)
|
||||||
om(3,3) = c(2)
|
om(3,3) = c(2)
|
||||||
|
|
||||||
where(abs(om)<1.0e-12_pReal) om = 0.0_pReal
|
where(abs(om)<1.0e-12_pReal) om = 0.0_pReal
|
||||||
|
|
Loading…
Reference in New Issue