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" <simon.am...@mpsa.com> To: "getfem-users" <getfem-users@nongnu.org> 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
