[dune-pdelab] Implementation of alpha_volume_method by direct multiplication of Az

[Ranjeet kumar]

Hello everyone,
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.
Is it correct to apply Az by direct matrix multiplication of element
stiffness matrix(Ae) and element vector ze like:

result_vector = Ae*ze;  // Ae.mv(ze, result_vector);

and then acclimate result_vector to *"residual r"* as per DoF pattern?

Thank You,

Regards,
Ranjeet
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20190910/b1e95e49/attachment.htm>