NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ distance_norm2_ret() [3/9]

template<typename T >
T cntr::distance_norm2_ret ( int  tstp,
herm_matrix< T > &  g1,
herm_matrix_timestep< T > &  g2 
)

Evaluate the Euclidean norm between the retarded components a herm_matrix and herm_matrix_timestep at a given time step.

Purpose

Evaluate the Euclidean norm between the retarded components of a herm_matrix \(g_1\) and a herm_matrix_timestep ( \(g_2\)) at a given time step (tstp). The norm is not normalized per elements, but it is the summention of all the elements.

Parameters
tstp

time step

g1

herm_matrix

g2

herm_matrix_timestep

Definition at line 824 of file cntr_utilities_impl.hpp.

References cntr::herm_matrix< T >::nt(), cntr::herm_matrix_timestep< T >::ntau(), cntr::herm_matrix< T >::ntau(), cntr::herm_matrix_timestep< T >::size1(), cntr::herm_matrix< T >::size1(), and cntr::herm_matrix_timestep< T >::tstp().

824  {
825  assert(g1.size1() == g2.size1());
826  assert(g1.ntau() == g2.ntau());
827  assert(g1.nt() >= tstp);
828  assert(g2.tstp() == tstp);
829  herm_matrix_timestep_view<T> g1_view(tstp, g1);
830  herm_matrix_timestep_view<T> g2_view(tstp, g2);
831  if (g1.size1() == 1)
832  return distance_norm2_ret_dispatch<T, 1>(tstp, g1_view, g2_view);
833  else
834  return distance_norm2_ret_dispatch<T, LARGESIZE>(tstp, g1_view, g2_view);
835 }
+ Here is the call graph for this function: