<div dir="ltr">Hello everyone,<div>Considering the system of linear equations in the form Az=b. The term Az is calculated in a matrix-free manner from a(u_h, phi_i) in the alpha_volume() method. </div><div>Is it correct to apply Az by direct matrix multiplication of element stiffness matrix(Ae) and element vector ze like: </div><div><br></div><div>result_vector = Ae*ze;  // Ae.mv(ze, result_vector);</div><div><br></div><div><div>and then acclimate result_vector to <b>"residual r"</b> as per DoF pattern? </div><div><br></div><div>Thank You,</div><div><br></div><div>Regards,</div><div>Ranjeet</div></div></div>