:
>
> On 01/30/2018 07:53 AM, Dulcimer0909 wrote:
> >
> > for(unsigned int i = 0; i < dofs_per_cell; i++)
> > {
> > for (unsigned int j = 0; j < dofs_per_cell; j++)
> >
Hello,
Perhaps this post is a repetition of a problem that I am facing. I am
trying to replace the use of laplace and mass matrices with cell wise
assembly as seen in step-4 to know if I have understood things properly,
which I feel, since I am posting this problem I have not.
What does the fo
atrix_norm_square (solution_v) +
laplace_matrix.matrix_norm_square (solution_u)) / 2
<< std::endl;
old_solution_u = solution_u;
old_solution_v = solution_v;
}
}
}
Thank you
On Sunday, January 28, 2018 at 4:31:17 PM UTC+1, Dulci
ion_v) +
laplace_matrix.matrix_norm_square (solution_u)) / 2
<< std::endl;
old_solution_u = solution_u;
old_solution_v = solution_v;
}
}
}
Dulcimer
On Tuesday, January 23, 2018 at 8:27:37 PM UTC+1, Wolfgang Bangerth wrote:
&
Hello all,
An additional question regarding this thread:
If I do go ahead and replace code, so that it does a cell by cell assembly,
I am a bit lost on how I would store the old_solution (U^(n-1)) for each
cell and retrieve it during the assembly for the Rhs.
grateful if anyone could help.
Du