old_contents: from pygbe.util import an_solution
  from convergence_lspr import (run_convergence, picklesave, pickleload,
                                report_results, mesh)
  
  def main():
      print('{:-^60}'.format('Running sphere_lspr test'))
      try:
          test_outputs = pickleload()
      except FileNotFoundError:
          test_outputs = {}
  
      problem_folder = 'input_files'
  
      # dirichlet_surface
      param = 'sphere_complex.param'
      test_name = 'sphere_complex'
      if test_name not in test_outputs.keys():
         N, iterations, expected_rate, Cext, Time = run_convergence(
              mesh, test_name, problem_folder, param)
          test_outputs[test_name] = [N, iterations, expected_rate, Cext[0], Time]
  
      picklesave(test_outputs)
  
      # load data for analysis
      N, iterations, expected_rate = test_outputs['sphere_complex'][:3]
      Cext = test_outputs['sphere_complex'][3]
      Time = test_outputs['sphere_complex'][-1]
  
      total_time = Time
      
      #This test is for 10 nm radius silver sphere in water, at wavelength 380 nm 
      radius = 10.
      wavelength = 380.
      diel_out = 1.7972083599999999 + 1j * 8.504766399999999e-09 #water value extrapolated
      diel_in = -3.3876520488233184 + 1j * 0.19220746083441781 #silver value extrapolated
  
      analytical = an_solution.Cext_analytical(radius, wavelength, diel_out, diel_in)
  
      error = abs(Cext - analytical) / abs(analytical)
  
      report_results(error,
                     N,
                     expected_rate,
                     iterations,
                     Cext,
                     analytical,
                     total_time,
                     test_name='sphere_complex')
  
  
  if __name__ == "__main__":
      from check_for_meshes import check_mesh
      check_mesh()
      main()
  
new_contents: from pygbe.util import an_solution
  from convergence_lspr import (run_convergence, picklesave, pickleload,
                                report_results, mesh)
  
  def main():
      print('{:-^60}'.format('Running sphere_lspr test'))
      try:
          test_outputs = pickleload()
      except FileNotFoundError:
          test_outputs = {}
  
      problem_folder = 'input_files'
  
      # dirichlet_surface
      param = 'sphere_complex.param'
      test_name = 'sphere_complex'
      if test_name not in test_outputs.keys():
         N, iterations, expected_rate, Cext, Time = run_convergence(
              mesh, test_name, problem_folder, param)
         test_outputs[test_name] = [N, iterations, expected_rate, Cext[0], Time]
  
      picklesave(test_outputs)
  
      # load data for analysis
      N, iterations, expected_rate = test_outputs['sphere_complex'][:3]
      Cext = test_outputs['sphere_complex'][3]
      Time = test_outputs['sphere_complex'][-1]
  
      total_time = Time
      
      #This test is for 10 nm radius silver sphere in water, at wavelength 380 nm 
      radius = 10.
      wavelength = 380.
      diel_out = 1.7972083599999999 + 1j * 8.504766399999999e-09 #water value extrapolated
      diel_in = -3.3876520488233184 + 1j * 0.19220746083441781 #silver value extrapolated
  
      analytical = an_solution.Cext_analytical(radius, wavelength, diel_out, diel_in)
  
      error = abs(Cext - analytical) / abs(analytical)
  
      report_results(error,
                     N,
                     expected_rate,
                     iterations,
                     Cext,
                     analytical,
                     total_time,
                     test_name='sphere_complex')
  
  
  if __name__ == "__main__":
      main()
  
subject: Fix indentation error