On 4/25/19 7:56 AM, Gabriel Peters wrote:
> I wanted solution_o to be vector-valued.
> And I want to multiply it with a gradient function. So I want
>
> local_matrix(i,j) += solution_o * fe_values[velocities].shape_gradient(i,q)
> *
> fe_values[velocities].shape_gradient(j,q);
>
> Does this
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 25.04.19 um 15:28 schrieb Wolfgang Bangerth
> On 4/25/19 4:51 AM, gabriel.pet...@koeln.de wrote:
> >
> > std::vector > solution_old_values(n_q_points);
> >
> > While assembling the system over all
On 4/25/19 4:51 AM, gabriel.pet...@koeln.de wrote:
>
> std::vector > 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
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
