diff --git a/tests/regression_tests/torus/large_major/__init__.py b/tests/regression_tests/torus/large_major/__init__.py
new file mode 100644
index 0000000000..e69de29bb2
diff --git a/tests/regression_tests/torus/large_major/inputs_true.dat b/tests/regression_tests/torus/large_major/inputs_true.dat
new file mode 100644
index 0000000000..ae3f72d0fc
--- /dev/null
+++ b/tests/regression_tests/torus/large_major/inputs_true.dat
@@ -0,0 +1,37 @@
+<?xml version='1.0' encoding='utf-8'?>
+<geometry>
+  <cell id="1" material="void" region="-1" universe="1" />
+  <cell id="2" material="1" region="1 -2" universe="1" />
+  <cell id="3" material="2" region="2 -3" universe="1" />
+  <surface coeffs="0.0 0.0 0.0 1000.0 30.0 30.0" id="1" type="z-torus" />
+  <surface coeffs="0.0 0.0 0.0 1000.0 35.0 35.0" id="2" type="z-torus" />
+  <surface boundary="vacuum" coeffs="0.0 0.0 0.0 1000.0 40.0 40.0" id="3" type="z-torus" />
+</geometry>
+<?xml version='1.0' encoding='utf-8'?>
+<materials>
+  <material id="1">
+    <density units="g/cm3" value="1.0" />
+    <nuclide ao="1.0" name="W184" />
+  </material>
+  <material id="2">
+    <density units="g/cm3" value="5.0" />
+    <nuclide ao="1.0" name="Fe56" />
+  </material>
+</materials>
+<?xml version='1.0' encoding='utf-8'?>
+<settings>
+  <run_mode>fixed source</run_mode>
+  <particles>1000</particles>
+  <batches>10</batches>
+  <source strength="1.0">
+    <space type="point">
+      <parameters>-1000.0 0 0</parameters>
+    </space>
+  </source>
+</settings>
+<?xml version='1.0' encoding='utf-8'?>
+<tallies>
+  <tally id="1">
+    <scores>flux</scores>
+  </tally>
+</tallies>
diff --git a/tests/regression_tests/torus/large_major/results_true.dat b/tests/regression_tests/torus/large_major/results_true.dat
new file mode 100644
index 0000000000..cb9f54970c
--- /dev/null
+++ b/tests/regression_tests/torus/large_major/results_true.dat
@@ -0,0 +1,3 @@
+tally 1:
+9.675396E+02
+9.363406E+04
diff --git a/tests/regression_tests/torus/large_major/test.py b/tests/regression_tests/torus/large_major/test.py
new file mode 100644
index 0000000000..256716f4e1
--- /dev/null
+++ b/tests/regression_tests/torus/large_major/test.py
@@ -0,0 +1,41 @@
+import openmc
+import pytest
+
+from tests.testing_harness import PyAPITestHarness
+
+
+@pytest.fixture
+def model():
+    model = openmc.Model()
+    tungsten = openmc.Material()
+    tungsten.set_density('g/cm3', 1.0)
+    tungsten.add_nuclide('W184', 1.0)
+    ss = openmc.Material()
+    ss.set_density('g/cm3', 5.0)
+    ss.add_nuclide('Fe56', 1.0)
+    model.materials.extend([tungsten, ss])
+
+    # Create nested torii with very large major radii
+    R = 1000.0
+    vacuum = openmc.ZTorus(a=R, b=30.0, c=30.0)
+    first_wall = openmc.ZTorus(a=R, b=35.0, c=35.0)
+    vessel = openmc.ZTorus(a=R, b=40.0, c=40.0, boundary_type='vacuum')
+    cell1 = openmc.Cell(region=-vacuum)
+    cell2 = openmc.Cell(fill=tungsten, region=+vacuum & -first_wall)
+    cell3 = openmc.Cell(fill=ss, region=+first_wall & -vessel)
+    model.geometry = openmc.Geometry([cell1, cell2, cell3])
+
+    model.settings.run_mode ='fixed source'
+    model.settings.particles = 1000
+    model.settings.batches = 10
+    model.settings.source = openmc.Source(space=openmc.stats.Point((-R, 0, 0,)))
+
+    tally = openmc.Tally()
+    tally.scores = ['flux']
+    model.tallies.append(tally)
+    return model
+
+
+def test_torus_large_major(model):
+    harness = PyAPITestHarness('statepoint.10.h5', model)
+    harness.main()