Hi Wolfgang, thank you for your reply.
I've had another look at the functions available and i'm stuck on how to
implement these in this case.
The whole point of not imposing zero for all components in a Dirichlet
sense is that I am testing for a case where I have inhomogeneous normal
Hi Wolfgang
Thanks for your reply. I carefully tested these two methods and found that
they produce two different rhs vectors given identical previous solution.
There is no non-zero boundary values in my case, and the constrained values
are indeed zero. It seems that the unconstrained
On 11/28/2017 02:33 PM, Jie Cheng wrote:
Because the global stiffness matrix is the sum of local matrices, ideally this
approach should work as the first one. But the solution turned out to be
garbage. Could anybody see why the second approach is wrong?
Can you be more specific what