Hello Alexander,
I am not sure about that. How is a vector-valued delta-function defined?
Best Regards,
Markus
On 25.04.2012 14:24, Alexander Grayver wrote:
Hello,
I realized that this function works only for scalar FE. Looking at the
implementation it seems the only thing I need to change when
considering vector FE is this:
for(unsigned int i=0; i<dofs_per_cell; i++)
rhs_vector(local_dof_indices[i]) = fe_values.shape_value(i,0);
Here, instead of /shape_value/ I should use /shape_value_component/
for each component and then put dot product to the rhs_vector?
Thanks.
On 17.04.2012 08:16, Alexander Grayver wrote:
I found VectorTools::create_point_source_vector
This should be what I need. I guess I can scale this vector by any
scalar /a/ and obtain /f(p) = a * \delta (p - p_s);
/Thanks/.
/
On 15.04.2012 14:50, Alexander Grayver wrote:
Hello,
What is the best way to represent delta source in deal.II?
More precisely I want my RHS to be
/f(p) = a * \delta (p - p_s);/
/f(p) != 0/ iif /|p - p_s| = 0/;
Doing this properly would mean that I have to define element within
this function happens not to be zero and then distribute
contribution among its q points?
--
Regards,
Alexander
--
Regards,
Alexander
_______________________________________________
dealii mailing listhttp://poisson.dealii.org/mailman/listinfo/dealii
--
Regards,
Alexander
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
