Hello all,

the initialize_info method in MeshWorker::Assembler::Functional is defined as


template<typename number  = double> 
template<int dim> 
void MeshWorker::Assembler::Functional< number >::initialize_info  ( DoFInfo< 
dim > &  info,  


 
 bool  interior_face  

 ) 
 

Shouldn't it be templated with 

    template <typename number = double>
    template <class DOFINFO>
    void initialize_info(DOFINFO& info, bool interior_face);

instead? This would make it consistent with the other similar initialize_info 
methods and get rid of some runtime link errors that I am seeing when 
assembling 
a DG functional with terms defined on inner element faces:

  MeshWorker::DoFInfo<dim> dof_info(dof_handler);
  assembler.initialize_info(dof_info, true);

and get at link time

/home/mihai/MyDocuments/autodiff/adaptive/optimality/code/c++/EllipticOptimalitySystem.cc:1213:
 undefined reference to `void 
dealii::MeshWorker::Assembler::Functional<double>::initialize_info<2>(dealii::MeshWorker::DoFInfo<2,
 2, double>&, bool)'

Or should I use intialize_info in some other way?

Thanks,
 -- Mihai

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to