[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