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.