On Monday January 16 2012 19:44:23 Marie E. Rognes wrote: > Diff works, thanks! > > While stumbling around with this, I'm encountering > some more stuff: > > For instance, with this: > -- > from dolfin import * > mesh = UnitSquare(2, 2) > V = VectorFunctionSpace(mesh, "CG", 2) > Q = FunctionSpace(mesh, "CG", 1) > W = V*Q > > W0 = W.sub(0) > print W0.component() > -- > > I get this: > > <Swig Object of type 'std::vector< dolfin::uint > *' at > 0x36641e0>swig/python detected a memory leak of type 'std::vector< > dolfin::uint > *', no destructor found. > > (1) Is the memory leak supposed to be there?
No. > (2) Could we perhaps also have an out typemap for the std::vector? Yes, will add out typemap macro for these guys. Johan > -- > Marie > > On 01/16/2012 06:36 PM, Johan Hake wrote: > > Apply this diff and it should work: > > > > ----------------------------------------------- > > > > === modified file 'dolfin/swig/std_vector_typemaps.i' > > --- dolfin/swig/std_vector_typemaps.i 2011-10-25 17:01:18 +0000 > > +++ dolfin/swig/std_vector_typemaps.i 2012-01-16 17:12:09 +0000 > > @@ -17,7 +17,7 @@ > > > > // along with DOLFIN. If not, see<http://www.gnu.org/licenses/>. > > // > > // First added: 2009-08-31 > > > > -// Last changed: 2011-10-09 > > +// Last changed: 2012-01-16 > > > > //===================================================================== > > ======== // In this file we declare what types that should be able to > > be passed using > > > > a > > @@ -464,8 +464,9 @@ > > > > ARGOUT_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(dolfin::uint, INT32, columns, > > > > NPY_INT) > > > > ARGOUT_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(double, DOUBLE, values, > > NPY_DOUBLE) > > > > -IN_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(double, DOUBLE, values, NPY_DOUBLE, > > double, d) > > -IN_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(dolfin::uint, INT32, columns, > > NPY_UINT, uint, I) > > +IN_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(double, DOUBLE, , NPY_DOUBLE, > > double, d) +IN_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(int, INT32, , NPY_INT, > > int, i) +IN_TYPEMAP_STD_VECTOR_OF_PRIMITIVES(dolfin::uint, INT32, , > > NPY_UINT, uint, I) > > > > PY_SEQUENCE_OF_SCALARS_TO_VECTOR_OF_PRIMITIVES(dolfin::uint, INT32, > > > > coloring_type, uint, -1) > > > > PY_SEQUENCE_OF_SCALARS_TO_VECTOR_OF_PRIMITIVES(dolfin::uint, INT32, > > > > value_shape, uint, -1) > > > > --------------------------------------------------------- > > But the returned SubSpace is not a dolfin.FunctionSpace. Maybe you can > > add that to the site-packages/dolfin/function/functionspace.py? > > > > Johan > > > > On Monday January 16 2012 17:47:36 Marie E. Rognes wrote: > >> How do one extract a sub^n-space using a component tuple? > >> I've tried various versions of this: > >> > >> import numpy > >> from dolfin import * > >> mesh = UnitSquare(2, 2) > >> V = VectorFunctionSpace(mesh, "CG", 2) > >> Q = FunctionSpace(mesh, "CG", 1) > >> W = V*Q > >> component = numpy.array((0, 0)) > >> W.extract_sub_space(component) > >> > >> but generally get this: > >> > >> TypeError: in method 'FunctionSpace_extract_sub_space', argument 2 of > >> type 'std::vector< dolfin::uint> const&' > >> > >> > >> -- > >> Marie > >> > >> _______________________________________________ > >> Mailing list: https://launchpad.net/~dolfin > >> Post to : dolfin@lists.launchpad.net > >> Unsubscribe : https://launchpad.net/~dolfin > >> More help : https://help.launchpad.net/ListHelp _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@lists.launchpad.net Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
