construct_m_matrix Subroutine

public subroutine construct_m_matrix(p, f, du, m)

this subroutine returns the matrix $$M_{ij c}= \sum_o^{\text{occ orbitals}} \sum_{k} \frac{d u_{ko} ({\bf R}) }{d {\bf R}} \left( 4 (i j | k o) - (i k | j o) - (i o | j k) \right)$$ it expects to find the MO-space twoelectron integrals in F%teint

Arguments

Type	Intent	Optional	Attributes		Name
type(param),	intent(in)			::	p
type(wf_data),	intent(in)			::	f
real(kind=long),	intent(in),		dimension(:, :, :)	::	du
real(kind=long),	intent(out),		dimension(:, :, :)	::	m