Hi Lawrence,
no worries, I can work with the 2+1 case, which is the really interesting case of course, just wanted to check.
Btw. the \hat{z} field works now, and I made some progress with the mixed solver, implemented the mass matrix \tilde{M}_u=M_u+\omega_N^2 Q M_b^{-1} Q and operators H and \hat{H} for the Helmholtz operator and multigrid preconditioner, so everything is there to at least test the pressure solve, which is the hard bit anyway. I guess everything that works in 2d will also work in 2+1d, so it should be a relatively straighforward rewrite of what I’ve got so far.
Thanks,
Eike
--
Dr Eike Hermann Mueller
Research Associate (PostDoc)
Department of Mathematical Sciences
University of Bath
Bath BA2 7AY, United Kingdom
+44 1225 38 5803
e.mueller@bath.ac.ukhttp://people.bath.ac.uk/em459/
On 7 Dec 2014, at 17:49, Eike Mueller <e.mueller@bath.ac.uk> wrote:Dear firedrakers,
do the hierarchical meshes and function spaces currently only work in 2+1 dimension and not for 1+1? If I run the code below it works for dimension=3, but if I replace this by dimension=2 it crashes with
[0]PETSC ERROR: DMPlexGetCellRefiner_Internal() line 6777 in /Users/eikemueller/PostDocBath/EllipticSolvers/petsc/src/dm/impls/plex/plexrefine.c
Unknown dimension 1 for cell refiner
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 62.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
I want to work in dimension=2+1, but having dimension=1+1 would be useful for testing.
There are two parts of this that won't work:1. DMPlex doesn't know how to refine intervals: that's the error above.I think this is easy to add, so can have a go.2. I haven't added the necessary numbering magic and so forth for the generated interval mesh hierarchies.If 1. is done, I think this should not be too difficult to add, but might take a little while.Lawrence_______________________________________________firedrake mailing listfiredrake@imperial.ac.ukhttps://mailman.ic.ac.uk/mailman/listinfo/firedrake