21 June 2024 01:26:01 PM

ppc_porous_medium_test ( ):
  C version.
  Test ppc_mesh().

pme1_solve():
  Solve porous medium problem #1:

  Integrate over time t for 0 <= t <= T = 1
  Take m = 40 equal steps in time
  dt = 0.025
  a = -1 <= x <= b = 1
  Use n = 20 subintervals in space 
  dx = 0.1
  r = dt / 2 / dx^2 = 1.25
  s = sqrt(r) = 1.11803
  Geomview script written to pme1.gv
  Maple script written to pme1.mpl
  Matlab script written to pme1.m

  Maximum difference between computed and exact: 0.197579

pme2_solve():
  Solve porous medium problem #2:

  Integrate over time t for 0 <= t <= T = 1
  Take m = 40 equal steps in time
  dt = 0.025
  a = -1 <= x <= b = 1
  Use n = 20 subintervals in space 
  dx = 0.1
  r = dt / 2 / dx^2 = 1.25
  s = sqrt(r) = 1.11803
  Geomview script written to pme2.gv
  Maple script written to pme2.mpl
  Matlab script written to pme2.m

ppc_porous_medium_test ( ):
  Normal end of execution.
21 June 2024 01:26:01 PM