(Deleted) Dirichlet conditions at corners of one simplex
3 months ago by
from dolfin import * mesh = UnitSquareMesh(10, 10) V = FunctionSpace(mesh, 'CG', 2) u = TrialFunction(V) v = TestFunction(V) A = assemble(dot(grad(u), grad(v))*dx)
I now need to take any cell, let's say `cell`, and replace the lines corresponding to the test functions at the corners of the cell by "Dirichlet" lines, i.e., zero all over and 1 at the node index `i`. (I wouldn't mind extending this to all test functions with nonvanishing support in that cell.)
Community: FEniCS Project
Please login to add an answer/comment or follow this question.
The thread is closed. No new answer/comment may be added.