getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Periodic problem


From: Bourriche Lahoussaine
Subject: Periodic problem
Date: Wed, 20 Jul 2022 09:40:24 +0200 (CEST)

Hy;

For periodic condition i do something like :

--------------------------------------------------------------------------------------------------------------------
ConstraintMatrix = gf.Spmat('empty',1 ,mfu.nb_basic_dof())
L=np.zeros((mfu.nb_basic_dof(),1))
   
for i in range(len(leftDof)):
      ConstraintMatrix[0,leftDof[i]]=1
      ConstraintMatrix[0,rightDof[i]]=1
      md.add_variable('mult_spec'+str(i),1)
      # print(i)
      md.add_constraint_with_multipliers( 'u','mult_spec'+str(i), ConstraintMatrix,0)

--------------------------------------------------------------------------------------------------------------------
And it worked perfectly;
Now i want to add another periodic condition: this new on derivatives 
The derivative of u in leftDof  must equal the derivative of T in rightDOF

Do you have ideas



Best regards
Lahoussaine

reply via email to

[Prev in Thread] Current Thread [Next in Thread]