Hello Alexander,
no, you have a Vector<double> as destination object. Thus, you have to
save a scalar, not a vector.
Best Regards,
Markus
On 26.04.2012 08:59, Alexander Grayver wrote:
Hallo Markus,
Orientation vector is defined based on experiment you want to reproduce.
For instance, modeling any finite length dipole of specified power can
be done as:
/*f*(*p*) = a * *n* * \delta (*p* - p_s);/
*n* = {n_x n_y n_z} -- orientation vector;
a = I * L; I -- current, L -- length
Same can be done for magnetic dipoles (antennas).
So this is quite convenient and many real sources can be well
approximated by this vector function.
However, I'm not sure now that the dot product of orientation vector
and shape function is correct representation.
The dot product gives us scalar, whereas we need vector with each
component assigned to individual field components of FE_Nedelec. Is
this still correct?
Thanks.
On 25.04.2012 18:32, Markus Bürg wrote:
Hello Alexander,
then you are right and inserting the dot product of orientation
vector and shape function should work. Where do you get the
orientation vector from? Should this be another argument of the function?
Best Regards,
Markus
On 25.04.2012 15:21, Alexander Grayver wrote:
Markus,
My previous email wasn't clear.
For EM, we can have some vector that defines orientation of the unit
dipole, then source is defined as following
/f(p) = n * \delta (p - p_s);/ n = {n_x n_y n_z} -- orientation vector
Elements of this vector should then be assigned to corresponding
cell's components (i.e. E{x,y,z}).
For instance, x-directed unit dipole would give us 3 component
vector with only 1st element == 1.
--
Regards,
Alexander
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
