MCTDH-X
v2.22
|
This module collects the analysis routines dedicated to two-dimensional R-MCTDHB simulations. More...
Public Member Functions | |
subroutine | get_correlationsskewdiagonal2d (time, PSI, mode, rho_jk, rho_ijkl, dilation) |
This routine computes the skew diagonal of 2D correlation functions G1(r,r'=-r), and G2(r,r'=-r), where r=(x,y), and writes them to a file. The first four columns contain x,y,x',y'. More... | |
subroutine | get_corr_slice (time, PSI, mode, rho_jk, rho_ijkl, x1const, x1slice, y1const, y1slice, x2const, x2slice, y2const, y2slice, dilation) |
This routine computes cuts of the correlation functions in two spatial dimensions G1(r,r), G1(r,r'), G1(r',r') and G2(r,r',r',r), where r=(x,y), and writes them to a file. The user has to choose two of the coordinates to stay constant, to obtain a 3D plot. More... | |
subroutine | get_vtrapprojection (PSI, Rho_IJ, time, dir) |
This routine calculates a projected potential where the second spatial dimension is integrated weighted with the 1-body-density. More... | |
subroutine | get_lz (PSI, NOcc, NatVec, time) |
This routine Calculates the angular momentum expectation and matrix elements of a given 2-dimensional wavefunction. More... | |
This module collects the analysis routines dedicated to two-dimensional R-MCTDHB simulations.
subroutine analysis_routines_2d::get_corr_slice | ( | real*8, intent(in) | time, |
complex*16, dimension(ndvr_x*ndvr_y*ndvr_z,morb) | PSI, | ||
integer, intent(in) | mode, | ||
complex*16, dimension(morb,morb), intent(in) | rho_jk, | ||
complex*16, dimension(morb,morb,morb,morb), intent(in) | rho_ijkl, | ||
logical, intent(in) | x1const, | ||
real*8, intent(in) | x1slice, | ||
logical, intent(in) | y1const, | ||
real*8, intent(in) | y1slice, | ||
logical, intent(in) | x2const, | ||
real*8, intent(in) | x2slice, | ||
logical, intent(in) | y2const, | ||
real*8, intent(in) | y2slice, | ||
integer, intent(in) | dilation | ||
) |
This routine computes cuts of the correlation functions in two spatial dimensions G1(r,r), G1(r,r'), G1(r',r') and G2(r,r',r',r), where r=(x,y), and writes them to a file. The user has to choose two of the coordinates to stay constant, to obtain a 3D plot.
References auxiliary_routines::get_ijk_from_m(), auxiliary_analysis_routines::get_maximum(), auxiliary_routines::get_mom(), auxiliary_routines::getdoubleasstring(), and auxiliary_analysis_routines::getsliceasstring().
Referenced by analysis_core::apply_analysis().
subroutine analysis_routines_2d::get_correlationsskewdiagonal2d | ( | real*8, intent(in) | time, |
complex*16, dimension(ndvr_x*ndvr_y*ndvr_z,morb) | PSI, | ||
integer, intent(in) | mode, | ||
complex*16, dimension(morb,morb), intent(in) | rho_jk, | ||
complex*16, dimension(morb,morb,morb,morb), intent(in) | rho_ijkl, | ||
integer, optional | dilation | ||
) |
This routine computes the skew diagonal of 2D correlation functions G1(r,r'=-r), and G2(r,r'=-r), where r=(x,y), and writes them to a file. The first four columns contain x,y,x',y'.
References auxiliary_routines::get_ijk_from_m(), auxiliary_analysis_routines::get_maximum(), auxiliary_routines::get_mom(), and auxiliary_routines::getdoubleasstring().
Referenced by analysis_core::apply_analysis().
subroutine analysis_routines_2d::get_lz | ( | complex*16, dimension(ndvr_x*ndvr_y*ndvr_z,morb), intent(in) | PSI, |
real*8, dimension(morb), intent(in) | NOcc, | ||
complex*16, dimension(morb,morb), intent(in) | NatVec, | ||
real*8, intent(in) | time | ||
) |
This routine Calculates the angular momentum expectation and matrix elements of a given 2-dimensional wavefunction.
References kineticenergyaction::get_kineticenergyaction_l_fft(), auxiliary_routines::get_nos(), and function_library::normvxz().
Referenced by analysis_core::apply_analysis().
subroutine analysis_routines_2d::get_vtrapprojection | ( | complex*16, dimension(ndvr_x*ndvr_y*ndvr_z,morb) | PSI, |
complex*16, dimension(morb,morb) | Rho_IJ, | ||
real*8 | time, | ||
character(len=1) | dir | ||
) |
This routine calculates a projected potential where the second spatial dimension is integrated weighted with the 1-body-density.
References auxiliary_routines::getdoubleasstring().
Referenced by analysis_core::apply_analysis().