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:

• stochastic_heat2d_test.f90, calls all the tests.
• stochastic_heat2d_test.sh, runs all the tests.
• stochastic_heat2d_test.txt, the output file.

• 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.