Commits

Lisandro Dalcin committed fe0c3cd

Fix boundary conditions in Elasticity example

  • Participants
  • Parent commits 9a4e911

Comments (0)

Files changed (1)

File demo/Elasticity.c

   for (i=0; i<dim; i++) {ierr = IGASetFieldName(iga,i,fieldname[i]);CHKERRQ(ierr);}
 
   // Set boundary conditions
+  for (i=0; i<dim; i++) {ierr = IGASetBoundaryValue(iga,0,0,i,0.0);CHKERRQ(ierr);}     // Dirichlet
   for (i=0; i<dim; i++) {ierr = IGASetBoundaryLoad (iga,0,1,i,load[i]);CHKERRQ(ierr);} // Neumann
 
   // Create linear system