06 October 2025   6:50:34.201 PM

stochastic_heat2d_test():
  Fortran90 version
  test  stochastic_heat2d().

TEST01:
  Consider the steady heat equation in the unit square,
  with 0 Dirichlet boundary conditions, 
  and a heat source term F that is a Gaussian centered at (0.60,0.80).

  Model the diffusivity coefficient as spatially varying,
  with a stochastic dependence on parameters OMEGA(1:4),
  as described in Babuska, Nobile, Tempone (BNT).

  Compute and display the solution U for a given choice
  of the parameters OMEGA.
 
  Sampled OMEGA values:
 
         1:   -2.5768189    
         2:   0.74496194    
         3:   0.89135992    
         4:    1.5486025    
 
  Created graphics data file "solution_data.txt".
  Created graphics command file "solution_commands.txt".
 
  Mean value of U is   0.101042    

TEST02:
  Fix OMEGA(3) = 4, OMEGA(4) = 0, and
  examine dependence of average temperature on OMEGA(1) and OMEGA(2)
  over the range [-10,+10].

  Omega(3) fixed at    4.00000    
  Omega(4) fixed at    0.00000    
 
  Created graphics data file "umean_data.txt".
  Created graphics command file "umean_commands.txt".

  U_Mean_Max =   0.641953    

stochastic_heat2d_test():
  Normal end of execution.

06 October 2025   6:50:36.166 PM