[CIG-LONG] Install Gale 2.0.0 with Mumps on Ubuntu

[Walter Landry]

Walter Landry <walter at geodynamics.org> wrote:
> nicolas riel <a.nicolas.riel at gmail.com> wrote:
>> The screenshot represents the thermal conditions at Step 1.
>> 
>> I'm using the extension.json I sent to you in a previous mail.
>> 
>> When I run Gale without mpirun it works very well, However when I try it in
>> parallel there is this error.
>> 
>> So maybe I should use mpich or another solver ?
>> 
>> Maybe my installation is wrong but, tried to install it on a  brand new
>> linux (ubuntu 11.10) to avoid conflicts...
> 
> That input file crashes for me if I run with more than one core, so it
> is not your installation.  I will look into it and get back to you.

I have figured out the problem.  The mesh for velocity and temperature
are different.  The velocity mesh is quadratic, so it can have curved
boundaries.  The temperature mesh is linear, so the nodes are
connected by straight lines.  This means that when the mesh is
distorted, points that are inside the velocity mesh may not be inside
the temperature mesh.

A temporary fix for this is to make temperature use the same mesh as
the velocity.  If you can compile code from the repository [1], I have
committed a fix.  The only change you might want to make is to remove
"T-mesh" from your EulerDeform struct, since it is no longer needed.

You can also fix it by adding the following items to the end of the
"components" struct

        ,
        "temperature":
        {
            "Type": "MeshVariable",
            "Rank": "Scalar",
            "DataType": "Double",
            "mesh": "v-mesh"
        },
        "temperatureBCs":
        {
            "Type": "CompositeVC",
            "Data": "v-mesh"
        },
        "temperatureICs":
        {
            "Type": "CompositeVC",
            "Data": "v-mesh"
        },
        "temperatureDofLayout":
        {
            "Type": "DofLayout",
            "mesh": "v-mesh",
            "BaseVariables": [
                "temperature"
            ]
        },
        "TemperatureField":
        {
            
            "Type": "FeVariable",
            "FEMesh": "v-mesh",
            "DofLayout": "temperatureDofLayout",
            "BC": "temperatureBCs",
            "IC": "temperatureICs",
            "LinkedDofInfo": "temperatureLinkedDofs"
        }

This overrides the defaults for these components.  In the
visualization output, you will still see "T-mesh", but it will not
have any data.  The temperature will instead be on "v-mesh".

I will probably make a new release in a few weeks that includes this
fix.

Cheers,
Walter Landry

[1] Note that the repository now requires Petsc 3.2.