Changeset 15037


Ignore:
Timestamp:
05/16/13 09:22:27 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: cosmetics (removed double blank lines)

Location:
issm/trunk-jpl/src/c
Files:
14 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/hydrology_core.cpp

    r15015 r15037  
    6363
    6464                        /*transfer water column thickness to old water column thickness: */
    65                        
     65
    6666                        InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,WatercolumnEnum,WaterColumnOldEnum);
    6767
    … …  
    8383
    8484                        if(VerboseSolution()) _pprintLine_("   computing water transfer");
    85                        
    8685
    8786                        if(VerboseSolution()) _pprintLine_("   computing water head");

  • issm/trunk-jpl/src/c/classes/Constraints/Constraints.h

    r15012 r15037  
    2828                /*}}}*/
    2929
    30 
    3130                /*numerics*/
    3231                int   NumberOfConstraints(void);

  • issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

    r15015 r15037  
    93569356}
    93579357/*}}}*/
    9358        
    93599358
    93609359#endif

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r15020 r15037  
    19141914                /*Get dof list: */
    19151915                GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    1916                
     1916
    19171917                /*Use the dof list to index into the vector: */
    19181918                for(int i=0;i<numdof;i++){
    … …  
    19221922                /*Add input to the element: */
    19231923                this->inputs->AddInput(new TriaP1Input(name,values));
    1924                
     1924
    19251925                /*Free ressources:*/
    19261926                xDelete<int>(doflist);
    … …  
    61736173        /* Start  looping on the number of gaussian points: */
    61746174        gauss=new GaussTria(2);
    6175        
     6175
    61766176        for(int ig=gauss->begin();ig<gauss->end();ig++){
    61776177
    … …  
    62456245                if(reCast<bool,IssmDouble>(dt)) scalar = scalar*dt;
    62466246                for(int i=0;i<numdof;i++) pe->values[i]+=scalar*basis[i];
    6247                
     6247
    62486248                /*Transient term*/
    62496249                if(reCast<bool,IssmDouble>(dt)){
    … …  
    63366336        GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    63376337
    6338 
    63396338        /*Use the dof list to index into the solution vector: */
    63406339        for(int i=0;i<numdof;i++){
    … …  
    63526351                this->parameters->FindParam(&penalty_factor,HydrologydcPenaltyFactorEnum);
    63536352                kappa=kmax*pow(10.,penalty_factor);
    6354                
     6353
    63556354                for(int i=0;i<NUMVERTICES;i++){
    63566355                        this->GetHydrologyDCInefficientHmax(&h_max,nodes[i]);
    … …  
    64746473/*FUNCTION Tria::GetHydrologyTransfer{{{*/
    64756474void  Tria::GetHydrologyTransfer(Vector<IssmDouble>* transfer){
    6476        
     6475
    64776476        const int  numdof         = NDOF1 *NUMVERTICES;
    64786477        int        *doflist       = NULL;
    … …  
    64856484        IssmDouble epl_head[numdof],sed_head[numdof];
    64866485
    6487 
    64886486        GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    6489        
     6487
    64906488        /*Get the flag to know if the efficient layer is present*/
    64916489        this->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
    6492        
     6490
    64936491        if(!isefficientlayer){
    64946492                transfer->SetValues(numdof,doflist,&wh_trans[0],ADD_VAL);
    … …  
    64976495        /*Also get the flag to the transfer method (TO DO)*/
    64986496        this->parameters->FindParam(&transfermethod,HydrologydcTransferFlagEnum);
    6499        
     6497
    65006498        /*Switch between the different transfer methods cases(TO DO)*/
    65016499        switch(transfermethod){
    … …  
    65056503        case 1:
    65066504                this->parameters->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    6507                
     6505
    65086506                sed_trans = matpar->GetSedimentTransmitivity();
    65096507                sed_thick = matpar->GetSedimentThickness();
    6510                
     6508
    65116509                GetInputListOnVertices(&sed_head[0],SedimentHeadEnum);
    65126510                GetInputListOnVertices(&epl_head[0],EplHeadEnum);
    6513                
     6511
    65146512                for(int i=0;i<numdof;i++){
    65156513                        if(sed_head[i]>epl_head[i]){
    … …  
    65266524                _error_("no case higher than 1 for the Transfer method");
    65276525        }
    6528        
     6526
    65296527        /*Assign output pointer*/
    65306528        transfer->SetValues(numdof,doflist,&wh_trans[0],ADD_VAL);

  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r15015 r15037  
    14581458#ifdef  _HAVE_DAKOTA_
    14591459void FemModel::DakotaResponsesx(double* d_responses,char** responses_descriptors,int numresponsedescriptors,int d_numresponses){/*{{{*/
    1460        
     1460
    14611461        int        i,j;
    14621462        int        my_rank;
    … …  
    15621562
    15631563  int      i;
    1564        
     1564
    15651565        /*intermediary: */
    15661566        Element *element     = NULL;
    1567        
     1567
    15681568        /*Go through elements, and add contribution from each element to the deflection vector wg:*/
    15691569        for (i=0;i<elements->Size();i++){
    … …  
    15781578
    15791579        Vector<IssmDouble>* basisg=NULL;
    1580        
     1580
    15811581        /*Vector allocation*/
    15821582        basisg=new Vector<IssmDouble>(nodes->NumberOfNodes());
    1583        
     1583
    15841584        for (int i=0;i<elements->Size();i++){
    15851585                Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
    … …  
    15881588        /*Assemble*/
    15891589        basisg->Assemble();
    1590        
     1590
    15911591        /*Update Inputs*/
    15921592        InputUpdateFromVectorx(elements,nodes,vertices,loads,materials,parameters,basisg,BasisIntegralEnum,NodesEnum);
    1593        
     1593
    15941594}
    15951595/*}}}*/
    … …  
    15981598
    15991599        Vector<IssmDouble>* transferg=NULL;
    1600        
     1600
    16011601        /*Vector allocation*/
    16021602        transferg=new Vector<IssmDouble>(nodes->NumberOfNodes());
    1603        
     1603
    16041604        for (int i=0;i<elements->Size();i++){
    16051605                Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
    … …  
    16081608        /*Assemble*/
    16091609        transferg->Assemble();
    1610        
     1610
    16111611        /*Update Inputs*/
    16121612        InputUpdateFromVectorx(elements,nodes,vertices,loads,materials,parameters,transferg,WaterTransferEnum,NodesEnum);
    1613        
    1614 }
    1615 /*}}}*/
     1613
     1614}
     1615/*}}}*/

  • issm/trunk-jpl/src/c/classes/Inputs/Input.h

    r15012 r15037  
    1717#include "../../shared/shared.h"
    1818/*}}}*/
    19 
    2019
    2120class Input: public Object{

  • issm/trunk-jpl/src/c/classes/Inputs/Inputs.cpp

    r15012 r15037  
    428428/*Methods relating to inputs: */
    429429void IsInputConverged(IssmDouble* peps, Input** new_inputs,Input** old_inputs,int num_inputs,int criterion_enum){ /*{{{*/
    430  
     430
    431431        /*figure out convergence at the input level.
    432432          We are given a list of inputs, new and old, and a criterion, and using the

  • issm/trunk-jpl/src/c/classes/Node.cpp

    r15036 r15037  
    926926/*}}}*/
    927927
    928 
    929928/*Methods inherent to Node: */
    930929int* GetLocalDofList(Node** nodes,int numnodes,int setenum,int approximation){ /*{{{*/
    … …  
    987986int* GetGlobalDofList(Node** nodes,int numnodes,int setenum,int approximation){/*{{{*/
    988987
    989 
    990988        int  i,numdof,count;
    991989        int* ndof_list=NULL;
    … …  
    10251023int GetNumberOfDofs(Node** nodes,int numnodes,int setenum,int approximation){/*{{{*/
    10261024
    1027 
    10281025        /*output: */
    10291026        int numberofdofs=0;
    … …  
    10381035#ifdef _HAVE_DIAGNOSTIC_
    10391036void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum){/*{{{*/
    1040 
    10411037
    10421038        int* cs_array=NULL;
    … …  
    10541050/*}}}*/
    10551051void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array){/*{{{*/
    1056 
    10571052
    10581053        int     i,j;
    … …  
    11051100void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int* cs_array){/*{{{*/
    11061101
    1107 
    11081102        int     i;
    11091103        int     numdofs   = 0;
    … …  
    11391133void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int cs_enum){/*{{{*/
    11401134
    1141 
    11421135        int* cs_array=NULL;
    11431136
    … …  
    11541147/*}}}*/
    11551148void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int* cs_array){/*{{{*/
    1156 
    11571149
    11581150        int     i;
    … …  
    11891181void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum){/*{{{*/
    11901182
    1191 
    11921183        int* cs_array=NULL;
    11931184
    … …  
    12041195/*}}}*/
    12051196void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array){/*{{{*/
    1206 
    12071197
    12081198        int     i,j;
    … …  
    12391229/*}}}*/
    12401230void CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array){/*{{{*/
    1241 
    12421231
    12431232        int     i,counter;

  • issm/trunk-jpl/src/c/classes/Nodes.h

    r15012 r15037  
    99class Nodes;
    1010class Materials;
    11 
    1211
    1312/*!\brief Declaration of Nodes class.

  • issm/trunk-jpl/src/c/classes/Params/Parameters.cpp

    r15012 r15037  
    571571         * and parameters
    572572         */
    573 
    574573
    575574        char* options=NULL;

  • issm/trunk-jpl/src/c/classes/Params/Parameters.h

    r15012 r15037  
    6262void        ToolkitsOptionsFromAnalysis(Parameters* parameters,int analysis_type);
    6363
    64 
    6564#endif //ifndef _PARAMETERS_H_

  • issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/GiaDeflectionCorex.cpp

    r15025 r15037  
    1616        double pset[7];
    1717};
    18        
     18
    1919struct blocko{
    2020        double rhoi;
    … …  
    154154        /*dmi: */
    155155        blockt_dmi=xNew<IssmDouble>(Ntimm);
    156        
     156
    157157        /*irate: */
    158158        blockn_.irate=irate;
    … …  
    187187
    188188}
    189      

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.h

    r14996 r15037  
    99
    1010#include "../../classes/classes.h"
    11 
    1211
    1312void ModelProcessorx(Elements** pelements, Nodes** pnodes, Vertices** pvertices, Materials** pmaterials, Constraints** pconstraints, Loads** ploads, Parameters** pparameters, FILE* iomodel_handle,char* rootpath,const int solution_type,const int nummodels,const int* analysis_type_listh);

  • issm/trunk-jpl/src/c/shared/Numerics/extrema.cpp

    r14998 r15037  
    3838}
    3939#endif
    40 
    41 
Note: See TracChangeset for help on using the changeset viewer.