27-Sep-2022 14:22:56

plasma_matrix_test():
  MATLAB/Octave version 4.2.2
  Test plasma_matrix()

plasma_matrix_test01():
  Linear order of spatial grid = 5
  Number of grid points is 25
  Approximate number of nonzero matrix entries is 125.
  Full storage matrix would require 625 entries.

  Norm of right hand side b is 6.1997

  Sparsity image saved as "plasma_5_sparsity.png"

  Norm of solution is 0.315693

plasma_matrix_test01():
  Linear order of spatial grid = 51
  Number of grid points is 2601
  Approximate number of nonzero matrix entries is 13005.
  Full storage matrix would require 6765201 entries.

  Norm of right hand side b is 19.9934

  Sparsity image saved as "plasma_51_sparsity.png"

  Norm of solution is 14.6975

plasma_matrix_test02():
  Linear order of spatial grid = 101
  Number of grid points is 10201
  Approximate number of nonzero matrix entries is 51005.
  Full storage matrix would require 104060401 entries.

  Norm of right hand side b is 27.7259

  Norm of backslash solution x is 29.2756
bicg stopped at iteration 20 without converging to the desired tolerance 1.000000e-06
because the maximum number of iterations was reached. The iterate returned (number 20) has relative residual 2.477131e-02

  Norm of bicg solution x is 24.4076
bicg converged at iteration 39 to a solution with relative residual 9.727749e-07

  Norm of bicg+ilu solution x is 29.2756

plasma_matrix_test():
  Normal end of execution.

27-Sep-2022 14:23:59