solve Subroutine

public subroutine solve(matrix_A, matrix_B, x, settings, snapshots, snap_times)

Solve the initial value problem

Arguments

Type IntentOptional Attributes Name
type(matrix_t) :: matrix_A

FEM matrix A

type(matrix_t) :: matrix_B

FEM matrix B

complex(kind=dp), intent(inout), dimension(:) :: x

initial condition, gets updated with final result

type(settings_t), intent(in) :: settings

settings

complex(kind=dp), intent(out), optional, dimension(:,:) :: snapshots

optionally save every n-th step in this 2D array

real(kind=dp), intent(out), optional, dimension(:) :: snap_times

snapshot times