NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ distance_norm2_ret() [2/9]

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

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

Purpose

Evaluate the Euclidean norm between the retarded components of a herm_matrix_timestep \(g_1\) and a herm_matrix ( \(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_timestep

g2

herm_matrix

Definition at line 792 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().

792  {
793  assert(g1.size1() == g2.size1());
794  assert(g1.ntau() == g2.ntau());
795  assert(g1.tstp() == tstp);
796  assert(g2.nt() >= tstp);
797  herm_matrix_timestep_view<T> g1_view(tstp, g1);
798  herm_matrix_timestep_view<T> g2_view(tstp, g2);
799  if (g1.size1() == 1)
800  return distance_norm2_ret_dispatch<T, 1>(tstp, g1_view, g2_view);
801  else
802  return distance_norm2_ret_dispatch<T, LARGESIZE>(tstp, g1_view, g2_view);
803 }
+ Here is the call graph for this function: