Can anyone comment on the extraction of the linear part in varf() ? There is an issue of sign, it seems that it is treated differently depending if it is an integral source or an on() boundary condition. The simple code int n=3; mesh Th=square(n,n); //plot(Th); real a=1.,b=1.; fespace Vh(Th,P1);…

solve supposes the linear and bilinear form represent the LHS and RHS of an equation system. varf defines a linear or bilinear form that is not necessarily part of an equation. Since you impose inhomogeneous BC, the sign flip in the latter is changing the BC, not just the forcing. varf test(u,v) =…

Thank you. It is quite puzzling that solve and varf do not use the same conventions. Additionally in varf the convention for int2d or int1d is opposite to that of on(). Then I have to remind of it when changing from one formulation to the other. Setting u=-b instead of u=b is very much counterintuit…

Additionally in varf the convention for int2d or int1d is opposite to that of on() You are thinking that there is a convention, but there is none. int2d or int1d are used to discretise linear or bilinear form (which may be related or not, as Chris highlighted). on() is used to impose Dirichlet bo…

Hi Chris, I noticed that A and tab are constructed from same varf, but A won’t include the linear part a*v, right?

To complete the post of Chris, a third option is to use varf test(u,v) =int2d(Th)(dx(u)*dx(v)+dy(u)*dy(v)) -int2d(Th)(a*v) +int1d(Th)(tgv*u*v) -int1d(Th)(tgv*b*v) ; This formulation is valid for solve and varf without change (the latter must use the line tab=-tab;). This formulation avoids the use…

Why would you want to avoid the use of on()? That allows more flexibility in terms of how Dirichlet BC’s are enforced via penalization tgv>0, or elimination (asymmetric tgv=-1 or symmetric tgv=-2). It also allows you to ignore BC with tgv=0. matrix A=test(Vh,Vh,tgv=...); tab=test(0,Vh,tgv=...);

It depends of course of what you want to do.

You are not answering to the question “Why would you want to avoid the use of on() ?”

Behavior of varf with rhs and on()

General Discussion

prj May 14, 2024, 12:25pm 3

Duplicate of Change of sign in linear part of variational problem (varf vs problem).

Topic		Replies	Views
Option "rhsvarf" for on() Feature Request	8	113	July 19, 2024
Change of sign in linear part of variational problem (varf vs problem) General Discussion	1	343	February 22, 2024
The elasticity problem with "varf" General Discussion	23	974	October 19, 2022
Variable formulation error General Discussion	3	84	June 22, 2025
Type mismatch when using varf General Discussion	0	300	August 20, 2021

Behavior of varf with rhs and on()

Related topics