Hello everyone I am having issues trying to solve a problem that has a matrix-vector product in its weak formulation. I defined a macro epsilon as follows : macro epsilon(u1, u2) [dx(u1), dy(u2), dy(u1)+dx(u2)] // and there is a term in the weak formulation that looks like this : int2d(Th)( …

Sadly, I don’t think it’s quite possible. However, you can trick the varf with a hand-defined macro. mesh Th; real[int, int] C(2,2); macro Cdense [[C(0,0), C(0,1)],[C(1,0), C(1,1)]]// varf vPb(u, v)= int2d(Th)( [2*v, v]'*Cdense*[3*u, u] );

Matrix-vector multiplication in weak formulation

General Discussion

prj July 28, 2020, 3:30pm 25

No, I mean that your script is wrong, I can’t even launch it on my machine.

Topic		Replies	Views
Call PETSc matrix operator MatDiagonalSet (…) General Discussion	38	1774	June 3, 2024
Poisson equation with lagrange multiplier using PETSc General Discussion	5	1011	February 27, 2023
Lagrange multiplier for Dirichlet boundary condition General Discussion	30	1055	June 21, 2023
Two Different PDE solvers give different answers General Discussion	17	1244	January 12, 2021
Oldroyd-B model freefem++ code bugs and issues: Assertion fail General Discussion	7	169	February 28, 2025

Matrix-vector multiplication in weak formulation

Related topics