PB3D
[2.45]
Ideal linear high-n MHD stability in 3-D
|
Calculate inverse of square matrix A
.
More...
Public Member Functions | |
integer function | calc_inv_0d (inv_0D, A) |
private constant version More... | |
integer function | calc_inv_3d (inv_3D, A, n) |
array version More... | |
Calculate inverse of square matrix A
.
This matrix can be defined on a 3-D grid or constant. The storage convention described in eq_vars.eq_2_type is used.
In the former case the size of the matrix (last two indices) should be small, as direct inversion is performed using Cramer's rule.
In the latter case, lapack routines are used.
Definition at line 81 of file num_utilities.f90.
integer function num_utilities::calc_inv::calc_inv_0d | ( | real(dp), dimension(:,:), intent(inout) | inv_0D, |
real(dp), dimension(:,:), intent(in) | A | ||
) |
private constant version
[in,out] | inv_0d | output |
[in] | a | input |
Definition at line 772 of file num_utilities.f90.
integer function num_utilities::calc_inv::calc_inv_3d | ( | real(dp), dimension(:,:,:,:), intent(inout) | inv_3D, |
real(dp), dimension(:,:,:,:), intent(in) | A, | ||
integer, intent(in) | n | ||
) |
array version
[in,out] | inv_3d | output |
[in] | a | input |
[in] | n | size of matrix |
Definition at line 676 of file num_utilities.f90.