[Libmesh-users] error compiling 'patch_recovery_error_estimator.C' when enabling complex

[Dirk C. <d.c...@tu...>]

Hi,

after configuring the latest CVS-snapshot with --enable-complex,
compilation fails for  src/numerics/patch_recovery_error_estimator.C

// Compute the gradient on the current patch element
// at the quadrature point
Gradient grad_u_h;

for (unsigned int i=3D0; i<n_dofs; i++)
  grad_u_h +=3D dphi[i][qp]*system.current_solution(dof_indices[i]);

There's a type mismatch between dphi and system.current_solution

   PatchRecoveryErrorEstimator::estimate_error(const SteadySystem&,
   std::vector<float, std::allocator<float> >&)':
src/numerics/patch_recovery_error_estimator.C:172: error: no match for
   'operator*' in '(+(+dphi)->std::vector<_Tp, _Alloc>::operator[] [with _T=
p =3D
   std::vector<RealGradient, std::allocator<RealGradient> >, _Alloc =3D
   std::allocator<std::vector<RealGradient, std::allocator<RealGradient> >
   >](i))->std::vector<_Tp, _Alloc>::operator[] [with _Tp =3D RealGradient,
   _Alloc =3D std::allocator<RealGradient>](qp) *
   System::current_solution(unsigned int)
   const((&dof_indices)->std::vector<_Tp, _Alloc>::operator[] [with _Tp =3D
   unsigned int, _Alloc =3D std::allocator<unsigned int>](i))'
/home/dc/src/libmesh.gnu/include/numerics/type_vector.h:491: error: candida=
tes
   are: TypeVector<T> TypeVector<T>::operator*(T) const [with T =3D Real]


Dirk

 --=20
   Dirk Clasen
   Technical University Braunschweig
   Institute of Applied Mechanics
   Spielmannstra=DFe 11
   D-38106 Braunschweig
   URL: http://www.infam.tu-braunschweig.de/~dc