Adding interface functions to LAW=2

Author: whaeck

ReleaseNotes.md

@@ -4,6 +4,9 @@ Given here are some release notes for ENDFtk.
 ## [ENDFtk v1.1.1](https://github.com/njoy/ENDFtk/pull/xxx)
 This update removes the regions() and pairs() interface functions on the TAB1 record interface functions that are unused. The removal of these interface functions has no impact on the Python interface as these interface functions were not included on the Python side. Miscellaneous documentation updates were made as well.
 
+This update also adds NBT(), INT(), boundaries(), interpolants(), NR(), numberInterpolationRegions() on a number of ENDFtk objects that mimic the behaviour of a TAB1 record but that are not actual TAB1 records:
+- TwoBodyScattering::TabulatedDistribution in MF6 and MF26
+
 In addition, the following issues were corrected:
   - A minor bug in the rectangular matrix covariance block was corrected. The values for the row and column energies are lifted out of a larger array using the std::ranges::take and std::ranges::drop function. For the column energies, we forgot to properly end the sequence. As a result, the end() iterator of the range did not point to the end of the column energies but to the end of the covariance values, which is now corrected.
   - In MF8 MT457 DiscreteSpectrum, NT=8 (which can occur for electrons, i.e. STYP=8), was explicitly disallowed by ENDFtk. This was corrected and a new constructor reflecting this usage was added as well.

python/src/section/6/LAW2/TabulatedDistribution.python.cpp

@@ -111,6 +111,7 @@ void wrapTabulatedDistribution( python::module& module, python::module& ) {
   );
 
   // add standard component definitions
+  addStandardInterpolationTableDefinitions< Component >( component );
   addStandardComponentDefinitions< Component >( component );
 }
 

python/test/MF6/LAW2/Test_ENDFtk_MF6_LAW2_TabulatedDistribution.py

@@ -50,6 +50,17 @@ def verify_chunk_lang12( self, chunk ) :
             self.assertAlmostEqual( 4., chunk.probabilities[1] )
             self.assertAlmostEqual( 6., chunk.probabilities[2] )
 
+            self.assertEqual( 1, chunk.NR )
+            self.assertEqual( 1, chunk.number_interpolation_regions )
+            self.assertEqual( 1, len( chunk.INT ) )
+            self.assertEqual( 1, len( chunk.NBT ) )
+            self.assertEqual( 2, chunk.INT[0] )
+            self.assertEqual( 3, chunk.NBT[0] )
+            self.assertEqual( 1, len( chunk.interpolants ) )
+            self.assertEqual( 1, len( chunk.boundaries ) )
+            self.assertEqual( 2, chunk.interpolants[0] )
+            self.assertEqual( 3, chunk.boundaries[0] )
+
             self.assertEqual( 2, chunk.NC )
 
             # verify string
@@ -82,6 +93,17 @@ def verify_chunk_lang14( self, chunk ) :
             self.assertAlmostEqual( 4., chunk.probabilities[1] )
             self.assertAlmostEqual( 6., chunk.probabilities[2] )
 
+            self.assertEqual( 1, chunk.NR )
+            self.assertEqual( 1, chunk.number_interpolation_regions )
+            self.assertEqual( 1, len( chunk.INT ) )
+            self.assertEqual( 1, len( chunk.NBT ) )
+            self.assertEqual( 4, chunk.INT[0] )
+            self.assertEqual( 3, chunk.NBT[0] )
+            self.assertEqual( 1, len( chunk.interpolants ) )
+            self.assertEqual( 1, len( chunk.boundaries ) )
+            self.assertEqual( 4, chunk.interpolants[0] )
+            self.assertEqual( 3, chunk.boundaries[0] )
+
             self.assertEqual( 2, chunk.NC )
 
             # verify string

src/ENDFtk/section/6/DiscreteTwoBodyScattering/TabulatedDistribution.hpp

@@ -90,6 +90,42 @@ class ENDFTK_PYTHON_EXPORT TabulatedDistribution : protected ListRecord {
     return this->F();
   }
 
+  /**
+   *  @brief Return the number of interpolation ranges
+   */
+  static constexpr long NR() { return 1; }
+
+  /**
+   *  @brief Return the number of interpolation ranges
+   */
+  static constexpr long numberInterpolationRegions() { return NR(); }
+
+  /**
+   *  @brief Return the interpolants
+   */
+  auto INT() const {
+
+    return ranges::cpp20::views::single( this->LANG() - 10 );
+  }
+
+  /**
+   *  @brief Return the interpolants
+   */
+  auto interpolants() const { return this->INT(); }
+
+  /**
+   *  @brief Return the boundaries
+   */
+  auto NBT() const {
+
+    return ranges::cpp20::views::single( this->numberCosineValues() );
+  }
+
+  /**
+   *  @brief Return the boundaries
+   */
+  auto boundaries() const { return this->NBT(); }
+
   using ListRecord::NC;
   using ListRecord::print;
 };

src/ENDFtk/section/6/DiscreteTwoBodyScattering/TabulatedDistribution/test/TabulatedDistribution.test.cpp

@@ -223,6 +223,17 @@ void verifyChunkLANG12( const TabulatedDistribution& chunk ) {
   CHECK_THAT( 4., WithinRel( chunk.probabilities()[1] ) );
   CHECK_THAT( 6., WithinRel( chunk.probabilities()[2] ) );
 
+  CHECK( 1 == chunk.NR() );
+  CHECK( 1 == chunk.numberInterpolationRegions() );
+  CHECK( 1 == chunk.INT().size() );
+  CHECK( 1 == chunk.interpolants().size() );
+  CHECK( 1 == chunk.NBT().size() );
+  CHECK( 1 == chunk.boundaries().size() );
+  CHECK( 2 == chunk.INT()[0] );
+  CHECK( 2 == chunk.interpolants()[0] );
+  CHECK( 3 == chunk.NBT()[0] );
+  CHECK( 3 == chunk.boundaries()[0] );
+
   CHECK( 2 == chunk.NC() );
 }
 
@@ -258,6 +269,17 @@ void verifyChunkLANG14( const TabulatedDistribution& chunk ) {
   CHECK_THAT( 4., WithinRel( chunk.probabilities()[1] ) );
   CHECK_THAT( 6., WithinRel( chunk.probabilities()[2] ) );
 
+  CHECK( 1 == chunk.NR() );
+  CHECK( 1 == chunk.numberInterpolationRegions() );
+  CHECK( 1 == chunk.INT().size() );
+  CHECK( 1 == chunk.interpolants().size() );
+  CHECK( 1 == chunk.NBT().size() );
+  CHECK( 1 == chunk.boundaries().size() );
+  CHECK( 4 == chunk.INT()[0] );
+  CHECK( 4 == chunk.interpolants()[0] );
+  CHECK( 3 == chunk.NBT()[0] );
+  CHECK( 3 == chunk.boundaries()[0] );
+
   CHECK( 2 == chunk.NC() );
 }