Using vectors, the dot function provides bad results in // while results are
good when using element_type
For example if Xh is the function space :
auto v1 = backend->newVector( Xh );
auto v2 = backend->newVector( Xh );
[do operations ....]
v1->close(); v2->close();
dot( *v1 , *v2 ); // wrong result
element_ptrtype ev1 ( new element_type ( Xh ) );
element_ptrtype ev2 ( new element_type ( Xh ) );
*ev1 = *v1;
*ev2 = *v2;
dot( *ev1 , *ev2 );// good result (and same result than in sequential)
Original issue reported on code.google.com by stephane...@gmail.com on 30 Aug 2012 at 7:53
Original issue reported on code.google.com by
stephane...@gmail.com
on 30 Aug 2012 at 7:53