[Dune] [#522] Jacobian type is FieldVector< FieldVector< K,m> >

Dune dune at dune-project.org
Fri Jun 26 12:14:58 CEST 2009


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

The following task has a new comment added:

FS#522 - Jacobian type is FieldVector< FieldVector< K,m> >
User who did this - Uli Sack (usack)

----------
In my opinion suggestion 2) is prefereable. Firstly because FieldVector is not supposed to be nested in the first place, and secondly because a FieldVector<FieldVector,double,N>,M> simply doesn't work as expected. A very simple example:

Dune::FieldVector<Dune::FieldVector<double,3>,1> a(Dune::FieldVector<double,3>(0.0)),b;
b=a;
a==b;

Compiler error:
test.cc: In function ‘int main()’:
test.cc:41: error: no match for ‘operator==’ in ‘a == b’

Lastly, if full jacobian matrices are supported (meaning for vector valued basis functions) I expect them to behave like matrices including linear algebra functionality. 

Therefore I agree with Carsten and support a switch to FieldMatrix in the JacobianType for the LocalBasisTraits.
----------

More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=522#comment989

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