new_contents: import numpy as np
  import unittest
  from desc.io import InputReader
  from desc.equilibrium import Equilibrium
  from desc.profiles import PowerSeriesProfile
  
  
  class TestProfiles(unittest.TestCase):
      def test_same_result(self):
          input_path = "examples/DESC/SOLOVEV"
          ir = InputReader(input_path)
  
          eq1 = Equilibrium(ir.inputs[-1])
          eq2 = eq1.copy()
          eq2.pressure = eq1.pressure.to_spline()
          eq2.iota = eq1.iota.to_spline()
  
          eq1.solve()
          eq2.solve()
  
          np.testing.assert_allclose(
              eq1.x,
              eq2.x,
              rtol=1e-05,
              atol=1e-08,
          )
  
      def test_close_values(self):
  
          pp = PowerSeriesProfile(np.array([0, 2, 4]), np.array([1, -2, 1]))
          sp = pp.to_spline()
          x = np.linspace(0, 1, 100)
  
          np.testing.assert_allclose(pp.compute(x), sp.compute(x), rtol=1e-5, atol=1e-3)
  
          pp1 = sp.to_powerseries(order=4)
          np.testing.assert_allclose(pp.coeffs, pp1.coeffs, rtol=1e-5, atol=1e-2)
  
subject: Add tests for profile classes