Re: [deal.II] Re: Proper use of CompositionManifold

2016-09-21 Thread luca.heltai
For future reference, if you are in doubt on how to use a single class/function, you can skim through the directory tests, i.e., tests/manifold/composition_manifold_* would have shown to you an example which is identical to what you are trying to achieve, and in general it contains very simple

[deal.II] Re: Proper use of CompositionManifold

2016-09-21 Thread Praveen C
I was making a mistake. Reading the code of CompositionManifold clarified it. This is correct static const PolarManifold<2,3> sphere; static const FunctionManifold<3,3,3> rotate("x;-z;y","x;z;-y"); static const CompositionManifold<2,3,3,3,2,3> rotated_sphere(sphere, rotate);

[deal.II] Re: Proper use of CompositionManifold

2016-09-21 Thread Praveen C
Hello I think I understand this better after reading the code. Here is my attempt static const PolarManifold<2,3> sphere; static const FunctionManifold<3,3,3> rotate("x,-z,y","x,z,-y"); static CompositionManifold<2,3,2,3,3,3> surface(sphere, rotate); I use dim=2 The dimensi

Re: [deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Praveen C
Having curved boundaries alone is not sufficient sometimes. One would like to use curved elements in the interior also. It would be very useful if one could read high order meshes, e.g., from gmsh. The extra nodes can be used to create an euler vector that can be used in a MappingFEField. There ha

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-21 Thread Daniel Jodlbauer
Sounds like we ran into the same problem https://github.com/dealii/dealii/pull/2536, although is also experienced this issue in serial. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You receiv

[deal.II] Trilinos SparseMatrix mmult bug

2016-09-21 Thread Ben Shields
Hello everyone, So firstly, I'm using version 8.3.0. If this is a real bug but has been addressed in a newer version please let me know, but I have not found any info about it in the change logs up to version 8.4.2 or the user group. I'm using a block matrix to solve a steady, variable density

Re: [deal.II] compiling error "No rule to make target"

2016-09-21 Thread Andrew Lambe
That seems to have been the issue. I had to install libblas-dev, liblapack-dev, and libatlas-base-dev but the program installs normally now. Thanks, Andrew On Wednesday, September 21, 2016 at 11:40:59 AM UTC-4, Matthias Maier wrote: > > > On Wed, Sep 21, 2016, at 10:08 CDT, Andrew Lambe > wrot

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Daniel Arndt
Chenchen, We only support to read in meshes that are not refined. Translated to your situation this means that you are only able to use Q1 meshes as initial grid. Of course, you can then use triangulation.set_manifold to define the correct behavior when you refine your mesh within dealii along

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Bruno Turcksin
Chenchen, On Wednesday, September 21, 2016 at 11:56:31 AM UTC-4, Daniel Arndt wrote: > > We only support to read in meshes that are not refined. Translated to your > situation this means that you are only able to use Q1 meshes as initial > grid. > Of course, you can then use triangulation.set_ma

Re: [deal.II] compiling error "No rule to make target"

2016-09-21 Thread Matthias Maier
On Wed, Sep 21, 2016, at 10:08 CDT, Andrew Lambe wrote: > This seems odd since there is a file liblapack.so in the /usr/lib/ > directory. liblapack.so is typically a symlink to a versioned shared library liblapack.so.X. It might simply be a stale symlink (pointing to a nonexistent file). Hav

[deal.II] Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Chenchen Liu
Hi all, I have to use quadratic shape function (Q2 elements) to do my project (2D). Previously the code works well with Q1 linear elements, when I want to input the gird with Q2 elements generated by Gmsh, it gives the following error: An error occurred in line <1418> of file <../source/gr

[deal.II] compiling error "No rule to make target"

2016-09-21 Thread Andrew Lambe
This is probably a really simple question to answer, but I haven't found the fix yet. I recently upgraded my Ubuntu distribution to 16.04 and I'm trying to recompile deal.II (version 8.3) for the new OS. Note that I'm using version 3.4 of CMake which comes packaged with Ubuntu 16.04. Everything

[deal.II] Proper use of CompositionManifold

2016-09-21 Thread Praveen C
Dear all I am trying to compose two manifolds. I want to apply a 90 deg rotation about x-axis to PolarManifold to obtain description of surface of sphere and thereby avoid the poles. So I try something like this static PolarManifold<2,3> sphere; static FunctionManifold<3,3,3> rota