Changeset 15643


Ignore:
Timestamp:
07/27/13 14:25:22 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: cosmetics, removed double blank lines

Location:
issm/trunk-jpl/src
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/Elements/Element.h

    r15613 r15643  
    133133                virtual void HydrologyEPLGetActive(Vector<IssmDouble>* active)=0;
    134134                #endif
    135                
     135
    136136                #ifdef _HAVE_GROUNDINGLINE_
    137137                virtual void   MigrateGroundingLine(IssmDouble* old_floating_ice,IssmDouble* sheet_ungrounding)=0;

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

    r15641 r15643  
    794794        area_init=fabs(xyz_list[1][0]*xyz_list[2][1] - xyz_list[1][1]*xyz_list[2][0] + xyz_list[0][0]*xyz_list[1][1] - xyz_list[0][1]*xyz_list[1][0] + xyz_list[2][0]*xyz_list[0][1] - xyz_list[2][1]*xyz_list[0][0])/2.;
    795795
    796 
    797796        /*Initialize xyz_list with original xyz_list of triangle coordinates*/
    798797        for(j=0;j<3;j++){
    … …  
    23422341                        /*Add input to the element: */
    23432342                        this->inputs->AddInput(new PentaInput(name,values,P1Enum));
    2344                        
     2343
    23452344                        /*Free ressources:*/
    23462345                        xDelete<int>(doflist);
    23472346                        return;
    2348                        
     2347
    23492348          default:
    23502349                        _error_("type " << type << " (" << EnumToStringx(type) << ") not implemented yet");
    … …  
    24222421                                name==HydrologydcMaskEplactiveEnum ||
    24232422                                name==WaterTransferEnum
    2424                                
     2423
    24252424                                ) {
    24262425                return true;
    … …  
    73257324        /*Transform Coordinate System*/
    73267325        TransformStiffnessMatrixCoord(Ke,nodes,NUMVERTICES,XYZPEnum);
    7327        
     7326
    73287327        /*Clean up and return*/
    73297328        delete gauss;
    … …  
    81738172                        _error_("Finite element" << fe_FS << " not supported yet");
    81748173        }
    8175 
    81768174
    81778175        /*clean-up and return*/
    … …  
    99859983
    99869984        if (!IsOnBed()) return;
    9987        
     9985
    99889986        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    99899987        tria->GetHydrologyDCInefficientHmax(ph_max,innode);
    … …  
    99959993
    99969994        if (!IsOnBed()) return;
    9997        
     9995
    99989996        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    99999997        tria->GetHydrologyTransfer(transfer);
    … …  
    1003610034
    1003710035        if (!IsOnBed())return;
    10038        
     10036
    1003910037        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    1004010038        tria->HydrologyEPLGetActive(active_vec);
    … …  
    1004710045
    1004810046        if (!IsOnBed())return;
    10049        
     10047
    1005010048        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    1005110049        tria->HydrologyEPLGetMask(vec_mask);
    1005210050        delete tria->material; delete tria;
    10053        
     10051
    1005410052}
    1005510053/*}}}*/
    … …  
    1008710085                this->parameters->FindParam(&kmax,HydrologySedimentKmaxEnum);
    1008810086                this->parameters->FindParam(&penalty_factor,HydrologydcPenaltyFactorEnum);
    10089                
     10087
    1009010088                kappa=kmax*pow(10.,penalty_factor);
    10091                
     10089
    1009210090                Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.       
    1009310091                for(int i=0;i<NUMVERTICES2D;i++){

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

    r15625 r15643  
    168168                void   InputControlUpdate(IssmDouble scalar,bool save_parameter);
    169169                #endif
    170                
     170
    171171                #ifdef _HAVE_GROUNDINGLINE_
    172172                void   MigrateGroundingLine(IssmDouble* old_floating_ice,IssmDouble* sheet_ungrounding);
    … …  
    316316
    317317                #ifdef _HAVE_HYDROLOGY_
    318                
     318
    319319                ElementMatrix* CreateKMatrixHydrologyDCInefficient(void);
    320320                ElementMatrix* CreateKMatrixHydrologyDCEfficient(void);
    321321                ElementVector* CreatePVectorHydrologyDCInefficient(void);
    322322                ElementVector* CreatePVectorHydrologyDCEfficient(void);
    323                
     323
    324324                void    GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode);
    325325                void    GetHydrologyTransfer(Vector<IssmDouble>* transfer);

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

    r15642 r15643  
    14571457        }
    14581458
    1459 
    14601459}
    14611460/*}}}*/

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

    r15641 r15643  
    2424/*FUNCTION Tria::Tria(){{{*/
    2525Tria::Tria(){
    26        
     26
    2727        int i;
    2828
    … …  
    34833483        xDelete<IssmDouble>(values);
    34843484        xDelete<int>(doflist);
    3485        
    34863485
    34873486}
    … …  
    45414540                                        pow(log((fabs(vx)+epsvel)/(fabs(vxobs)+epsvel)),2) +
    45424541                                        pow(log((fabs(vy)+epsvel)/(fabs(vyobs)+epsvel)),2) );
    4543 
    45444542
    45454543                /*Add to cost function*/
    … …  
    56285626/*FUNCTION Tria::AllActive{{{*/
    56295627bool Tria::AllActive(void){
    5630        
     5628
    56315629        /*Intermediaries*/
    56325630        const int  numnodes = NUMVERTICES;
    5633        
     5631
    56345632        for(int i=0;i<numnodes;i++){
    56355633                if(!this->nodes[i]->IsActive()) return false;
    56365634        }
    5637        
     5635
    56385636        return true;
    56395637}
    … …  
    56415639/*FUNCTION Tria::AnyActive{{{*/
    56425640bool Tria::AnyActive(void){
    5643        
     5641
    56445642        /*Intermediaries*/
    56455643        const int  numnodes = NUMVERTICES;
    5646        
     5644
    56475645        for(int i=0;i<numnodes;i++){
    56485646                if(this->nodes[i]->IsActive()) return true;
    56495647        }
    5650        
     5648
    56515649        return false;
    56525650}/*}}}*/
    … …  
    60966094        }
    60976095        delete gauss;
    6098        
     6096
    60996097        /*      Add residual if necessary*/
    61006098        gauss=new GaussTria();
    … …  
    61136111/*FUNCTION Tria::GetSolutionFromInputsOneDof{{{*/
    61146112void  Tria::GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution, int enum_type){
    6115        
     6113
    61166114        const int    numdof=NDOF1*NUMVERTICES;
    61176115
    … …  
    62006198                this->parameters->FindParam(&kmax,HydrologySedimentKmaxEnum);
    62016199                this->parameters->FindParam(&penalty_factor,HydrologydcPenaltyFactorEnum);
    6202                
     6200
    62036201                kappa=kmax*pow(10.,penalty_factor);
    6204                
     6202
    62056203                for(int i=0;i<NUMVERTICES;i++){
    62066204                        this->GetHydrologyDCInefficientHmax(&h_max,nodes[i]);
    … …  
    62726270
    62736271        GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    6274        
     6272
    62756273        /*Get the flag to know if the efficient layer is present*/
    62766274        this->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
    6277        
     6275
    62786276        if(isefficientlayer){
    62796277                /*Also get the flag to the transfer method*/
    … …  
    62906288                        GetInputListOnVertices(&sed_head[0],SedimentHeadEnum);
    62916289                        GetInputListOnVertices(&epl_head[0],EplHeadEnum);
    6292                        
     6290
    62936291                        this->parameters->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    62946292
    62956293                        sed_trans = matpar->GetSedimentTransmitivity();
    62966294                        sed_thick = matpar->GetSedimentThickness();
    6297                        
     6295
    62986296                        for(int i=0;i<numdof;i++){
    62996297                                /*No transfer if the EPL is not active*/
    … …  
    63046302                                        eplstoring[i]=matpar->GetEplStoring();         
    63056303                                        sedstoring[i]=matpar->GetSedimentStoring();
    6306                                        
     6304
    63076305                                        /*EPL head higher than sediment head, transfer from the epl to the sediment*/
    63086306                                        if(epl_head[i]>sed_head[i]){
    63096307                                                wh_trans=eplstoring[i]*sed_trans*(epl_head[i]-sed_head[i])/(leakage*sed_thick);                         
    6310                                                
     6308
    63116309                                                /*No transfer if the sediment head is allready at the maximum*/
    63126310                                                this->GetHydrologyDCInefficientHmax(&h_max,nodes[i]);
    … …  
    63376335        IssmDouble flag     = 0.;
    63386336        IssmDouble active[numnodes];
    6339        
     6337
    63406338        GetInputListOnVertices(&active[0],HydrologydcMaskEplactiveEnum);
    63416339
    63426340        for(int i=0;i<numnodes;i++) flag+=active[i];
    6343        
     6341
    63446342        if(flag>0.){
    63456343                for(int i=0;i<numnodes;i++){
    … …  
    63506348                /*Do not do anything: at least one node is active for this element but this element is not solved for*/
    63516349        }
    6352        
     6350
    63536351}
    63546352/*}}}*/
    … …  
    63746372        sedheadmin=sedhead[0];
    63756373        for(i=1;i<numdof;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i];
    6376        
     6374
    63776375        for(i=0;i<numdof;i++){
    63786376                /*Activate EPL if residual is >0 */
    … …  
    63806378                        vec_mask->SetValue(nodes[i]->Sid(),1.,INS_VAL);
    63816379                }
    6382                
     6380
    63836381                /*If mask was alread one, keep one*/
    63846382                else if(old_active[i]>0.){

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

    r15559 r15643  
    345345                _printf0_("Core solution number of flops : " << solution_flops  << " Flops\n");
    346346                _printf0_("Core solution memory used     : " << solution_memory << " Bytes\n");
    347                
     347
    348348                /*Add to results: */
    349349                results->AddObject(new GenericExternalResult<IssmDouble>(results->Size()+1, ProfilingSolutionTimeEnum, solution_time, 1, 0));
    … …  
    359359                _printf0_("AD Solution number of flops : " << solution_flops  << " Flops\n");
    360360                _printf0_("AD Solution memory used     : " << solution_memory << " Bytes\n");
    361 
    362361
    363362                #endif
    … …  
    884883                element->CreateKMatrix(Kff,Kfs);
    885884        }
    886        
     885
    887886        for (i=0;i<this->loads->Size();i++){
    888887                load=dynamic_cast<Load*>(this->loads->GetObjectByOffset(i));

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

    r15322 r15643  
    181181                iomodel->data[name]=matrix;
    182182                iomodel->independents[name]=true;
    183                
     183
    184184                //Finally, record the number of independents created:
    185185                this->numberofindependents=M*N;

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

    r15625 r15643  
    288288        _assert_(this->NumberofNodes()==6); //Check Tria too
    289289
    290 
    291290        /*Here, we are computing the strain rate of (vx,0)*/
    292291        for(i=0;i<numnodes;i++){
    … …  
    313312        GetBHO(&B[0][0], xyz_list, gauss);
    314313        _assert_(this->NumberofNodes()==6); //Check Tria too
    315 
    316314
    317315        /*Here, we are computing the strain rate of (0,vy)*/

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

    r15557 r15643  
    6262        if (yposting<0) yflip=1;
    6363        else yflip=0;
    64 
    65 
    6664
    6765        /*Get extreme coordinates of the grid*/

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateLoadsHydrologyDCInefficient.cpp

    r15465 r15643  
    1616        iomodel->Constant(&hydrology_model,HydrologyModelEnum);
    1717        if(hydrology_model!=HydrologydcEnum) return;
    18        
     18
    1919        iomodel->FetchData(1,MeshVertexonbedEnum);
    2020

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/NodesPartitioning.cpp

    r15611 r15643  
    4343         * continuous galerkin, but the nodes partitioning involves faces, which
    4444         * messes up sorting of ids. */
    45 
    4645
    4746        /*Intermediaries*/

  • issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp

    r15412 r15643  
    1919        Vector<IssmDouble>* uf_epl_sub_iter=NULL;
    2020        Vector<IssmDouble>* ug_epl_main_iter=NULL;
    21  
     21
    2222        Vector<IssmDouble>* ys=NULL;
    2323        Vector<IssmDouble>* dug=NULL;
    24        
     24
    2525        Matrix<IssmDouble>* Kff=NULL;
    2626        Matrix<IssmDouble>* Kfs=NULL;
    … …  
    5151        //femmodel->SetCurrentConfiguration(HydrologyDCInefficientAnalysisEnum);
    5252        GetSolutionFromInputsx(&ug_sed, femmodel->elements, femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials, femmodel->parameters);
    53        
     53
    5454        if(isefficientlayer) {
    5555                femmodel->SetCurrentConfiguration(HydrologyDCEfficientAnalysisEnum);
    … …  
    137137                                ConstraintsStatex(&constraints_converged,&num_unstable_constraints, femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
    138138                                femmodel->HydrologyEPLupdateDomainx();                 
    139                                
     139
    140140                                if (!eplconverged){
    141141                                        if(VerboseConvergence()) _printf0_("   # EPL unstable constraints = " << num_unstable_constraints << "\n");
    … …  
    206206                if(hydroconverged)break;
    207207        }
    208        
     208
    209209        InputUpdateFromSolutionx(femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,ug_sed);
    210210        if(isefficientlayer)InputUpdateFromSolutionx(femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,ug_epl);

  • issm/trunk-jpl/src/m/classes/model/model.m

    r15614 r15643  
    12071207                eval(['delete ' id]);
    12081208
    1209 
    12101209                end % }}}
    12111210                function md=download(md) % {{{
    1212                
     1211
    12131212                %the goal of this routine is to download the internals of the current model from a server, because
    12141213                %this model is empty, except for the settings which tell us where to go and find this model!
    1215                
     1214
    12161215                %Download the file:
    12171216                issmscpin(md.settings.upload_server, md.settings.upload_login, md.settings.upload_port, md.settings.upload_path, {md.settings.upload_filename});

  • issm/trunk-jpl/src/m/classes/settings.m

    r15316 r15643  
    4646                        %0 to deactivate
    4747                        obj.waitonlock=Inf;
    48                
     48
    4949                        %upload options:
    5050                        upload_port         = 0;

  • issm/trunk-jpl/src/m/coordsystems/lambert2xy.m

    r15139 r15643  
    6767        y=(B/D)*((cos(b0)*sin(b))-(sin(b0)*cos(b)*cos(lam-lam0)));
    6868end
    69 

  • issm/trunk-jpl/src/m/coordsystems/xy2lambert.m

    r15139 r15643  
    6767        lon=lon_rad*180/pi;
    6868end
    69 

  • issm/trunk-jpl/src/m/plot/plot_icefront.m

    r15612 r15643  
    3535        %Plot zero icelevelset line
    3636
    37 
    3837else
    3938
Note: See TracChangeset for help on using the changeset viewer.