[deal.II] obtaining an ellipsoid grid and manifold from GridGenerator::hyper_sphere() and a subclass of ChartManifold

I am trying to obtain the mesh for a codimension-1 ellipsoid and attach an ellipsoidal manifold to it in order to refine the mesh. My strategy is failin

Re: [deal.II] Re: global dof renumbering with parallel triangulation and trilinos direct solver (Bug ?)

On 08/08/2016 02:18 PM, Denis Davydov wrote: However, when I changed the partition into something like { 0, 2, 4, ... } , { 1, 3, 5, ... } it fails, which makes me believe that non-contiguous partitions are not (completely) supported by dealii or trilinos should be fine for Tr

[deal.II] Re: global dof renumbering with parallel triangulation and trilinos direct solver (Bug ?)

On Monday, August 8, 2016 at 9:38:20 PM UTC+2, Daniel Jodlbauer wrote: > > - Debug is enabled (at least for dealii, I will have to rebuild trilinos > with debug later) > i mean deal.II. No need to recompile trilinos. > - I am not sure if I got you correctly, but If I use a regular > Triangul

[deal.II] Re: Embedding a Python script

Hi Ode, First, Spack is a tool to compile and install packages, which I don't think you are using as you said that you play around with Luca's build. So the link to this wiki should not be relevant here. Second, you CMake input is wrong. you need to drop DEAL_II_INVOKE_AUTOPILOT() see the ver

[deal.II] Re: global dof renumbering with parallel triangulation and trilinos direct solver (Bug ?)

Hi, - Debug is enabled (at least for dealii, I will have to rebuild trilinos with debug later) - I am not sure if I got you correctly, but If I use a regular Triangulation, then every rank owns all dofs and finally the initialization of the distributed vectors fails (as expected) What I additi

[deal.II] Re: global dof renumbering with parallel triangulation and trilinos direct solver (Bug ?)

Hi Daniel, First, did you try running in debug? If not,please do so as more checks will be triggered and it would be easier to debug. Second, can you please try it with usual Triangulation, not parallel::shared? If the problem persists then at least we know that it's not related to the triangul

[deal.II] global dof renumbering with parallel triangulation and trilinos direct solver (Bug ?)

Hi all ! I am trying to renumber the degrees of freedom globally using code like this: vector new_number(dof_handler.n_dofs()); for (unsigned int i = 0; i < dof_handler.n_dofs(); i++) new_number[i] = dof_handler.n_dofs() - i - 1; // simple example vector local_new_number; for

[deal.II] Embedding a Python script

Dear group, I apologize if my question below has already been addressed in the one of deal.ii’s documentation resources, but I didn't understand how to answer it using the pieces of documentation that I have read such as “How to use CMake to configure your projects with deal.II” https://

Re: [deal.II] Stresses at support point locations

Indeed, one commonly for stress recovery used approach is nodal averaging. See for example Zienkiewicz and Zhu "A Simple Estimator and Adaptive Procedure for Practical Engineering Analysis" http://onlinelibrary.wiley.com/doi/10.1002/nme.1620240206/abstract You can obtain the best fit in L2 norm