Re: [Getfem-users] Sliding angle for Dirichlet condition

[Yves Renard]

Dear Simon,

A priori, H is not necessarilly a hat matrix. In that case, you can juste take 
an orthognal vector, for instance [-2 1] and set
H = [-2 1;0 0]], this should work, I think.

Best regards,

Yves


----- Original Message -----
From: "SIMON AMEYE" <address@hidden>
To: "getfem-users" <address@hidden>
Sent: Tuesday, October 30, 2018 1:52:54 PM
Subject: [Getfem-users] Sliding angle for Dirichlet condition

C1-Non sensitive

________________________________
Hi all,

It has been a long time I have tried to understand the H matrix for Dirichlet 
conditions.
I work with a 2D mesh, and I would like my boundary to slide freely according 
to a vector, let's say [1 2].
I think I understood that H needs to be a Hat matrix.
How to do so ?

I tried to construct a hat matrix this way, but the "Angle" is not respected :

x4 = 1;
x1 = x4/(pi/2/Angle);
k = x4+x1;
x2 = sqrt(x1)*(-sqrt(k-x1));
x3 = x2;
gf_model_set(md, 'add initialized data', 'H_TOP', [x1 x2;x3 x4]);
gf_model_set(md, 'add initialized data', 'VECTOR_TOP', [0;0]);
gf_model_set(md, 'add generalized Dirichlet condition with multipliers', mim, 
'u', mfu, TOP_BOUND,'VECTOR_TOP', 'H_TOP');

Simon


[cid:image001.png@01D193D3.0D03A420]



SIMON AMEYE
DQI/DRIA/DSTF/SEPC
Apprenti IFP
________________________________
CENTRE TECHNIQUE VELIZY A