Commit 97119e4d authored by Håkan Wennlöf's avatar Håkan Wennlöf
Browse files

Fixed the issue of scalar field rotation in a more elegant way. 

The issue appeared when the x-coordinate was no longer in the first position in the Point.

Also fixed the unit of the plots in mesh_plotter, so that the user can define the desired one.
parent 3da00663

tools/mesh_converter/MeshConverter.cpp

+15 −10
Original line number Diff line number Diff line
@@ -240,25 +240,20 @@ int main(int argc, char** argv) { 
            throw std::runtime_error("Field and grid file do not match, found " + std::to_string(points.size()) + " and " +

                                     std::to_string(field.size()) + " data points, respectively.");

        }



        auto points_temp = points;

        auto field_temp = field;

        if(rot.at(0) == "-y" || rot.at(0) == "y") {

            for(size_t i = 0; i < points.size(); ++i) {

                points_temp[i].x = points[i].y;

                if(vector_field) {

                field_temp[i].x = field[i].y;

            }

        }

        }

        if(rot.at(0) == "-z" || rot.at(0) == "z") {

            for(size_t i = 0; i < points.size(); ++i) {

                points_temp[i].x = points[i].z;

                if(vector_field) {

                field_temp[i].x = field[i].z;

            }

        }

        }

        if(rot.at(1) == "-x" || rot.at(1) == "x") {

            for(size_t i = 0; i < points.size(); ++i) {

                points_temp[i].y = points[i].x;
@@ -500,13 +495,23 @@ int main(int argc, char** argv) { 
            for(unsigned int j = 0; j < divisions.y(); ++j) {

                for(unsigned int k = 0; k < divisions.z(); ++k) {

                    auto& point = e_field_new_mesh[i * divisions.y() * divisions.z() + j * divisions.z() + k];

                    // We need to convert to framework-internal units:

                    data->push_back(Units::get(point.x, units));

                    // For a vector field, we push three values:

                    LOG(DEBUG) << "Values of data point (" << i << ", " << j << ", " << k << "): " << point;

                    // For a vector field, we push three values:

                    if(quantity == FieldQuantity::VECTOR) {

                        // We need to convert to framework-internal units:

                        data->push_back(Units::get(point.x, units));

                        data->push_back(Units::get(point.y, units));

                        data->push_back(Units::get(point.z, units));

                    } else {

                        // For a scalar field, we need to push only one value, but which one depends on the field rotation.

                        // We need the original x-position, as that is the only filled one in the parsed field

                        if(rot.at(1) == "-x" || rot.at(1) == "x") {

                            data->push_back(Units::get(point.y, units));

                        } else if(rot.at(2) == "-x" || rot.at(2) == "x") {

                            data->push_back(Units::get(point.z, units));

                        } else {

                            data->push_back(Units::get(point.x, units));

                        }

                    }

                }

            }

tools/mesh_converter/MeshPlotter.cpp

+7 −4
Original line number Diff line number Diff line
@@ -261,13 +261,16 @@ int main(int argc, char** argv) { 
                            plot_x,

                            plot_y,

                            sqrt(pow(data->at(base + 0), 2) + pow(data->at(base + 1), 2) + pow(data->at(base + 2), 2)));

                        exfield_map->Fill(plot_x, plot_y, data->at(base + 0));

                        eyfield_map->Fill(plot_x, plot_y, data->at(base + 1));

                        ezfield_map->Fill(plot_x, plot_y, data->at(base + 2));

                        exfield_map->Fill(plot_x, plot_y, static_cast<double>(Units::convert(data->at(base + 0), units)));

                        eyfield_map->Fill(plot_x, plot_y, static_cast<double>(Units::convert(data->at(base + 1), units)));

                        ezfield_map->Fill(plot_x, plot_y, static_cast<double>(Units::convert(data->at(base + 2), units)));



                    } else {

                        // Fill one map with the scalar quantity

                        efield_map->Fill(plot_x, plot_y, data->at(x * ydiv * zdiv + y * zdiv + z));

                        efield_map->Fill(

                            plot_x,

                            plot_y,

                            static_cast<double>(Units::convert(data->at(x * ydiv * zdiv + y * zdiv + z), units)));

                    }

                }

            }