The following minimal test for Function in c++ reveals some bugs. I
guess this example can help me with dealing with the current issues of
ImageFunction.
(1) interpolate.py does not work when a Function is created in c++ and
wrapped (see comment [2]). It seems that the bug is originated from
the copy constructor (see comment [3])
(2) In order to perform the interpolation, why is it necessary to
create another Function and then copy it?
(3) Signature checkes seem not working properly (see comment [1]). The
signature-based assignments are error-prone anyway, why the
object-oriented approach is not used?
-Ali
-----------------------
// FunctionTest.cpp
#include <dolfin/function/Function.h>
#include <dolfin/function/FunctionSpace.h>
#include <dolfin/mesh/Mesh.h>
#include <dolfin/fem/FiniteElement.h>
#include <dolfin/fem/DofMap.h>
#include <dolfin/mesh/UnitSquare.h>
#include <boost/shared_ptr.hpp>
using namespace dolfin;
int main()
{
typedef boost::shared_ptr< const Mesh > MeshConstPointerType;
typedef boost::shared_ptr< const FiniteElement >
ElementConstPointerType;
typedef boost::shared_ptr< const DofMap > DofMapConstPointerType;
std::cout << "- Creating mesh...";
int msize[2];
msize[0] = 9; msize[1] = 9;
UnitSquare mesh(msize[0], msize[1]);
std::cout << "done." << std::endl;
std::string elemSig("FiniteElement('Lagrange', 'triangle', 1)");
// std::string elemSig("dummy"); // [1] This gives segmentation fault.
std::string dofSig("FFC dof map for FiniteElement('Lagrange',
'triangle', 1)");
std::cout << "- Creating FunctionSpace...";
FunctionSpace V(MeshConstPointerType(&mesh, NoDeleter<const
dolfin::Mesh>()),
ElementConstPointerType(new
dolfin::FiniteElement(elemSig)),
DofMapConstPointerType(new
dolfin::DofMap(dofSig, mesh)));
std::cout << "done." << std::endl;
std::cout << "- Creating Function...";
Function v(V);
std::cout << "done." << std::endl;
/*
// [2] This gives segmentation fault.
std::cout << "- Interpolating as interpolate.py...";
Function Pv(v); // [3] The copy ctor doesn't work.
v.interpolate(Pv.vector(), v.function_space());
std::cout << "done." << std::endl;
*/
// [4] This seems to work fine.
std::cout << "- Interpolating...";
v.interpolate(v.vector(), v.function_space());
std::cout << "done." << std::endl;
};
-----------------------
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
