test_gfield_dipole.cc

Go to the documentation of this file.

 
// gfields
#include "gfield_options.h"
#include "gmagneto.h"
 
// c++
#include <iostream>
 
using namespace std;
 
int main(int argc, char* argv[]) {
 
    // Initialize GOptions (parsed from YAML or another source)
    auto gopts = std::make_shared<GOptions>(argc, argv, gfields::defineOptions());
 
    if (gfields::runFieldQueries(gopts)) { return EXIT_SUCCESS; }
 
    // Create a GMagneto instance to manage fields
    auto magneto = std::make_shared<GMagneto>(gopts);
 
    // Probe whichever fields the YAML configured (e.g. "dipole", "solenoid", "torus"), instead of a
    // hardcoded name, so the example works for every field definition.
    const auto field_names = magneto->getFieldNames();
    if (field_names.empty()) {
        cout << "No fields were configured." << endl;
        return EXIT_SUCCESS;
    }
 
    for (const auto& field_name : field_names) {
        auto field         = magneto->getField(field_name);
        auto field_manager = magneto->getFieldMgr(field_name);
 
        // Non-Doxygen summary:
        // Demonstrate repeated evaluation of the field; this intentionally does not attach the field manager
        // to any volume, as the goal is only to show the API surface.
 
        // pos is a double[3] representing x,y,z
        // iterate through 100 pos that span random positions between 0 and 100 in each coordinate
        for (int i = 0; i < 100; i++) {
            double x = rand() % 100;
            double y = rand() % 100;
            double z = rand() % 100;
 
            double pos[3] = {x, y, z};
            double bfield[3];
            field->GetFieldValue(pos, bfield);
        }
    }
 
    return EXIT_SUCCESS;
}