Hi, when running DOLFIN in parallel, data reduction would be a nice thing to have. Typically, it would be good to have something like this:
mesh = UnitSquare(2000, 2000) reduced_mesh = UnitSquare(200, 200) V = FunctionSpace(mesh, "CG", 1) reduced_V = FunctionSpace(reduced_mesh, "CG", 1) u = Function(V) ... reduced_u = interpolate(u, reduced_V) f << reduced_u However, the problem is that the automatic partitioning of the meshes constructs two different non-overlapping meshes on each local node. I really don't want to allow extrapolation in this case. I would be happy to contribute fixing this problem; does anybody see how to attack it? Ola -- Ola Skavhaug Research Programmer Simula Research Laboratory _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@lists.launchpad.net Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
