Hi Bruno, Thank you for the reference.
>> If I understand correctly, your right-hand-side is similar to a Dirac delta. So when you write your weak form, your right-hand-side is something like \int \int \phi(x,y) >> \delta(x_0,y_0) dx dy = \phi(x_0,y_0) So to build your right-hand-side, you just need to evaluate the test functions at one point. Yes, as mentioned in the cited thread with the dirac delta function, my \phi(x,y) is just a constant. So my question is, shouldn't the value on the graph attached show a value of 1 without multiplying this constant \phi value? The function is as follows: template <int dim> void AdvectionProblem<dim>::create_point_source() { Point<dim> sourcePoint; for(unsigned int i = 0; i < dim; ++i ) { sourcePoint(i) = 0.6; } VectorTools::create_point_source_vector(dof_handler, sourcePoint, right_hand_side); // right_hand_side /= 800; //80; //10; } //==================== // In all static functions //==================== beta(0) = 1.; beta(1) = 1.; //==================== // function name : BoundaryValues<dim>::value_list(...) //==================== values[i] = 0; On Friday, May 12, 2017 at 2:19:05 PM UTC+2, Bruno Turcksin wrote: > > Kartik, > > On Thursday, May 11, 2017 at 7:05:35 PM UTC-4, Kartik Jujare wrote: >> >> I added a create_point_source function after the assemble function in >> step 12. I have also changed the beta value to (1, 1) and g=0 everywhere on >> the boundary. Although the solution is similar to what I expected. I am not >> able to understand why the magnitude is scaled by an arbitrary number(or >> are they arbitrary?), in my case 800. Inside this >> create_point_source_vector function, only the shape values get added to the >> rhs then where does this magnitude come from and how can one control it? >> Would be grateful to anyone who can shed some light on this. >> > We need to see the function but you may want to look here for example on > how to do a point source > https://groups.google.com/forum/#!searchin/dealii/dirac$20delta%7Csort:relevance/dealii/HTcSxhnzc5s/jDTMErSICAAJ > > 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.