Hey everyone,
I want to set up a BlockMatrix, with entries, which are enforced by a
velocity components of aBlockVector "solution_old".
The Vector contains velocity entries in the first dim terms and pressure
entries in the last component.
Similar to the 'assemble_system' function of Step-21 tutorial, I use
std::vector<Vector<double> > solution_old_values(n_q_points);
While assembling the system over all cells I call
fe_values.get_function_values (solution_old, solution_old_values);
In the Step-21 tutorial, they call
const double old_s = old_solution_values[q](dim+1);
to get the dim+2 component of the vector.
Up to here, everything compiles fine, but I dont know how to apply the
values correctly
Therefore I have the following two questions:
Do I simply have to call: const double solution_o =
solution_old_values[q](0);
to get the velocity component? In fact they are vector-valued, so a double
doesn't make sense that much.
And thus my second question is: Do I apply the value correctly to the local
matrix
by calling:
local_matrix(i,j) += solution_o * fe_values[velocities].shape_value(i,q) *
fe_values[velocities].shape_value(j,q);
As always, thanks a lot anf best regards
Gabriel
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
