Dear Kyusik, I gather that what you are trying to implement here is the equation that you present in this post <https://groups.google.com/d/msg/dealii/7NaagkXrevs/YAMg3QPBAgAJ>. I concur with Bruno's assessment of the problem - unfortunately your implementation is incorrect (compare what you've written as contributions to the cell_rhs to the equation you listed), so you probably want to think a little more carefully about what you're doing here.
Regards, Jean-Paul On Friday, May 19, 2017 at 2:33:18 PM UTC+2, Bruno Turcksin wrote: > > Kyusik, > > On Friday, May 19, 2017 at 8:27:39 AM UTC-4, hanks0...@gmail.com wrote: >> >> >> std::vector<Tensor<1,dim> > 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_values2.JxW(q_index); >> > > >> /user2/hanks318/dealii/dealii_pre/examples/step-6/step-6.cc:665: error: >> no match for ‘operator+=’ in ‘cell_rhs.dealii::Vector<number>::operator() >> [with Number = double](i) += dealii::operator*(const dealii::Tensor<1, dim, >> double>&, double) [with int dim = 2](fe_values2.dealii::FEValues<2, >> 2>::<anonymous>.dealii::FEValuesBase<dim, spacedim>::JxW [with int dim = 2, >> int spacedim = 2](q_index))’ >> > sol_grad[q_index] is a Tensor<1,dim,double> so it has dim components. You > cannot add it to a double. > > Best, > > Bruno > -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
