On 11/29/2017 09:39 PM, Jie Cheng wrote:
Here is what I suspect: Method 1 uses the stiffness_matrix that has been
modified during the assembly, while method 2 is using the original
stiffness. Although method 2 will eventually distribute the rhs vector
as well as the matrix, it is distributed
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 component
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 "garbage
Hi everyone
I have a time-dependent linear elasticity solver which is different from
step-17 in which I assume small-deformation all the time. Due to the time
discretization, I have a -K*u^n term at the rhs of the system where u^n is
the known displacement (at previous time step) and K is the s
