convolution_omp() [1/2]

template<typename T >

void cntr::convolution_omp	(	int	omp_num_threads,
		herm_matrix< T > &	C,
		herm_matrix< T > &	A,
		herm_matrix< T > &	Acc,
		function< T > &	ft,
		herm_matrix< T > &	B,
		herm_matrix< T > &	Bcc,
		T	beta,
		T	h,
		int	SolveOrder
	)

Returns convolution \(C = A\ast f \ast B\) for all timesteps.

Purpose

Computes contour convolution \(C = A\ast f \ast B\) of the objects with class 'herm_matrix' for all time steps. Works for a scalar and square matrices. openMP parallelized version.

Parameters

omp_num_threads	[int] The number of openMP threads to be used. Set to the number of threads in the current team.
C	[herm_matrix] Matrix to which the result of the convolution on Matsubara axis is given
A	[herm_matrix] contour Green's function
Acc	[herm_matrix] hermitian conjugate of A
ft	[function] function \(f\)
B	[herm_matrix] contour Green's function
Bcc	[herm_matrix] hermitian conjugate of B
beta	inversed temperature
h	time step interval
SolveOrder	[int] integrator order

Definition at line 5548 of file cntr_convolution_impl.hpp.

                                                                       {
    int tstp;
    for (tstp = -1; tstp <= C.nt(); tstp++)
        convolution_timestep_omp<T>(omp_num_threads, tstp, C, A, Acc, ft, B, Bcc, beta,
                                    h, SolveOrder);
}