Field Examples

Field types provide grid storage for PDE and vector-calculus examples.

Scalar Field

num::Grid3D grid{64, 64, 64, 1.0 / 63.0};
 
num::ScalarField3D rho(grid);
rho.fill([&](double x, double y, double z) {
    return std::exp(-40.0 * ((x - 0.5) * (x - 0.5)
                           + (y - 0.5) * (y - 0.5)
                           + (z - 0.5) * (z - 0.5)));
});

Vector Field

num::VectorField3D v(grid);
v.x.fill(0.0);
v.y.fill(0.0);
v.z.fill(1.0);

Poisson Utility

num::ScalarField3D phi(grid);
 
num::SolverResult info =
    num::FieldSolver::solve_poisson(phi, rho, 1e-8, 1000);

Internally this wraps the finite-difference operator with num::operators::make_op and solves with conjugate gradients.

Magnetic Field Utility

num::VectorField3D J(grid);
load_current_density(J);
 
num::VectorField3D B = num::MagneticSolver::solve_magnetic_field(J);

This computes vector potential components from Poisson solves and returns \(\mathbf{B}=\nabla\times\mathbf{A}\).