def METHOD_NAME(self,beta1,beta2,k,index,Cdtype):
    A = self.space.DDG_stiff_matrix(self.qf, self.pde.diffusion_coefficient)
    P,S,H = self.space.edge_matrix(self.IntervalQuadrature, self.pde.diffusion_coefficient,index,Cdtype)
    AD = A-S-S.T+beta1*P + beta2*H
    return AD