stochastic_heat2d_test
stochastic_heat2d_test,
a Fortran90 code which
calls stochastic_heat2d(), which
solves the steady state heat equation in a 2D
rectangular region with a stochastic heat diffusivity,
using the finite difference method (FDM), and stochastic model
from Babuska, Nobile and Tempone, using GNUPLOT to illustrate
the results.
Licensing:
The information on this web page is distributed under the MIT license.
Related Data and codes:
stochastic_heat2d,
a Fortran90 code which
solves the steady state heat equation in a 2D
rectangular region with a stochastic heat diffusivity,
using the finite difference method (FDM), and stochastic model
from Babuska, Nobile and Tempone, using GNUPLOT to illustrate
the results.
Source Code:

solution_data.txt
graphics data for an image of the solution corresponding to an example selection of OMEGA.

solution_commands.txt
GNUPLOT commands to create an image of the solution corresponding to an example selection of OMEGA.

solution.png
an image of the solution corresponding to an example selection of OMEGA.

umean_data.txt
graphics data to create an image of mean temperature as a function of OMEGA(1) and OMEGA(2).

umean_commands.txt
GNUPLOT commands to create an image of mean temperature as a function of OMEGA(1) and OMEGA(2).

umean.png
an image of mean temperature as a function of OMEGA(1) and OMEGA(2).
List of Routines:

DIFFUSIVITY_2D_BNT evaluates a 2D stochastic diffusivity function.

GET_UNIT returns a free Fortran unit number.

INTERIOR sets up the matrix and right hand side at interior nodes.

R8_UNIFORM_01 returns a unit pseudorandom R8.

R8MAT_FS factors and solves a system with one right hand side.

R8VEC_LINSPACE creates a vector of linearly spaced values.

R8VEC_MESH_2D creates a 2D mesh from X and Y vectors.

R8VEC_NORMAL_01 returns a unit pseudonormal R8VEC.

R8VEC_PRINT prints an R8VEC.

R8VEC_UNIFORM_01 returns a unit pseudorandom R8VEC.

STOCHASTIC_HEAT2D solves the steady 2D heat equation.

TIMESTAMP prints the current YMDHMS date as a time stamp.
You can go up one level to
the Fortran90 source codes.
Last revised on 25 August 2013.