Implement unit tests for solovev_free_bdy

Author: jonathanschilling

src/vmecpp/cpp/vmecpp/vmec/output_quantities/BUILD.bazel

@@ -41,6 +41,7 @@ cc_test(
     data = [
         "//vmecpp/test_data:solovev",
         "//vmecpp/test_data:solovev_no_axis",
+        "//vmecpp/test_data:solovev_free_bdy",
         "//vmecpp/test_data:cth_like_fixed_bdy",
         "//vmecpp/test_data:cth_like_fixed_bdy_nzeta_37",
         "//vmecpp/test_data:cma",

src/vmecpp/cpp/vmecpp/vmec/output_quantities/output_quantities_test.cc

@@ -415,6 +415,7 @@ INSTANTIATE_TEST_SUITE_P(
     TestOutputQuantities, WOutFileContentsTest,
     Values(DataSource{.identifier = "solovev", .tolerance = 5.0e-7},
            DataSource{.identifier = "solovev_no_axis", .tolerance = 5.0e-7},
+           DataSource{.identifier = "solovev_free_bdy", .tolerance = 5.0e-7},
            DataSource{.identifier = "cth_like_fixed_bdy", .tolerance = 1.0e-6},
            DataSource{.identifier = "cth_like_fixed_bdy_nzeta_37",
                       .tolerance = 1.0e-6},

tests/test_free_boundary.py

@@ -68,6 +68,21 @@ def test_run_free_boundary_from_response_table():
     )
 
 
+@pytest.mark.skip(
+    reason="Axisymmetric free-boundary code path not yet fully implemented (segfaults)"
+)
+def test_run_free_boundary_solovev():
+    makegrid_params = vmecpp.MakegridParameters.from_file(
+        TEST_DATA_DIR / "makegrid_parameters_solovev.json"
+    )
+    response = vmecpp.MagneticFieldResponseTable.from_coils_file(
+        TEST_DATA_DIR / "coils.solovev", makegrid_params
+    )
+    vmec_input = vmecpp.VmecInput.from_file(TEST_DATA_DIR / "solovev_free_bdy.json")
+    vmec_output = vmecpp.run(vmec_input, response, verbose=False)
+    assert vmec_output.wout.volume == pytest.approx(128.8377, rel=1e-3, abs=1e-3)
+
+
 def test_raise_invalid_nzeta():
     makegrid_params = vmecpp.MakegridParameters.from_file(
         TEST_DATA_DIR / "makegrid_parameters_cth_like.json"