[Dune] [#1011] operator* in densevector used as dot product in dune-istl is not suitable for complex vectors
Dune
flyspray at dune-project.org
Thu Jun 21 19:08:41 CEST 2012
THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.
The following task has a new comment added:
FS#1011 - operator* in densevector used as dot product in dune-istl is not suitable for complex vectors
User who did this - Jö Fahlke (joe)
----------
There was a decision about this in the Dune-Dev-Meeting 2012 in Münster (see
http://users.dune-project.org/projects/dev-meeting-2012/wiki/Topics).
In essence:
- We want a free function dot(a,b), where a and b have the same recursive
blocking structure, but may have different field types.
+ This function should work for all meaningful combinations of fundamental
types and complex<fundamental type> [1].
+ This functions should work for scalars, dense vectors and ISTL vectors
[2].
- On the vector classes we want a template method a.dot(b), with the same
limitations for a and b as above.
- The complex conjugate should be done in the same argument as it is done by
matlab [3].
[1] This requires promotion traits to determine the return type. Markus had
implemented promotion traits (also in the context of complex) at some
point, but had to revert because of high-level problems in PDELab. This
could however be used as a starting point.
[2] The destinction between vector on one hand and scalars and complex on the
other hand can be done by checking whether the type A of the argument has
a member type A::FieldType.
[3] We found that pythons dot() is simply operator*(). Matlab does have dot()
with complex conjugate, but I don't remember which of the two arguments is
conjugated.
----------
More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=1011#comment3745
You are receiving this message because you have requested it from the Flyspray bugtracking system. If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above.
More information about the Dune
mailing list