auto op_b10 = linear_operator(a.block(1, 0));
auto op_b11 = linear_operator(a.block(1, 1));
- std::array<std::array<decltype(op_b00), 2>, 2> temp{
- {op_b00, op_b01, op_b10, op_b11}};
+ std::array<std::array<decltype(op_b00), 2>, 2> temp
+ {
+ {op_b00, op_b01, op_b10, op_b11}
+ };
auto op_b = block_operator<2, 2, BlockVector<double>>(temp);
{
// And finally complicated block structures:
- std::array<std::array<decltype(op_b00), 3>, 3> temp2{
- {op_b00, op_b01, op_b00, op_b10, op_b11, op_b10, op_b10, op_b11, op_b10}};
+ std::array<std::array<decltype(op_b00), 3>, 3> temp2
+ {
+ {op_b00, op_b01, op_b00, op_b10, op_b11, op_b10, op_b10, op_b11, op_b10}
+ };
auto op_upp_x_upu = block_operator<3, 3, BlockVector<double>>(temp2);
op_upp_x_upu.reinit_domain_vector(u, false);
op_upp_x_upu.vmult_add(v, u);
PRINTME("v", v);
- std::array<std::array<decltype(op_b01), 1>, 3> temp3{
- {op_b01, op_b11, op_b11}};
+ std::array<std::array<decltype(op_b01), 1>, 3> temp3
+ {
+ {op_b01, op_b11, op_b11}
+ };
auto op_upp_x_p = block_operator<3, 1, BlockVector<double>>(temp3);
- std::array<std::array<decltype(op_b01), 3>, 1> temp4{
- {op_b00, op_b01, op_b00}};
+ std::array<std::array<decltype(op_b01), 3>, 1> temp4
+ {
+ {op_b00, op_b01, op_b00}
+ };
auto op_u_x_upu = block_operator<1, 3, BlockVector<double>>(temp4);
auto op_long = op_u_x_upu * transpose_operator(op_upp_x_upu) * op_upp_x_p;
auto op_b1 = linear_operator(a.block(1, 1));
auto op_b2 = linear_operator(a.block(2, 2));
- std::array<decltype(op_b0), 3> temp{{op_b0, op_b1, op_b2}};
+ std::array<decltype(op_b0), 3> temp {{op_b0, op_b1, op_b2}};
auto op_b = block_diagonal_operator<3, BlockVector<double>>(temp);
// And finally the other block_diagonal_operator variant:
- std::array<decltype(op_b0), 5> temp2{{op_b0, op_b0, op_b0, op_b0, op_b0}};
+ std::array<decltype(op_b0), 5> temp2 {{op_b0, op_b0, op_b0, op_b0, op_b0}};
auto op_c = block_diagonal_operator<5, BlockVector<double>>(temp2);
auto op_d = block_diagonal_operator<5, BlockVector<double>>(op_b0);