Hello,
#include <base/quadrature_lib.h>
I am trying to use the *1d* specialization of the
KellErrorEstimate::estimate. I've blindly copied the example from
step-17 which does not work in 1d, but modified parts of the code such
that it should; I hope :-)
Question: Why does this not work where dim=1?
KellyErrorEstimator<dim>::estimate (dof_handler,
QGauss<dim>(2),
typename FunctionMap<dim>::type(),
localized_solution,
local_error_per_cell,
std::vector<bool>(),
0,
multithread_info.n_default_threads,
this_mpi_process);
Instead a find an error like this:
============================ Remaking Makefile.dep
==============debug========= step-2.cc
step-2.cc: In member function 'void
SchrodingerProblem<dim>::refine_grid() [with int dim = 1]':
step-2.cc:795: instantiated from 'void SchrodingerProblem<dim>::run()
[with int dim = 1]' step-2.cc:824: instantiated from here
step-2.cc:643: error: no matching function for call to
'dealii::KellyErrorEstimator<1>::estimate(dealii::DoFHandler<1>&,
dealii::QGauss<1>, std::map<unsigned char, const dealii::Function<1>*,
std::less<unsigned char>, std::allocator<std::pair<const unsigned char,
const dealii::Function<1>*> > >, dealii::PETScWrappers::Vector&,
dealii::Vector<float>&, std::vector<bool, std::allocator<bool> >, int,
unsigned int&, const unsigned int&)' make: *** [step-2.g.o] Error 1
I'm fairly sure the function *is* there in quadrature_lib.h
What have I missed here? I'm baffled by this? Has anyone any idea how to
run adaptive meshing for a 1d problem or know what I have overlooked?.
Thanks!
Best,
Toby
--
Toby D. Young - Adiunkt (Assistant Professor)
Department of Computational Science
Institute of Fundamental Technological Research
Polish Academy of Sciences
Room 206, ul. Swietokrzyska 21
00-049 Warszawa, Polska
+48 22 826 12 81 ext. 184
http://rav.ippt.gov.pl/~tyoung
_______________________________________________
