dolfin team mailing list archive

[Ola Skavhaug <skavhaug@xxxxxxxxx>]

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