shells/fem/selement.h

Go to the documentation of this file.

/* -*- C -*- */
/*
 * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 *
 *                                   Fehmi Cirak
 *                        California Institute of Technology
 *                           (C) 2001 All Rights Reserved
 *
 * <LicenseText>
 *
 * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*/
#ifndef SELEMENT_H
#define SELEMENT_H

#if defined(HAVE_CONFIG_H)
/* autoconf */
#include <config.h>
#elif defined(BLD_PROCEDURE)
/* build procedure */
#include <portinfo.h>
#endif

#include <assert.h>
#include <stdlib.h>

#include "svertex.h"

#ifdef __cplusplus
namespace shells {
#define DECLARE_EXTERN extern "C"
#else
#define DECLARE_EXTERN extern    
#endif
    
    typedef struct SElementS {
        struct SElementS        *adt[3];     /* three neigbours of the element */
        SVertex                 *vtx[3];     /* three vertices of the element*/
        
        struct InterpolationS   *interpol;   /* shape functions */ 
        
        struct InternalStorageS *matStorage; /* storage for material */ 
    } SElement; 
    
    
    DECLARE_EXTERN SElement *createSElement();
    DECLARE_EXTERN SElement *createSElement2(SVertex *vtx0, 
                                             SVertex *vtx1, 
                                             SVertex *vtx2);
    DECLARE_EXTERN SElement *duplicateSElement(SElement *);
    DECLARE_EXTERN void deleteSElement(SElement *);
    
    DECLARE_EXTERN int selementsOneRingSVertex(SElement *element, 
                                               int inode, 
                                               SElement* oneRing[MAXVAL]);
    
    DECLARE_EXTERN int sverticesOneRingSVertex(SElement *element,
                                               int inode,
                                               SVertex* oneRing[MAXVAL]);

    DECLARE_EXTERN void selementPtrToHistory(SElement *element,
                                             double **data, 
                                             int *size);

    DECLARE_EXTERN void selementResetHistory(SElement  *element,
                                             double *data);
    
    static int localSEdgeNumber(SElement *t, SVertex *v1, SVertex *v2) 
        {
            int result;
            assert( (v1 != NULL) && (v2 != NULL) ); 
            assert( v1 != v2 );
        
            result = ( v1 == t->vtx[0] ? ( v2 == t->vtx[1] ? 0 : 2 )
                       : ( v1 == t->vtx[1] ? ( v2 == t->vtx[2] ? 1 : 0 ) :
                           ( v2 == t->vtx[0] ? 2 : 1 ) ) );
            return result;
        }

    static int localSEdgeNumber2(SElement *t, SElement *neighbor) 
        {
            int j;
            assert(t != NULL);
            assert(neighbor != NULL);
            j = (t->adt[0] == neighbor) ? 0 : ((t->adt[1] == neighbor) ? 1 : 
                                               ((t->adt[2] == neighbor) ? 2 : 3));
            assert(j!=3);
            
            return j;
        } 
        
    static int localSVtxNumber(SVertex * vtx, SElement * t) 
        {
            int j;
            assert(t != NULL);
            assert(vtx != NULL);
            j = (t->vtx[0] == vtx) ? 0 : ((t->vtx[1] == vtx) ? 1 : 
                                          ((t->vtx[2] == vtx) ? 2 : 3));
            
            return j;
        }
    
    static void edgeTangentInParamSpace(double t[3], unsigned iedg)
        {
            switch (iedg) {
            case 0:
                t[0] = 1.0;   
                t[1] = 0.0;                   
                t[2] = 0.0;
                break;
            case 1:
                t[0] =-1.0;                 
                t[1] = 1.0;            
                t[2] = 0.0;
                break;
            case 2:
                t[0] = 0.0;                 
                t[1] =-1.0;                    
                t[2] = 0.0;
                break;  
            default:
                assert(0);
                break;
            }
        }
    
    static void resetSElementSVertices(SElement *element, SVertex *vtx0, 
                                       SVertex *vtx1, SVertex *vtx2) 
        {
            element->vtx[0] = vtx0;
            element->vtx[1] = vtx1;
            element->vtx[2] = vtx2;
        }
    
    static void resetSElementNeighbors(SElement *element, SElement *adt0, 
                                       SElement *adt1, SElement *adt2) 
        {
            element->adt[0] = adt0;
            element->adt[1] = adt1;
            element->adt[2] = adt2;
        }

    
#ifdef __cplusplus
}
#endif
#undef DECLARE_EXTERN

#endif
/* End of file */

---