17 December  2023   7:56:55.746 AM      

stochastic_heat2d_test():
  Fortran77 version
  Test stochastic_heat2d().

TEST01:
  Consider the steady heat equation in the unit square,
  with 0 Dirichlet boundary conditions, 
  and a Gaussian heat source term F 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:    3.3580805    
         2:  -0.94553760    
         3:   -1.1321196    
         4:  -0.46224826    
 
  Created graphics data file "solution_data.txt".
  Created graphics command file "solution_commands.txt".
 
  Mean value of U is   0.337698    

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.

17 December  2023   7:56:59.066 AM