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] );

Dear prj I have implemented my complete problem using TAO, but the due to I can not find a correct optimizing function J for my DJ and HJ, almost each Newton iteration needs very long CG iteration computing J and DJ and can not reach convergence. In my original IPOPT program, I turned off the line…

I’m not a TAO expert, you’ll have to try different methods. Also, I still don’t know what you are trying to solve exactly, so it’s difficult for me to give you precise pointers. SNESSolve either solves F(x) = 0, or F(x) = b. You are free to put what you want in the output vector, but this acts as an…

Thank you for your reply, prj. Here I uploaded my problem using SNESSolve, all the macro and constants are deleted to reduce the program from 2300 lines to 270+ lines. This is a nonlinear problem F(x)=0, the F is given in funcRes part and the correct Jacobian matrix is obtained in funcJ part. Thei…

Your script doesn’t even run.

Sorry, I don’t have a simplified version that can run, sorry.

Is the way of adding boundary condition in SNESSolve the same with TAO? If the initial value of boundary nodes are set to meet the requirements of boundary condition, can it be maintained using the constrain of upper and lower bound or using the form as on(1,2,3,4, du1=0,du2=0,du3=0) in varf functi…

The boundary conditions have to be set in your initial guess.

Thans for your timely reply. I have found that in the definition of funcRes and funcJ, “tgv” is used in the vector and matrix forming such as “Loc = vPb(Wh, Wh, tgv = -2)” in SNESSolve. The tgv in FreeFem++ is always represent 1e30 to implement Dirichlet boundary conditions and in most situation c…

I always advise against using the default value of 10e+30. -1 (or -2 if your problem has homogeneous boundary conditions and is symmetric) is always better when you do iterative algorithms, such as a Newton loop in SNESSolve/TaoSolve.

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	2082	June 3, 2024
Poisson equation with lagrange multiplier using PETSc General Discussion	5	1061	February 27, 2023
Lagrange multiplier for Dirichlet boundary condition General Discussion	30	1233	June 21, 2023
Two Different PDE solvers give different answers General Discussion	17	1326	January 12, 2021
Oldroyd-B model freefem++ code bugs and issues: Assertion fail General Discussion	7	249	February 28, 2025

Matrix-vector multiplication in weak formulation

Related topics