Commit 0f65821d authored by Laurent Belcour's avatar Laurent Belcour

[Debug] Correction of the NEUMANN_3D transformation. Need to accound for

below hemisphere configuration in the data2data code.
parent 6e4e9461
......@@ -279,8 +279,6 @@ void params::to_cartesian(const double* invec, params::input intype,
outvec[5] = cos(invec[0]);
}
break;
// TODO: Add the phi into the reparametrization.
//
case STARK_3D:
{
const double Hx = invec[0];
......@@ -311,7 +309,7 @@ void params::to_cartesian(const double* invec, params::input intype,
outvec[0] = invec[0] + cosPhi*invec[1];
outvec[1] = sinPhi*invec[1];
outvec[2] = sqrt(1.0 - outvec[0]*outvec[0] - outvec[1]*outvec[1]);
outvec[3] = invec[0] + cosPhi*invec[1];
outvec[3] = invec[0] - cosPhi*invec[1];
outvec[4] = - sinPhi*invec[1];
outvec[5] = sqrt(1.0 - outvec[3]*outvec[3] - outvec[4]*outvec[4]);
}
......
......@@ -366,9 +366,9 @@ class params
const double cost = cos(theta);
const double sint = sin(theta);
const double temp = cost * vec[0] + sint * vec[1];
const double temp = cost * vec[0] - sint * vec[1];
vec[1] = cost * vec[1] - sint * vec[0];
vec[1] = cost * vec[1] + sint * vec[0];
vec[0] = temp;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment