Changeset 5908


Ignore:
Timestamp:
09/20/10 13:50:57 (14 years ago)
Author:
Mathieu Morlighem
Message:

NewElementMatrix and NewElementVector are now in shared, common to all elements and loads

Location:
issm/trunk/src/c
Files:
5 added
10 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/Makefile.am

    r5906 r5908  
    269269                                        ./shared/Elements/Paterson.cpp\
    270270                                        ./shared/Elements/GetVerticesCoordinates.cpp\
     271                                        ./shared/Elements/GetLocalDofList.cpp\
     272                                        ./shared/Elements/GetGlobalDofList.cpp\
     273                                        ./shared/Elements/GetNumberOfDofs.cpp\
     274                                        ./shared/Elements/NewElementMatrix.cpp\
     275                                        ./shared/Elements/NewElementVector.cpp\
    271276                                        ./shared/String/DescriptorIndex.cpp\
    272277                                        ./shared/String/sharedstring.h\
    … …  
    824829                                        ./shared/Elements/Paterson.cpp\
    825830                                        ./shared/Elements/GetVerticesCoordinates.cpp\
     831                                        ./shared/Elements/GetLocalDofList.cpp\
     832                                        ./shared/Elements/GetGlobalDofList.cpp\
     833                                        ./shared/Elements/GetNumberOfDofs.cpp\
     834                                        ./shared/Elements/NewElementMatrix.cpp\
     835                                        ./shared/Elements/NewElementVector.cpp\
    826836                                        ./shared/String/DescriptorIndex.cpp\
    827837                                        ./shared/String/sharedstring.h\

  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r5887 r5908  
    20482048        /*Initialize Element matrix and return if necessary*/
    20492049        if(IsOnWater()) return NULL;
    2050         ElementMatrix* Ke1=pentabase->NewElementMatrix(MacAyealApproximationEnum);
    2051         ElementMatrix* Ke2=this->NewElementMatrix(PattynApproximationEnum);
     2050        ElementMatrix* Ke1=NewElementMatrix(pentabase->nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
     2051        ElementMatrix* Ke2=NewElementMatrix(this->nodes     ,NUMVERTICES,this->parameters,PattynApproximationEnum);
    20522052        ElementMatrix* Ke=new ElementMatrix(Ke1,Ke2);
    20532053        delete Ke1; delete Ke2;
    … …  
    21962196        /*Initialize Element matrix and return if necessary*/
    21972197        if(IsOnWater()) return NULL;
    2198         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2198        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    21992199
    22002200        /*Spawn 3 beam elements: */
    … …  
    23192319        /*Initialize Element matrix and return if necessary*/
    23202320        if(IsOnWater()) return NULL;
    2321         ElementMatrix* Ke=tria->NewElementMatrix(MacAyealApproximationEnum);
     2321        ElementMatrix* Ke=NewElementMatrix(tria->nodes,NUMVERTICES2D,this->parameters,MacAyealApproximationEnum);
    23222322
    23232323        /*Retrieve all inputs and parameters*/
    … …  
    24472447        /*Initialize Element matrix and return if necessary*/
    24482448        if(IsOnWater()) return NULL;
    2449         ElementMatrix* Ke=this->NewElementMatrix(PattynApproximationEnum);
     2449        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,PattynApproximationEnum);
    24502450
    24512451        /*Retrieve all inputs and parameters*/
    … …  
    25672567        inputs->GetParameterValue(&approximation,ApproximationEnum);
    25682568        if(IsOnWater() || approximation!=StokesApproximationEnum) return NULL;
    2569         ElementMatrix* Ke=this->NewElementMatrix(StokesApproximationEnum);
     2569        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,StokesApproximationEnum);
    25702570
    25712571        /*Retrieve all inputs and parameters*/
    … …  
    26482648        inputs->GetParameterValue(&approximation,ApproximationEnum);
    26492649        if(IsOnWater() || IsOnShelf() || !IsOnBed() || approximation!=StokesApproximationEnum) return NULL;
    2650         ElementMatrix* Ke=this->NewElementMatrix(StokesApproximationEnum);
     2650        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,StokesApproximationEnum);
    26512651
    26522652        /*Retrieve all inputs and parameters*/
    … …  
    27392739        /*Initialize Element matrix and return if necessary*/
    27402740        if(IsOnWater()) return NULL;
    2741         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2741        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    27422742
    27432743        /*Retrieve all inputs and parameters*/
    … …  
    28852885        /*Initialize Element matrix and return if necessary*/
    28862886        if(IsOnWater()) return NULL;
    2887         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2887        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    28882888
    28892889        /*Retrieve all inputs and parameters*/
    … …  
    55155515}
    55165516/*}}}*/
    5517 /*FUNCTION Penta::NewElementMatrix{{{1*/
    5518 ElementMatrix* Penta::NewElementMatrix(int approximation){
    5519 
    5520         /*parameters: */
    5521         bool kff=false;
    5522 
    5523         /*output: */
    5524         ElementMatrix* Ke=NULL;
    5525         int gsize;
    5526         int fsize;
    5527         int ssize;
    5528         int* gglobaldoflist=NULL;
    5529         int* flocaldoflist=NULL;
    5530         int* fglobaldoflist=NULL;
    5531         int* slocaldoflist=NULL;
    5532         int* sglobaldoflist=NULL;
    5533         bool square=true;
    5534 
    5535         /*retrieve some parameters: */
    5536         this->parameters->FindParam(&kff,KffEnum);
    5537 
    5538         /*get number of dofs in sets g,f and s: */
    5539         gsize=this->GetNumberOfDofs(approximation,GsetEnum);
    5540         if(kff){
    5541                 fsize=this->GetNumberOfDofs(approximation,FsetEnum);
    5542                 ssize=this->GetNumberOfDofs(approximation,SsetEnum);
    5543         }
    5544 
    5545         /*get dof lists for f and s set: */
    5546         gglobaldoflist=this->GetGlobalDofList(approximation,GsetEnum);
    5547         if(kff){
    5548                 flocaldoflist=this->GetLocalDofList(approximation,FsetEnum);
    5549                 fglobaldoflist=this->GetGlobalDofList(approximation,FsetEnum);
    5550                 slocaldoflist=this->GetLocalDofList(approximation,SsetEnum);
    5551                 sglobaldoflist=this->GetGlobalDofList(approximation,SsetEnum);
    5552         }
    5553 
    5554         /*Use square constructor for ElementMatrix: */
    5555         if(!kff) Ke=new ElementMatrix(square,gglobaldoflist,gsize);
    5556         else     Ke=new ElementMatrix(square,gglobaldoflist,gsize,flocaldoflist,fglobaldoflist,fsize,slocaldoflist,sglobaldoflist,ssize);
    5557 
    5558         /*Free ressources and return:*/
    5559         xfree((void**)&gglobaldoflist);
    5560         xfree((void**)&flocaldoflist);
    5561         xfree((void**)&fglobaldoflist);
    5562         xfree((void**)&slocaldoflist);
    5563         xfree((void**)&sglobaldoflist);
    5564 
    5565         return Ke;
    5566 }
    5567 /*}}}*/
    5568 /*FUNCTION Penta::NewElementVector{{{1*/
    5569 ElementVector* Penta::NewElementVector(int approximation){
    5570 
    5571         /*parameters: */
    5572         bool kff=false;
    5573 
    5574         /*output: */
    5575         ElementVector* pe=NULL;
    5576         int gsize;
    5577         int fsize;
    5578         int* gglobaldoflist=NULL;
    5579         int* flocaldoflist=NULL;
    5580         int* fglobaldoflist=NULL;
    5581 
    5582         /*retrieve some parameters: */
    5583         this->parameters->FindParam(&kff,KffEnum);
    5584 
    5585         /*get number of dofs in sets g,f and s: */
    5586         gsize=this->GetNumberOfDofs(approximation,GsetEnum);
    5587         if(kff)fsize=this->GetNumberOfDofs(approximation,FsetEnum);
    5588 
    5589         /*get dof lists for f and s set: */
    5590         if(!kff){
    5591                 gglobaldoflist=this->GetGlobalDofList(approximation,GsetEnum);
    5592         }
    5593         else{
    5594                 flocaldoflist=this->GetLocalDofList(approximation,FsetEnum);
    5595                 fglobaldoflist=this->GetGlobalDofList(approximation,FsetEnum);
    5596         }
    5597 
    5598         /*constructor for ElementVector: */
    5599         if(!kff)pe=new ElementVector(gsize,gglobaldoflist);
    5600         else    pe=new ElementVector(gsize,flocaldoflist,fglobaldoflist,fsize);
    5601 
    5602         /*Free ressources and return:*/
    5603         xfree((void**)&gglobaldoflist);
    5604         xfree((void**)&flocaldoflist);
    5605         xfree((void**)&fglobaldoflist);
    5606 
    5607         return pe;
    5608 }
    5609 /*}}}*/
    56105517/*FUNCTION Penta::ReduceMatrixStokes {{{1*/
    56115518void Penta::ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp){

  • issm/trunk/src/c/objects/Elements/Penta.h

    r5879 r5908  
    204204                bool    IsOnShelf(void);
    205205                bool    IsOnWater(void);
    206                 ElementMatrix* NewElementMatrix(int approximation);
    207                 ElementVector* NewElementVector(int approximation);
    208206                void      ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
    209207                void      ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);

  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r5881 r5908  
    24292429        /*Initialize Element matrix and return if necessary*/
    24302430        if(IsOnWater()) return NULL;
    2431         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2431        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    24322432
    24332433        /*Retrieve all Inputs and parameters: */
    … …  
    25342534        /*Initialize Element matrix and return if necessary*/
    25352535        if(IsOnWater()) return NULL;
    2536         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2536        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    25372537
    25382538        /*Retrieve all inputs and parameters*/
    … …  
    26012601        /*Initialize Element matrix and return if necessary*/
    26022602        if(IsOnWater()) return NULL;
    2603         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     2603        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    26042604
    26052605        /*Retrieve all inputs and parameters*/
    … …  
    27332733        /*Initialize Element matrix and return if necessary*/
    27342734        if(IsOnWater()) return NULL;
    2735         ElementMatrix* Ke=this->NewElementMatrix(MacAyealApproximationEnum);
     2735        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
    27362736
    27372737        /*Retrieve all inputs and parameters*/
    … …  
    28012801        /*Initialize Element matrix and return if necessary*/
    28022802        if(IsOnWater() || IsOnShelf()) return NULL;
    2803         ElementMatrix* Ke=this->NewElementMatrix(MacAyealApproximationEnum);
     2803        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
    28042804
    28052805        /*Retrieve all inputs and parameters*/
    … …  
    28722872        /*Initialize Element matrix and return if necessary*/
    28732873        if(IsOnWater() || IsOnShelf()) return NULL;
    2874         ElementMatrix* Ke1=this->NewElementMatrix(MacAyealApproximationEnum);
    2875         ElementMatrix* Ke2=this->NewElementMatrix(PattynApproximationEnum);
     2874        ElementMatrix* Ke1=NewElementMatrix(nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
     2875        ElementMatrix* Ke2=NewElementMatrix(nodes,NUMVERTICES,this->parameters,PattynApproximationEnum);
    28762876        ElementMatrix* Ke=new ElementMatrix(Ke1,Ke2);
    28772877        delete Ke1; delete Ke2;
    … …  
    29612961        /*Initialize Element matrix and return if necessary*/
    29622962        if(IsOnWater() || IsOnShelf()) return NULL;
    2963         ElementMatrix* Ke=this->NewElementMatrix(PattynApproximationEnum);
     2963        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,PattynApproximationEnum);
    29642964
    29652965        /*Retrieve all inputs and parameters*/
    … …  
    30213021        /*Initialize Element matrix and return if necessary*/
    30223022        if(IsOnWater()) return NULL;
    3023         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3023        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    30243024
    30253025        /*Create Element matrix*/
    … …  
    30523052        /*Initialize Element matrix and return if necessary*/
    30533053        if(IsOnWater()) return NULL;
    3054         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3054        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    30553055
    30563056        /*Retrieve all inputs and parameters*/
    … …  
    31023102        /*Initialize Element matrix and return if necessary*/
    31033103        if(IsOnWater()) return NULL;
    3104         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3104        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    31053105
    31063106        /*Retrieve all inputs and parameters*/
    … …  
    31733173        /*Initialize Element matrix and return if necessary*/
    31743174        if(IsOnWater()) return NULL;
    3175         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3175        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    31763176
    31773177        /*Retrieve all inputs and parameters*/
    … …  
    32933293        /*Initialize Element matrix and return if necessary*/
    32943294        if(IsOnWater()) return NULL;
    3295         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3295        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    32963296
    32973297        /*Retrieve all inputs and parameters*/
    … …  
    33683368        /*Initialize Element matrix and return if necessary*/
    33693369        if(IsOnWater()) return NULL;
    3370         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3370        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    33713371
    33723372        GetVerticesCoordinates(&xyz_list[0][0],nodes,NUMVERTICES);
    … …  
    34193419        /*Initialize Element matrix and return if necessary*/
    34203420        if(IsOnWater() || !IsOnShelf()) return NULL;
    3421         ElementMatrix* Ke=this->NewElementMatrix(NoneApproximationEnum);
     3421        ElementMatrix* Ke=NewElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
    34223422
    34233423        /*Retrieve all inputs and parameters*/
    … …  
    37443744
    37453745        /*Initialize element vector: */
    3746         pe=this->NewElementVector(MacAyealApproximationEnum);
     3746        pe=NewElementVector(nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
    37473747
    37483748        /*Add pe_g values to pe element stifness load: */
    … …  
    54535453}
    54545454/*}}}*/
    5455 /*FUNCTION Tria::NewElementMatrix{{{1*/
    5456 ElementMatrix* Tria::NewElementMatrix(int approximation){
    5457 
    5458         /*parameters: */
    5459         bool kff=false;
    5460 
    5461         /*output: */
    5462         ElementMatrix* Ke=NULL;
    5463         int gsize;
    5464         int fsize;
    5465         int ssize;
    5466         int* gglobaldoflist=NULL;
    5467         int* flocaldoflist=NULL;
    5468         int* fglobaldoflist=NULL;
    5469         int* slocaldoflist=NULL;
    5470         int* sglobaldoflist=NULL;
    5471         bool square=true;
    5472 
    5473         /*retrieve some parameters: */
    5474         this->parameters->FindParam(&kff,KffEnum);
    5475 
    5476         /*get number of dofs in sets g,f and s: */
    5477         gsize=this->GetNumberOfDofs(approximation,GsetEnum);
    5478         if(kff){
    5479                 fsize=this->GetNumberOfDofs(approximation,FsetEnum);
    5480                 ssize=this->GetNumberOfDofs(approximation,SsetEnum);
    5481         }
    5482 
    5483         /*get dof lists for f and s set: */
    5484         gglobaldoflist=this->GetGlobalDofList(approximation,GsetEnum);
    5485         if(kff){
    5486                 flocaldoflist=this->GetLocalDofList(approximation,FsetEnum);
    5487                 fglobaldoflist=this->GetGlobalDofList(approximation,FsetEnum);
    5488                 slocaldoflist=this->GetLocalDofList(approximation,SsetEnum);
    5489                 sglobaldoflist=this->GetGlobalDofList(approximation,SsetEnum);
    5490         }
    5491 
    5492         /*Use square constructor for ElementMatrix: */
    5493         if(!kff) Ke=new ElementMatrix(square,gglobaldoflist,gsize);
    5494         else     Ke=new ElementMatrix(square,gglobaldoflist,gsize,flocaldoflist,fglobaldoflist,fsize,slocaldoflist,sglobaldoflist,ssize);
    5495 
    5496         /*Free ressources and return:*/
    5497         xfree((void**)&gglobaldoflist);
    5498         xfree((void**)&flocaldoflist);
    5499         xfree((void**)&fglobaldoflist);
    5500         xfree((void**)&slocaldoflist);
    5501         xfree((void**)&sglobaldoflist);
    5502 
    5503         return Ke;
    5504 }
    5505 /*}}}*/
    5506 /*FUNCTION Tria::NewElementVector{{{1*/
    5507 ElementVector* Tria::NewElementVector(int approximation){
    5508 
    5509         /*parameters: */
    5510         bool kff=false;
    5511 
    5512         /*output: */
    5513         ElementVector* pe=NULL;
    5514         int gsize;
    5515         int fsize;
    5516         int* gglobaldoflist=NULL;
    5517         int* flocaldoflist=NULL;
    5518         int* fglobaldoflist=NULL;
    5519 
    5520         /*retrieve some parameters: */
    5521         this->parameters->FindParam(&kff,KffEnum);
    5522 
    5523         /*get number of dofs in sets g,f and s: */
    5524         gsize=this->GetNumberOfDofs(approximation,GsetEnum);
    5525         if(kff)fsize=this->GetNumberOfDofs(approximation,FsetEnum);
    5526 
    5527         /*get dof lists for f and s set: */
    5528         if(!kff){
    5529                 gglobaldoflist=this->GetGlobalDofList(approximation,GsetEnum);
    5530         }
    5531         else{
    5532                 flocaldoflist=this->GetLocalDofList(approximation,FsetEnum);
    5533                 fglobaldoflist=this->GetGlobalDofList(approximation,FsetEnum);
    5534         }
    5535 
    5536         /*constructor for ElementVector: */
    5537         if(!kff)pe=new ElementVector(gsize,gglobaldoflist);
    5538         else    pe=new ElementVector(gsize,flocaldoflist,fglobaldoflist,fsize);
    5539 
    5540         /*Free ressources and return:*/
    5541         xfree((void**)&gglobaldoflist);
    5542         xfree((void**)&flocaldoflist);
    5543         xfree((void**)&fglobaldoflist);
    5544        
    5545         return pe;
    5546 }
    5547 /*}}}*/
    55485455/*FUNCTION Tria::SetClone {{{1*/
    55495456void  Tria::SetClone(int* minranks){

  • issm/trunk/src/c/objects/Elements/Tria.h

    r5881 r5908  
    172172                void      InputUpdateFromSolutionOneDof(double* solution,int enum_type);
    173173                bool      IsInput(int name);
    174                 ElementMatrix* NewElementMatrix(int approximation);
    175                 ElementVector* NewElementVector(int approximation);
    176174                void      SetClone(int* minranks);
    177175                void      SurfaceNormal(double* surface_normal, double xyz_list[3][3]);

  • issm/trunk/src/c/objects/Loads/Icefront.cpp

    r5808 r5908  
    496496
    497497        /*Initialize element matrix: */
    498         pe=this->NewElementVector(MacAyealApproximationEnum);
     498        pe=NewElementVector(nodes,NUMVERTICESSEG,this->parameters,MacAyealApproximationEnum);
    499499
    500500        /*Add pe_g values to pe element stifness load: */
    … …  
    910910}
    911911/*}}}*/
    912 /*FUNCTION Icefront::NewElementVector{{{1*/
    913 ElementVector* Icefront::NewElementVector(int approximation){
    914 
    915         /*parameters: */
    916         bool kff=false;
    917 
    918         /*output: */
    919         ElementVector* pe=NULL;
    920         int gsize;
    921         int fsize;
    922         int* gexternaldoflist=NULL;
    923         int* finternaldoflist=NULL;
    924         int* fexternaldoflist=NULL;
    925 
    926         /*retrieve some parameters: */
    927         this->parameters->FindParam(&kff,KffEnum);
    928 
    929         /*get number of dofs in sets g and f: */
    930         gsize=this->GetNumberOfDofs(approximation,GsetEnum);
    931         if(kff)fsize=this->GetNumberOfDofs(approximation,FsetEnum);
    932 
    933         /*get dof lists for f and s set: */
    934         if(!kff){
    935                 gexternaldoflist=this->GetExternalDofList(approximation,GsetEnum);
    936         }
    937         else{
    938                 finternaldoflist=this->GetInternalDofList(approximation,FsetEnum);
    939                 fexternaldoflist=this->GetExternalDofList(approximation,FsetEnum);
    940         }
    941 
    942         /*Use symmetric constructor for ElementVector: */
    943         if(!kff)pe=new ElementVector(gsize,gexternaldoflist);
    944         else    pe=new ElementVector(gsize,finternaldoflist,fexternaldoflist,fsize);
    945 
    946         /*Free ressources and return:*/
    947         xfree((void**)&gexternaldoflist);
    948         xfree((void**)&finternaldoflist);
    949         xfree((void**)&fexternaldoflist);
    950        
    951         return pe;
    952 }
    953 /*}}}*/

  • issm/trunk/src/c/objects/Loads/Icefront.h

    r5808 r5908  
    8787                void GetSegmentNormal(double* normal,double xyz_list[2][3]);
    8888                void GetQuadNormal(double* normal,double xyz_list[4][3]);
    89                 ElementVector* NewElementVector(int approximation);
    9089                /*}}}*/
    9190};

  • issm/trunk/src/c/objects/objects.h

    r5772 r5908  
    7777#include "./Numerics/ElementVector.h"
    7878
    79 
    8079/*Params: */
    8180#include "./Params/BoolParam.h"

  • issm/trunk/src/c/shared/Elements/GetVerticesCoordinates.cpp

    r5906 r5908  
    33 */
    44
    5 #include "../../objects/objects.h"
     5#include "./elements.h"
    66
    7 int GetVerticesCoordinates( double* xyz,  Node** nodes, int numvertices){
     7int GetVerticesCoordinates(double* xyz,  Node** nodes, int numvertices){
    88
    99        /*In debugging mode, check that nodes is not a NULL pointer*/

  • issm/trunk/src/c/shared/Elements/elements.h

    r5906 r5908  
    88#include "../../objects/objects.h"
    99#include "../../Container/Container.h"
     10class ElementMatrix;
     11class ElementVector;
    1012
    11 double Paterson(double temperature);
    12 int    GetVerticesCoordinates(double* xyz,  Node** nodes, int numvertices);
     13double         Paterson(double temperature);
     14int            GetVerticesCoordinates(double* xyz,  Node** nodes, int numvertices);
     15int            GetNumberOfDofs( Node** nodes,int numnodes,int setenum,int approximation_enum);
     16int*           GetLocalDofList( Node** nodes,int numnodes,int setenum,int approximation_enum);
     17int*           GetGlobalDofList(Node** nodes,int numnodes,int setenum,int approximation_enum);
     18ElementMatrix* NewElementMatrix(Node** nodes,int numnodes,Parameters* parameters,int approximation=NoneApproximationEnum);
     19ElementVector* NewElementVector(Node** nodes,int numnodes,Parameters* parameters,int approximation=NoneApproximationEnum);
     20//ElementVector* NewElementVector(Nodes** nodes,int numnodes,Parameters* parameters,int approximation=NoneApproximationEnum);
    1321
    1422inline void printarray(double* array,int lines,int cols=1){
Note: See TracChangeset for help on using the changeset viewer.