[petsc-users] question about the most effective implementation

[Michael Povolotskyi]

Dear Petsc team,

what would be the most efficient way to compute the following:

V += A*B*C,

where V is a dense square matrix,

B is a square diagonal matrix, A and C are dense rectangular matrices?

I would like to benefit from the fact that B is diagonal.

thank you,

Michael.