From 5283ba55468a71c407e783cf8a140d642056dbd5 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Mon, 31 May 2010 07:12:45 +0000 Subject: [PATCH] if it is pReal it should be 1.0 not just 1 --- code/math.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/code/math.f90 b/code/math.f90 index 49a3eaeb2..abb0ab759 100644 --- a/code/math.f90 +++ b/code/math.f90 @@ -1305,10 +1305,10 @@ pure function math_transpose3x3(A) math_RtoQuaternion = 0.0_pReal - absQ(1) = 0.5_pReal * dsqrt(1_pReal+R(1,1)+R(2,2)+R(3,3)) - absQ(2) = 0.5_pReal * dsqrt(1_pReal+R(1,1)-R(2,2)-R(3,3)) - absQ(3) = 0.5_pReal * dsqrt(1_pReal-R(1,1)+R(2,2)-R(3,3)) - absQ(4) = 0.5_pReal * dsqrt(1_pReal-R(1,1)-R(2,2)+R(3,3)) + absQ(1) = 0.5_pReal * dsqrt(1.0_pReal+R(1,1)+R(2,2)+R(3,3)) + absQ(2) = 0.5_pReal * dsqrt(1.0_pReal+R(1,1)-R(2,2)-R(3,3)) + absQ(3) = 0.5_pReal * dsqrt(1.0_pReal-R(1,1)+R(2,2)-R(3,3)) + absQ(4) = 0.5_pReal * dsqrt(1.0_pReal-R(1,1)-R(2,2)+R(3,3)) largest = maxloc(absQ) select case(largest(1))