left_multiply_hermconj()

template<typename T >

void cntr::herm_matrix_timestep< T >::left_multiply_hermconj	(	int	tstp,
		function< T > &	ft,
		T	weight = 1.0
	)

Left-multiplies the herm_matrix_timestep with the hermitian conjugate of a contour function.

Purpose

Performs the operation \(C(t,t^\prime) \rightarrow w F^\ddagger(t)C(t,t^\prime)\), where \(C(t,t^\prime)\) is the herm_matrix_timestep, \(F(t)\) is a contour function and \(w\) is a real weight. The operation is performed at given time step tstp.

Parameters

tstp	[int] The time step at which \(F^\ddagger(t)\) and \(C(t,t^\prime)\) are multiplied. This is a dummy argument in release mode.
ft	[function] The contour function.
weight	[T] The weight as above.

Definition at line 1630 of file cntr_herm_matrix_timestep_impl.hpp.

References cntr::function< T >::nt(), cntr::function< T >::ptr(), cntr::function< T >::size1(), and cntr::function< T >::size2().

                                                      {
    assert(tstp == tstp_);
    int m;
    cplx *xtemp, *ftemp, *x0, *f0, *fcc;
    xtemp = new cplx[element_size_];
    fcc = new cplx[element_size_];
    assert(tstp >= -1 && ft.nt() >= tstp &&
           ft.size1() == size1_ && ft.size2() == size2_ &&
      "tstp >= -1 && ft.nt() >= tstp && ft.size1() == size1_ && ft.size2() == size2_");

    f0 = ft.ptr(-1);
    element_conj<T, LARGESIZE>(size1_, fcc, f0);
    if (tstp == -1) {
        x0 = data_;
        for (m = 0; m <= ntau_; m++) {
            element_mult<T, LARGESIZE>(size1_, xtemp, fcc,
                                       x0 + m * element_size_);
            element_smul<T, LARGESIZE>(size1_, xtemp, weight);
            element_set<T, LARGESIZE>(size1_, x0 + m * element_size_, xtemp);
        }
    } else {
        ftemp = ft.ptr(tstp);
        element_conj<T, LARGESIZE>(size1_, fcc, ftemp);
        x0 = data_;
        for (m = 0; m <= tstp; m++) {
            element_mult<T, LARGESIZE>(size1_, xtemp, fcc,
                                       x0 + m * element_size_);
            element_smul<T, LARGESIZE>(size1_, xtemp, weight);
            element_set<T, LARGESIZE>(size1_, x0 + m * element_size_, xtemp);
        }
        x0 = data_ + (tstp_ + 1) * element_size_;
        for (m = 0; m <= ntau_; m++) {
            element_mult<T, LARGESIZE>(size1_, xtemp, fcc,
                                       x0 + m * element_size_);
            element_smul<T, LARGESIZE>(size1_, xtemp, weight);
            element_set<T, LARGESIZE>(size1_, x0 + m * element_size_, xtemp);
        }
        x0 = data_ + (tstp_ + 1 + ntau_ + 1) * element_size_;
        for (m = 0; m <= tstp; m++) {
            element_conj<T, LARGESIZE>(size1_, fcc, ft.ptr(m));
            element_mult<T, LARGESIZE>(size1_, xtemp, fcc,
                                       x0 + m * element_size_);
            element_smul<T, LARGESIZE>(size1_, xtemp, weight);
            element_set<T, LARGESIZE>(size1_, x0 + m * element_size_, xtemp);
        }
    }
    delete[] xtemp;
    delete[] fcc;
}

Here is the call graph for this function: