~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to src/kernel/grid/Tetrahedron.cpp

Committer: logg
Date: 2003-02-06 14:40:25 UTC
Revision ID: devnull@localhost-20030206144025-bmqurqq95f073l40

Tailorized "2003-02-06 08:40:24 by logg"
Fixes...

files removed:
src/kernel/fem/FunctionSpace.h

files modified:
ChangeLog

README

configure

configure.in

doc/policy

src/demo/Makefile.am

src/demo/Makefile.in

src/demo/fem/main.cpp

src/demo/solvers/Makefile.am

src/demo/solvers/navierstokes/Makefile.in

src/demo/solvers/poisson/main.cpp

src/demo/solvers/template/main.cpp

src/kernel/Makefile.am

src/kernel/Makefile.in

src/kernel/common/Makefile.am

src/kernel/common/Makefile.in

src/kernel/common/ParameterList.cpp

src/kernel/common/dolfin/List.h

src/kernel/common/dolfin/Makefile.am

src/kernel/common/dolfin/Makefile.in

src/kernel/common/dolfin/ShortList.h

src/kernel/common/dolfin/dolfin_common.h

src/kernel/common/utils.cpp

src/kernel/fem/ElementFunction.cpp

src/kernel/fem/Equation.cpp

src/kernel/fem/EquationSystem.cpp

src/kernel/fem/FiniteElement.cpp

src/kernel/fem/Function.cpp

src/kernel/fem/FunctionList.cpp

src/kernel/fem/FunctionSpace.cpp

src/kernel/fem/Galerkin.cpp

src/kernel/fem/Makefile.am

src/kernel/fem/Makefile.in

src/kernel/fem/Product.cpp

src/kernel/fem/ShapeFunction.cpp

src/kernel/fem/dolfin/ElementFunction.h

src/kernel/fem/dolfin/Function.h

src/kernel/fem/dolfin/FunctionList.h

src/kernel/fem/dolfin/FunctionSpace.h

src/kernel/fem/dolfin/Makefile.am

src/kernel/fem/dolfin/Makefile.in

src/kernel/fem/dolfin/Product.h

src/kernel/fem/dolfin/ShapeFunction.h

src/kernel/fem/dolfin/TriLinSpace.h

src/kernel/fem/dolfin/dolfin_fem.h

src/kernel/grid/Cell.cpp

src/kernel/grid/CellIterator.cpp

src/kernel/grid/Grid.cpp

src/kernel/grid/GridData.cpp

src/kernel/grid/GridData.h

src/kernel/grid/InitGrid.cpp

src/kernel/grid/Node.cpp

src/kernel/grid/NodeIterator.cpp

src/kernel/grid/Point.cpp

src/kernel/grid/Tetrahedron.cpp

src/kernel/grid/Triangle.cpp

src/kernel/grid/dolfin/Cell.h

src/kernel/grid/dolfin/CellIterator.h

src/kernel/grid/dolfin/GenericCell.h

src/kernel/grid/dolfin/Grid.h

src/kernel/grid/dolfin/Node.h

src/kernel/grid/dolfin/NodeIterator.h

src/kernel/grid/dolfin/Point.h

src/kernel/grid/dolfin/Tetrahedron.h

src/kernel/grid/dolfin/Triangle.h

src/kernel/io/DataInfo.cpp

src/kernel/io/Makefile.am

src/kernel/io/Makefile.in

src/kernel/io/XMLFile.cpp

src/kernel/io/XMLFile.h

src/kernel/io/XMLGrid.cpp

src/kernel/io/XMLGrid.h

src/kernel/io/XMLObject.h

src/kernel/io/XMLSparseMatrix.cpp

src/kernel/io/XMLVector.cpp

src/kernel/io/dolfin/Makefile.am

src/kernel/io/dolfin/Makefile.in

src/kernel/io/dolfin/dolfin_io.h

src/kernel/la/DenseMatrix.cpp

src/kernel/la/DirectSolver.cpp

src/kernel/la/KrylovSolver.cpp

src/kernel/la/Makefile.am

src/kernel/la/Makefile.in

src/kernel/la/SISolver.cpp

src/kernel/la/SparseMatrix.cpp

src/kernel/la/Vector.cpp

src/kernel/la/dolfin/Makefile.am

src/kernel/la/dolfin/Makefile.in

src/kernel/la/dolfin/SparseMatrix.h

src/kernel/la/dolfin/Vector.h

src/kernel/la/dolfin/dolfin_la.h

src/kernel/main/Problem.cpp

src/kernel/main/dolfin.h

src/kernel/main/dolfin/Problem.h

src/kernel/main/dolfin/constants.h

src/modules/Makefile.am

src/modules/Makefile.in

src/modules/poisson/Poisson.h

