Skip to content
Snippets Groups Projects
Commit 2a87c980 authored by Bas van der Tol's avatar Bas van der Tol
Browse files

Use option instead of commenting out code

parent febb86d1
No related branches found
No related tags found
1 merge request!59Resolve "Add stationresponse comparison for oskar"
......@@ -13,15 +13,15 @@ matrix22c_t Element::LocalResponse(real_t time, real_t freq,
id, freq, thetaphi[0], thetaphi[1],
reinterpret_cast<std::complex<double>(&)[2][2]>(result));
// if (options.rotate) {
// vector3r_t up = {0.0, 0.0, 1.0};
// vector3r_t e_phi = normalize(cross(up, direction));
// vector3r_t e_theta = cross(e_phi, direction);
// matrix22r_t rotation;
// rotation[0] = {dot(e_theta, options.north), dot(e_theta, options.east)};
// rotation[1] = {dot(e_phi, options.north), dot(e_phi, options.east)};
// result = result * rotation;
// }
if (options.rotate) {
vector3r_t up = {0.0, 0.0, 1.0};
vector3r_t e_phi = normalize(cross(up, direction));
vector3r_t e_theta = cross(e_phi, direction);
matrix22r_t rotation;
rotation[0] = {dot(e_theta, options.north), dot(e_theta, options.east)};
rotation[1] = {dot(e_phi, options.north), dot(e_phi, options.east)};
result = result * rotation;
}
return result;
}
......
......@@ -89,7 +89,7 @@ int main(int argc, char** argv){
};
real_t freq0 = 50e6;
auto result = station->Response(time, freq, direction, freq0, station0, tile0);
auto result = station->Response(time, freq, direction, freq0, station0, tile0, false);
result_arr[i][j][0][0] = result[0][0];
result_arr[i][j][0][1] = result[0][1];
result_arr[i][j][1][0] = result[1][0];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment