[deal.II] Re: "operator+=" error in the assemble_system()

Dear Jean-Paul, > Firstly (and just so that we're all clear on what you're trying to > implement), was I correct that your comment that I linked contained the > description of the integral that you are trying to compute? > > Yes it is. that is something I have to do. > No, this is not true

[deal.II] Re: "operator+=" error in the assemble_system()

Dear Kyusik, Firstly (and just so that we're all clear on what you're trying to implement), was I correct that your comment that I linked contained the description of the integral that you are trying to compute? > As far as I know fe_values2.shape_value(i,q_index) is vector-valued shape > fu

[deal.II] Re: "operator+=" error in the assemble_system()

Dear Bruuno and Jean-Paul First of all thank you for your reply I'm sorry, I should have been more specific about my question As far as I know fe_values2.shape_value(i,q_index) is vector-valued shape function(fe2 is distributed by "dof_handler2.distribute_dofs (fe2)") and sol_grad is also

[deal.II] Re: "operator+=" error in the assemble_system()

Dear Kyusik, I gather that what you are trying to implement here is the equation that you present in this post . I concur with Bruno's assessment of the problem - unfortunately your implementation is incorrect (compare what you'v

[deal.II] Re: "operator+=" error in the assemble_system()

Kyusik, On Friday, May 19, 2017 at 8:27:39 AM UTC-4, hanks0...@gmail.com wrote: > > > std::vector > sol_grad(n_q_points); > > cell_rhs(i) += fe_values2.shape_value(i,q_index) * > sol_grad[q_index] *(F/p(0)/grad_sol_sq)* > fe_