src/modules/poisson/PoissonSettings.h

src/modules/poisson/PoissonSolver.cpp

src/modules/poisson/PoissonSolver.h

src/modules/scanmodules.sh

src/modules/template/Template.h

src/modules/template/TemplateSettings.h

src/modules/template/TemplateSolver.cpp

src/modules/template/TemplateSolver.h

Show diffs side-by-side

added added

removed removed

src/kernel/grid/Tetrahedron.cpp

// Licensed under the GNU GPL Version 2.

#include <dolfin/Node.h>

#include <dolfin/Point.h>

#include <dolfin/Cell.h>

#include <dolfin/Tetrahedron.h>

using namespace dolfin;

//-----------------------------------------------------------------------------

int Tetrahedron::noNodes()

int Tetrahedron::noNodes() const

{

return 4;

}

//-----------------------------------------------------------------------------

int Tetrahedron::noEdges()

int Tetrahedron::noEdges() const

{

return 6;

}

//-----------------------------------------------------------------------------

int Tetrahedron::noFaces()

int Tetrahedron::noFaces() const

{

return 4;

}

//-----------------------------------------------------------------------------

int Tetrahedron::noBoundaries()

int Tetrahedron::noBound() const

{

return noFaces();

}

//-----------------------------------------------------------------------------

Cell::Type Tetrahedron::type()

Cell::Type Tetrahedron::type() const

{

return Cell::TETRAHEDRON;

}

//-----------------------------------------------------------------------------

real Tetrahedron::ComputeVolume(Grid *grid)

{

// Get the coordinates

Point *A = nodes[0]->GetCoord();

Point *B = nodes[1]->GetCoord();

Point *C = nodes[2]->GetCoord();

Point *D = nodes[3]->GetCoord();

// Make sure we get full precision

real x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4;

x1 = real(A->x); y1 = real(A->y); z1 = real(A->z);

x2 = real(B->x); y2 = real(B->y); z2 = real(B->z);

x3 = real(C->x); y3 = real(C->y); z3 = real(C->z);

x4 = real(D->x); y4 = real(D->y); z4 = real(D->z);

// Formula for volume from http://mathworld.wolfram.com

real v = ( x1 * ( y2*z3 + y4*z2 + y3*z4 - y3*z2 - y2*z4 - y4*z3 ) -

x2 * ( y1*z3 + y4*z1 + y3*z4 - y3*z1 - y1*z4 - y4*z3 ) +

x3 * ( y1*z2 + y4*z1 + y2*z4 - y2*z1 - y1*z4 - y4*z2 ) -

x4 * ( y1*z2 + y2*z3 + y3*z1 - y2*z1 - y3*z2 - y1*z3 ) );

//display->Message(0,"A = (%f,%f,%f)\nB = (%f,%f,%f)\nC = (%f,%f,%f)\nD = (%f,%f,%f)\n",

// x1,y1,z1,

// x2,y2,z2,

// x3,y3,z3,

// x4,y4,z4);

v = ( v >= 0.0 ? v : -v );

return ( v );

}

//-----------------------------------------------------------------------------

real Tetrahedron::ComputeCircumRadius(Grid *grid)

{

// Compute volume

real volume = ComputeVolume(grid);

// Compute radius

real radius = ComputeCircumRadius(grid,volume);

//display->Message(0,"v = %f R = %f",volume,radius);

// Return radius

return ( radius );

}

//-----------------------------------------------------------------------------

real Tetrahedron::ComputeCircumRadius(Grid *grid, real volume)

{

// Get the coordinates

Point *A = nodes[0]->GetCoord();

Point *B = nodes[1]->GetCoord();

Point *C = nodes[2]->GetCoord();

Point *D = nodes[3]->GetCoord();

// Compute side lengths

real a = B->Distance(*C);

real b = A->Distance(*C);

real c = A->Distance(*B);

real aa = A->Distance(*D);

real bb = B->Distance(*D);

real cc = C->Distance(*D);

// Compute "area" of triangle with strange side lengths

100

real l1 = a*aa;

101

real l2 = b*bb;

102

real l3 = c*cc;

103

real s = 0.5*(l1+l2+l3);

104

real area = sqrt(s*(s-l1)*(s-l2)*(s-l3));

105

106

// Formula for radius from http://mathworld.wolfram.com

107

real R = area/(6.0*volume);

108

109

return ( R );

110

}

111

112

//-----------------------------------------------------------------------------

113

bool Tetrahedron::neighbor(ShortList<Node *> &cn, Cell &cell)

bool Tetrahedron::neighbor(ShortList<Node *> &cn, Cell &cell) const

114

{

115

// Two tetrahedrons are neighbors if they have a common face or if they are

116

// the same tetrahedron, i.e. if they have 3 or 4 common nodes.

Older »