distance_norm2_eigen()

template<typename T >

T cntr::distance_norm2_eigen	(	int	tstp,
		herm_matrix_timestep< T > &	g1,
		herm_matrix< T > &	g2
	)

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

Purpose

Evaluate the Euclidean norm between 'herm_matrix_timestep' ( \(g_1 \)) and 'herm_matrix' ( \( g_2\)) at a given time step 'tstp' using the standard 'norm'-routine. To evaluate the norm, the elements of retarded, lesser and left-mixing components at the time step is used. The norm is not normalized per elements, but it is the summention of all the elements.

Parameters

tstp	time step
g1	object of the class 'herm_matrix_timestep'
g2	object of the class 'herm_matrix'

Definition at line 1960 of file cntr_utilities_impl.hpp.

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

                                                                               {
   int size1=g2.size1(),size2=g2.size2(),ntau=g2.ntau();
   T err=0.0;

   assert(g1.tstp_==tstp && g1.tstp_<=g2.nt());
   assert(g1.ntau_==g2.ntau() && g1.size1_==g2.size1() && g1.size2_==g2.size2() );
   cdmatrix matg1(size1,size2);
   cdmatrix matg2(size1,size2);

   if(tstp==-1){
     for(int i=0;i<=ntau;i++){
        g1.get_mat(i,matg1);
        g2.get_mat(i,matg2);
        err+=(matg1-matg2).norm();
     }
   }else{
     // Ret
     for(int i=0;i<=tstp;i++){
        g1.get_ret_tstp_t(i,matg1);
        g2.get_ret(tstp,i,matg2);
        err+=(matg1-matg2).norm();
     }
     // tv
     for(int i=0;i<=ntau;i++){
        g1.get_tv(i,matg1);
        g2.get_tv(tstp,i,matg2);
        err+=(matg1-matg2).norm();
     }
     // Les
     for(int i=0;i<=tstp;i++){
        g1.get_les_t_tstp(i,matg1);
        g2.get_les(i,tstp,matg2);
        err+=(matg1-matg2).norm();
     }
   }
}

Here is the call graph for this function: