Changeset 15104


Ignore:
Timestamp:
05/24/13 19:11:47 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: simplifying prints

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

Legend:

Unmodified
Added
Removed

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

    r15100 r15104  
    7171
    7272        /*Determine solution sequence: */
    73         if(VerboseQmu()) _printf0_("Starting " << EnumToStringx(solution_type) << " core:" << "\n");
     73        if(VerboseQmu()) _printf0_("Starting " << EnumToStringx(solution_type) << " core:\n");
    7474        WrapperCorePointerFromSolutionEnum(&solutioncore,femmodel->parameters,solution_type,nodakotacore);
    7575
    … …  
    7878
    7979        /*compute responses: */
    80         if(VerboseQmu()) _printf0_("compute dakota responses:" << "\n");
     80        if(VerboseQmu()) _printf0_("compute dakota responses:\n");
    8181        femmodel->DakotaResponsesx(d_responses,responses_descriptors,numresponsedescriptors,d_numresponses);
    8282

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

    r15100 r15104  
    1313        Nodes *nodes = NULL;
    1414
    15         if(VerboseSolution()) _printf0_("   updating boundary conditions..." << "\n");
     15        if(VerboseSolution()) _printf0_("   updating boundary conditions...\n");
    1616
    1717        /*set current analysis: */

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

    r15100 r15104  
    4949                        if(!(num_dependents*num_independents)) return;
    5050
    51                         if(VerboseAutodiff())_printf0_("   start ad core" << "\n");
     51                        if(VerboseAutodiff())_printf0_("   start ad core\n");
    5252
    5353                        /*retrieve state variable: */
    … …  
    286286                        tapestats(1,tape_stats); //reading of tape statistics
    287287                        if(VerboseAutodiff()){
    288                                 _printf0_("   ADOLC statistics: " << "\n");
     288                                _printf0_("   ADOLC statistics: \n");
    289289                                _printf0_("   "<<setw(45)<<left<<"Number of independents: " <<tape_stats[0] << "\n");
    290290                                _printf0_("   "<<setw(45)<<left<<"Number of dependents: " <<tape_stats[1] << "\n");
    … …  
    294294                                _printf0_("   "<<setw(45)<<left<<"Total number of operations recorded: " <<tape_stats[5] << "\n");
    295295                        }
    296                         if(VerboseAutodiff())_printf0_("   end AD core" << "\n");
     296                        if(VerboseAutodiff())_printf0_("   end AD core\n");
    297297
    298298                        /*Free resources: */

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

    r15100 r15104  
    1919
    2020        /*compute thickness */
    21         if(VerboseSolution()) _printf0_("   computing thickness" << "\n");
     21        if(VerboseSolution()) _printf0_("   computing thickness\n");
    2222        femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum);
    2323        solutionsequence_linear(femmodel);
    … …  
    2727
    2828        /*compute adjoint*/
    29         if(VerboseSolution()) _printf0_("   computing adjoint" << "\n");
     29        if(VerboseSolution()) _printf0_("   computing adjoint\n");
    3030        femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum,AdjointBalancethicknessAnalysisEnum);
    3131        solutionsequence_adjoint_linear(femmodel);
    … …  
    3333        /*Save results*/
    3434        if(save_results){
    35                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     35                if(VerboseSolution()) _printf0_("   saving results\n");
    3636                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointEnum);
    3737        }

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

    r15100 r15104  
    2222
    2323        /*Compute velocities*/
    24         if(VerboseSolution()) _printf0_("   computing velocities" << "\n");
     24        if(VerboseSolution()) _printf0_("   computing velocities\n");
    2525        femmodel->SetCurrentConfiguration(DiagnosticHorizAnalysisEnum);
    2626        solutionsequence_nonlinear(femmodel,conserve_loads);
    … …  
    3030
    3131        /*Compute adjoint*/
    32         if(VerboseSolution()) _printf0_("   computing adjoint" << "\n");
     32        if(VerboseSolution()) _printf0_("   computing adjoint\n");
    3333        femmodel->SetCurrentConfiguration(DiagnosticHorizAnalysisEnum,AdjointHorizAnalysisEnum);
    3434        solutionsequence_adjoint_linear(femmodel);
    … …  
    3636        /*Save results*/
    3737        if(save_results){
    38                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     38                if(VerboseSolution()) _printf0_("   saving results\n");
    3939                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointxEnum);
    4040                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointyEnum);

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

    r15100 r15104  
    2121        femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
    2222
    23         if(VerboseSolution()) _printf0_("call computational core:" << "\n");
     23        if(VerboseSolution()) _printf0_("call computational core:\n");
    2424        solutionsequence_linear(femmodel);
    2525
    2626        if(save_results){
    27                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     27                if(VerboseSolution()) _printf0_("   saving results\n");
    2828                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ThicknessEnum);
    2929        }

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

    r15100 r15104  
    1818        femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
    1919
    20         if(VerboseSolution()) _printf0_("   computing slope" << "\n");
     20        if(VerboseSolution()) _printf0_("   computing slope\n");
    2121
    2222        /*Call on core computations: */
    … …  
    2727
    2828        if(save_results){
    29                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     29                if(VerboseSolution()) _printf0_("   saving results\n");
    3030                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BedSlopeXEnum);
    3131                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BedSlopeYEnum);

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

    r15100 r15104  
    6666
    6767        /*Launch once a complete solution to set up all inputs*/
    68         if(VerboseControl()) _printf0_("   preparing initial solution" << "\n");
     68        if(VerboseControl()) _printf0_("   preparing initial solution\n");
    6969        if(isstokes) solutioncore(femmodel);
    7070
    … …  
    8888                if(solution_type==SteadystateSolutionEnum) solutioncore(femmodel);
    8989
    90                 if(VerboseControl()) _printf0_("   compute adjoint state:" << "\n");
     90                if(VerboseControl()) _printf0_("   compute adjoint state:\n");
    9191                adjointcore(femmodel);
    9292                gradient_core(femmodel,n,search_scalar==0);
    … …  
    9898                }
    9999
    100                 if(VerboseControl()) _printf0_("   optimizing along gradient direction" << "\n");
     100                if(VerboseControl()) _printf0_("   optimizing along gradient direction\n");
    101101                optpars.maxiter=reCast<int,IssmDouble>(maxiter[n]); optpars.cm_jump=cm_jump[n];
    102102                BrentSearch(&search_scalar,J+n,&optpars,&objectivefunction,&optargs);
    103103
    104                 if(VerboseControl()) _printf0_("   updating parameter using optimized search scalar" << "\n"); //true means update save controls
     104                if(VerboseControl()) _printf0_("   updating parameter using optimized search scalar\n"); //true means update save controls
    105105                InputControlUpdatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,search_scalar,true);
    106106
    … …  
    108108        }
    109109
    110         if(VerboseControl()) _printf0_("   preparing final solution" << "\n");
     110        if(VerboseControl()) _printf0_("   preparing final solution\n");
    111111        femmodel->parameters->SetParam(true,SaveResultsEnum);
    112112        solutioncore(femmodel);

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

    r15100 r15104  
    5454        /*Initialize TAO*/
    5555        TaoCreate(IssmComm::GetComm(),&tao);
    56         if(VerboseControl()) _printf0_("   Initializing the Toolkit for Advanced Optimization (TAO)" << "\n");
     56        if(VerboseControl()) _printf0_("   Initializing the Toolkit for Advanced Optimization (TAO)\n");
    5757        TaoSetFromOptions(tao);
    5858        TaoSetType(tao,"tao_blmvm");
    … …  
    7979
    8080        /*Solver optimization problem*/
    81         if(VerboseControl()) _printf0_("   Starting optimization" << "\n");
     81        if(VerboseControl()) _printf0_("   Starting optimization\n");
    8282        TaoSolve(tao);
    8383        TaoView(tao,PETSC_VIEWER_STDOUT_WORLD);
    … …  
    9090
    9191        /*Finalize*/
    92         if(VerboseControl()) _printf0_("   preparing final solution" << "\n");
     92        if(VerboseControl()) _printf0_("   preparing final solution\n");
    9393        femmodel->parameters->SetParam(true,SaveResultsEnum);
    9494        void (*solutioncore)(FemModel*)=NULL;
    … …  
    166166
    167167        TaoGetSolutionStatus(tao, &its, &f, &gnorm, &cnorm, &xdiff, NULL);
    168         if(its==0) _printf0_("Iter       Function      Residual  |  List of contributions" << "\n");
    169         if(its==0) _printf0_("-----------------------------------+-----------------------" << "\n");
     168        if(its==0) _printf0_("Iter       Function      Residual  |  List of contributions\n");
     169        if(its==0) _printf0_("-----------------------------------+-----------------------\n");
    170170        _printf0_(setw(4)<<its<<"   "<<setw(12)<<setprecision(7)<<f<<"  "<<setw(12)<<setprecision(7)<<gnorm<<"  | ");
    171171        user->J[its]=f;
    … …  
    176176                _printf0_(" "<<setw(12)<<setprecision(7)<<f);
    177177        }
    178         _printf0_("" << "\n");
     178        _printf0_("\n");
    179179
    180180        /*Clean-up and return*/

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

    r15100 r15104  
    5858        if(ishutter){
    5959
    60                 if(VerboseSolution()) _printf0_("   computing hutter velocities" << "\n");
     60                if(VerboseSolution()) _printf0_("   computing hutter velocities\n");
    6161
    6262                //Take the last velocity into account so that the velocity on the MacAyeal domain is not zero
    … …  
    7171        if ((ismacayealpattyn || isl1l2) ^ isstokes){ // ^ = xor
    7272
    73                 if(VerboseSolution()) _printf0_("   computing velocities" << "\n");
     73                if(VerboseSolution()) _printf0_("   computing velocities\n");
    7474                femmodel->SetCurrentConfiguration(DiagnosticHorizAnalysisEnum);
    7575                if(newton>0)
    … …  
    8181        if (ismacayealpattyn && isstokes){
    8282
    83                 if(VerboseSolution()) _printf0_("   computing coupling macayealpattyn and stokes velocities and pressure " << "\n");
     83                if(VerboseSolution()) _printf0_("   computing coupling macayealpattyn and stokes velocities and pressure \n");
    8484                solutionsequence_stokescoupling_nonlinear(femmodel,conserve_loads);
    8585        }
    … …  
    8787        if (dim==3 & (ishutter || ismacayealpattyn)){
    8888
    89                 if(VerboseSolution()) _printf0_("   computing vertical velocities" << "\n");
     89                if(VerboseSolution()) _printf0_("   computing vertical velocities\n");
    9090                femmodel->SetCurrentConfiguration(DiagnosticVertAnalysisEnum);
    9191                solutionsequence_linear(femmodel);
    … …  
    9393
    9494        if(save_results){
    95                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     95                if(VerboseSolution()) _printf0_("   saving results\n");
    9696                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VxEnum);
    9797                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VyEnum);

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

    r15100 r15104  
    1818        femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
    1919
    20         if(VerboseSolution()) _printf0_("   computing enthalpy" << "\n");
     20        if(VerboseSolution()) _printf0_("   computing enthalpy\n");
    2121        femmodel->SetCurrentConfiguration(EnthalpyAnalysisEnum);
    2222        solutionsequence_nonlinear(femmodel,true);
    … …  
    2626
    2727        if(save_results){
    28                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     28                if(VerboseSolution()) _printf0_("   saving results\n");
    2929                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum);
    3030                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,EnthalpyEnum);

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

    r15100 r15104  
    2626        femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
    2727
    28         if(VerboseSolution()) _printf0_("   computing GIA" << "\n");
     28        if(VerboseSolution()) _printf0_("   computing GIA\n");
    2929
    3030        /*Call on core computations: */
    … …  
    5050
    5151        if(save_results){
    52                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     52                if(VerboseSolution()) _printf0_("   saving results\n");
    5353                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,GiaWEnum);
    5454                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,GiadWdtEnum);

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

    r15100 r15104  
    2020
    2121        /*Compute gradient*/
    22         if(VerboseControl()) _printf0_("   compute cost function gradient" << "\n");
     22        if(VerboseControl()) _printf0_("   compute cost function gradient\n");
    2323        Gradjx(&gradient,&norm_list,femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters);
    2424
    2525        if (orthogonalize){
    26                 if(VerboseControl()) _printf0_("   orthogonalization" << "\n");
     26                if(VerboseControl()) _printf0_("   orthogonalization\n");
    2727                ControlInputGetGradientx(&old_gradient,femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters);
    2828                Orthx(&new_gradient,gradient,old_gradient); delete old_gradient; delete gradient;

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

    r15100 r15104  
    5858
    5959                if (hydrology_model==HydrologyshreveEnum){
    60                         if(VerboseSolution()) _printf0_("   computing water column" << "\n");
     60                        if(VerboseSolution()) _printf0_("   computing water column\n");
    6161                        femmodel->SetCurrentConfiguration(HydrologyShreveAnalysisEnum);
    6262                        solutionsequence_nonlinear(femmodel,modify_loads);
    … …  
    6767
    6868                        if(save_results && ((i+1)%output_frequency==0 || (i+1)==nsteps)){
    69                                 if(VerboseSolution()) _printf0_("   saving results " << "\n");
     69                                if(VerboseSolution()) _printf0_("   saving results \n");
    7070                                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,WatercolumnEnum);
    7171                                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,HydrologyWaterVxEnum);
    … …  
    7373
    7474                                /*unload results*/
    75                                 if(VerboseSolution()) _printf0_("   saving temporary results" << "\n");
     75                                if(VerboseSolution()) _printf0_("   saving temporary results\n");
    7676                                OutputResultsx(femmodel->elements, femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,femmodel->results);
    7777                        }
    … …  
    8282                        femmodel->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
    8383
    84                         if(VerboseSolution()) _printf0_("   computing water transfer" << "\n");
     84                        if(VerboseSolution()) _printf0_("   computing water transfer\n");
    8585
    86                         if(VerboseSolution()) _printf0_("   computing water head" << "\n");
     86                        if(VerboseSolution()) _printf0_("   computing water head\n");
    8787                        solutionsequence_hydro_nonlinear(femmodel);
    8888                        if(save_results && ((i+1)%output_frequency==0 || (i+1)==nsteps)){
    89                                 if(VerboseSolution()) _printf0_("   saving results " << "\n");
     89                                if(VerboseSolution()) _printf0_("   saving results \n");
    9090                                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SedimentHeadEnum);
    9191                                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SedimentHeadResidualEnum);
    … …  
    9494                                }
    9595                                /*unload results*/
    96                                 if(VerboseSolution()) _printf0_("   saving temporary results" << "\n");
     96                                if(VerboseSolution()) _printf0_("   saving temporary results\n");
    9797                                OutputResultsx(femmodel->elements, femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,femmodel->results);
    9898                        }

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

    r15100 r15104  
    3232
    3333        if(issmbgradients){
    34           if(VerboseSolution)_printf_(" call smb gradients module\n" << "\n");
     34          if(VerboseSolution)_printf_(" call smb gradients module\n\n");
    3535          SmbGradientsx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
    3636        }
    3737        if(ispdd){
    3838                if(isdelta18o){
    39                         if(VerboseSolution()) _printf0_("   call Delta18oParametrization module" << "\n");
     39                        if(VerboseSolution()) _printf0_("   call Delta18oParametrization module\n");
    4040                        Delta18oParameterizationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
    4141                }
    42                 if(VerboseSolution()) _printf0_("   call positive degree day module" << "\n");
     42                if(VerboseSolution()) _printf0_("   call positive degree day module\n");
    4343                PositiveDegreeDayx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
    4444        }
    45         if(VerboseSolution()) _printf0_("   call computational core" << "\n");
     45        if(VerboseSolution()) _printf0_("   call computational core\n");
    4646        solutionsequence_linear(femmodel);
    4747
    4848        if(save_results){
    49                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     49                if(VerboseSolution()) _printf0_("   saving results\n");
    5050                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ThicknessEnum);
    5151                femmodel->RequestedOutputsx(requested_outputs,numoutputs);

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

    r15100 r15104  
    6363                #endif
    6464
    65                 if(VerboseSolution()) _printf0_("   computing new velocity" << "\n");
     65                if(VerboseSolution()) _printf0_("   computing new velocity\n");
    6666                diagnostic_core(femmodel);
    6767                GetSolutionFromInputsx(&ug,femmodel->elements, femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials, femmodel->parameters);
    6868
    6969                if(step>1){
    70                         if(VerboseSolution()) _printf0_("   checking steadystate convergence" << "\n");
     70                        if(VerboseSolution()) _printf0_("   checking steadystate convergence\n");
    7171                        if(steadystateconvergence(tg,tg_old,ug,ug_old,reltol)) break;
    7272                }
    7373                if(step>maxiter){
    74                         if(VerboseSolution()) _printf0_("   maximum number steadystate iterations " << maxiter << " reached" << "\n");
     74                        if(VerboseSolution()) _printf0_("   maximum number steadystate iterations " << maxiter << " reached\n");
    7575                        break;
    7676                }
    … …  
    8383
    8484        if(save_results){
    85                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     85                if(VerboseSolution()) _printf0_("   saving results\n");
    8686                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VxEnum);
    8787                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VyEnum);
    … …  
    117117        if (xIsNan<IssmDouble>(ndu) || xIsNan<IssmDouble>(nu)) _error_("convergence criterion is NaN!");
    118118        if((ndu/nu)<reltol){
    119                 if(VerboseConvergence()) _printf0_("\n"<<setw(50)<<left<<"   Velocity convergence: norm(du)/norm(u)"<<ndu/nu*100<<" < "<<reltol*100<<" %" << "\n");
     119                if(VerboseConvergence()) _printf0_("\n"<<setw(50)<<left<<"   Velocity convergence: norm(du)/norm(u)"<<ndu/nu*100<<" < "<<reltol*100<<" %\n");
    120120        }
    121121        else{
    122                 if(VerboseConvergence()) _printf0_("\n"<<setw(50)<<left<<"   Velocity convergence: norm(du)/norm(u)"<<ndu/nu*100<<" > "<<reltol*100<<" %" << "\n");
     122                if(VerboseConvergence()) _printf0_("\n"<<setw(50)<<left<<"   Velocity convergence: norm(du)/norm(u)"<<ndu/nu*100<<" > "<<reltol*100<<" %\n");
    123123                converged=false;
    124124        }
    … …  
    129129        if (xIsNan<IssmDouble>(ndt) || xIsNan<IssmDouble>(nt)) _error_("convergence criterion is NaN!");
    130130        if((ndt/nt)<reltol){
    131                 if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   Temperature convergence: norm(dt)/norm(t)"<<ndt/nt*100<<" < "<<reltol*100<<" %" << "\n");
     131                if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   Temperature convergence: norm(dt)/norm(t)"<<ndt/nt*100<<" < "<<reltol*100<<" %\n");
    132132        }
    133133        else{
    134                 if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   Temperature convergence: norm(dt)/norm(t)"<<ndt/nt*100<<" > "<<reltol*100<<" %" << "\n");
     134                if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   Temperature convergence: norm(dt)/norm(t)"<<ndt/nt*100<<" > "<<reltol*100<<" %\n");
    135135                converged=false;
    136136        }

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

    r15100 r15104  
    1818        femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
    1919
    20         if(VerboseSolution()) _printf0_("computing slope..." << "\n");
     20        if(VerboseSolution()) _printf0_("computing slope...\n");
    2121
    2222        /*Call on core computations: */
    … …  
    2727
    2828        if(save_results){
    29                 if(VerboseSolution()) _printf0_("saving results:" << "\n");
     29                if(VerboseSolution()) _printf0_("saving results:\n");
    3030                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SurfaceSlopeXEnum);
    3131                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SurfaceSlopeYEnum);

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

    r15100 r15104  
    2525        }
    2626
    27         if(VerboseSolution()) _printf0_("   computing temperatures" << "\n");
     27        if(VerboseSolution()) _printf0_("   computing temperatures\n");
    2828        femmodel->SetCurrentConfiguration(ThermalAnalysisEnum);
    2929        solutionsequence_thermal_nonlinear(femmodel);
    3030
    31         if(VerboseSolution()) _printf0_("   computing melting" << "\n");
     31        if(VerboseSolution()) _printf0_("   computing melting\n");
    3232        femmodel->SetCurrentConfiguration(MeltingAnalysisEnum);
    3333        solutionsequence_linear(femmodel);
    3434
    3535        if(save_results){
    36                 if(VerboseSolution()) _printf0_("   saving results" << "\n");
     36                if(VerboseSolution()) _printf0_("   saving results\n");
    3737                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum);
    3838                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BasalforcingsMeltingRateEnum);

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

    r15100 r15104  
    8585                femmodel->parameters->SetParam(step,StepEnum);
    8686
    87                 if(VerboseSolution()) _printf0_("iteration " << step << "/" << floor((finaltime-time)/dt)+step << "  time [yr]: " << time/yts << " (time step: " << dt/yts << ")" << "\n");
     87                if(VerboseSolution()) _printf0_("iteration " << step << "/" << floor((finaltime-time)/dt)+step << "  time [yr]: " << time/yts << " (time step: " << dt/yts << ")\n");
    8888                if(step%output_frequency==0 || time==finaltime)
    8989                 save_results=true;
    … …  
    9393
    9494                if(isthermal && dim==3){
    95                         if(VerboseSolution()) _printf0_("   computing temperatures" << "\n");
     95                        if(VerboseSolution()) _printf0_("   computing temperatures\n");
    9696                        #ifdef _HAVE_THERMAL_
    9797                        if(isenthalpy==0){
    … …  
    107107
    108108                if(isdiagnostic){
    109                         if(VerboseSolution()) _printf0_("   computing new velocity" << "\n");
     109                        if(VerboseSolution()) _printf0_("   computing new velocity\n");
    110110                        #ifdef _HAVE_DIAGNOSTIC_
    111111                        diagnostic_core(femmodel);
    … …  
    116116
    117117                if(isprognostic){
    118                         if(VerboseSolution()) _printf0_("   computing new thickness" << "\n");
     118                        if(VerboseSolution()) _printf0_("   computing new thickness\n");
    119119                        prognostic_core(femmodel);
    120                         if(VerboseSolution()) _printf0_("   updating vertices positions" << "\n");
     120                        if(VerboseSolution()) _printf0_("   updating vertices positions\n");
    121121                        femmodel->UpdateVertexPositionsx();
    122122                }
    123123
    124124                if(isgroundingline){
    125                         if(VerboseSolution()) _printf0_("   computing new grounding line position" << "\n");
     125                        if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
    126126                        #ifdef _HAVE_GROUNDINGLINE_
    127127                        GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
    … …  
    131131                }
    132132                if(isgia){
    133                         if(VerboseSolution()) _printf0_("   computing glacial isostatic adjustment" << "\n");
     133                        if(VerboseSolution()) _printf0_("   computing glacial isostatic adjustment\n");
    134134                        #ifdef _HAVE_GIA_
    135135                        gia_core(femmodel);
    … …  
    142142                /*unload results*/
    143143                if(save_results){
    144                         if(VerboseSolution()) _printf0_("   saving transient results" << "\n");
     144                        if(VerboseSolution()) _printf0_("   saving transient results\n");
    145145                        InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SurfaceEnum);
    146146                        InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,BedEnum);
    … …  
    158158                        }
    159159
    160                         if(VerboseSolution()) _printf0_("   saving temporary results" << "\n");
     160                        if(VerboseSolution()) _printf0_("   saving temporary results\n");
    161161                        OutputResultsx(femmodel->elements, femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,femmodel->results);
    162162                }

  • issm/trunk-jpl/src/c/bamg/BamgVertex.cpp

    r15100 r15104  
    1515        void BamgVertex::Echo(void){
    1616
    17                 _printf_("Vertex:" << "\n");
     17                _printf_("Vertex:\n");
    1818                _printf_("  integer   coordinates i.x: " << i.x << ", i.y: " << i.y << "\n");
    1919                _printf_("  Euclidean coordinates r.x: " << r.x << ", r.y: " << r.y << "\n");

  • issm/trunk-jpl/src/c/bamg/Edge.cpp

    r15100 r15104  
    2727        /*FUNCTION Edge::Echo {{{*/
    2828        void Edge::Echo(void){
    29                 _printf_("Edge:" << "\n");
     29                _printf_("Edge:\n");
    3030                _printf_("   pointers towards two vertices: " << v[0] << " " << v[1] << "\n");
    3131                _printf_("   ReferenceNumber = " << ReferenceNumber << "\n");

  • issm/trunk-jpl/src/c/bamg/EigenMetric.cpp

    r15100 r15104  
    114114        void EigenMetric::Echo(void){
    115115
    116                 _printf_("EigenMetric:" << "\n");
     116                _printf_("EigenMetric:\n");
    117117                _printf_("   lambda1: " << lambda1 << "\n");
    118118                _printf_("   lambda2: " << lambda2 << "\n");

  • issm/trunk-jpl/src/c/bamg/Geometry.cpp

    r15100 r15104  
    8080                //Vertices
    8181                if (bamggeom->Vertices){
    82                         if(verbose>5) _printf_("      processing Vertices" << "\n");
     82                        if(verbose>5) _printf_("      processing Vertices\n");
    8383                        if (bamggeom->VerticesSize[1]!=3) _error_("Vertices should have 3 columns");
    8484                        vertices = new GeomVertex[nbv];
    … …  
    125125                        double* verticeslength=NULL;
    126126
    127                         if(verbose>5) _printf_("      processing Edges" << "\n");
     127                        if(verbose>5) _printf_("      processing Edges\n");
    128128                        if (bamggeom->EdgesSize[1]!=3) _error_("Edges should have 3 columns");
    129129                        edges = new GeomEdge[nbe];
    … …  
    180180                //hVertices
    181181                if(bamgopts->hVertices && bamgopts->hVerticesSize[0]==nbv){
    182                         if(verbose>5) _printf_("      processing hVertices" << "\n");
     182                        if(verbose>5) _printf_("      processing hVertices\n");
    183183                        for (i=0;i< nbv;i++){
    184184                                if (!xIsNan<IssmPDouble>(bamgopts->hVertices[i])){
    … …  
    190190                //MetricVertices
    191191                if(bamgopts->metric && bamgopts->metric[0]==nbv){
    192                         if(verbose>5) _printf_("      processing MetricVertices" << "\n");
     192                        if(verbose>5) _printf_("      processing MetricVertices\n");
    193193                        for (i=0;i< nbv;i++) {
    194194                                vertices[i].m = Metric((double)bamgopts->metric[i*3+0],(double)bamgopts->metric[i*3+1],(double)bamgopts->metric[i*3+2]);
    … …  
    198198                //MaxCornerAngle
    199199                if (bamgopts->MaxCornerAngle){
    200                         if(verbose>5) _printf_("      processing MaxCornerAngle" << "\n");
     200                        if(verbose>5) _printf_("      processing MaxCornerAngle\n");
    201201                        MaxCornerAngle=bamgopts->MaxCornerAngle*Pi/180;
    202202                }
    … …  
    237237                //RequiredVertices
    238238                if(bamggeom->RequiredVertices){
    239                         if(verbose>5) _printf_("      processing RequiredVertices" << "\n");
     239                        if(verbose>5) _printf_("      processing RequiredVertices\n");
    240240                        if (bamggeom->RequiredVerticesSize[1]!=1) _error_("RequiredVertices should have 1 column");
    241241                        n=bamggeom->RequiredVerticesSize[0];
    … …  
    249249                //RequiredEdges
    250250                if(bamggeom->RequiredEdges){
    251                         if(verbose>5) _printf_("      processing RequiredEdges" << "\n");
     251                        if(verbose>5) _printf_("      processing RequiredEdges\n");
    252252                        if (bamggeom->RequiredEdgesSize[1]!=1) _error_("RequiredEdges should have 1 column");
    253253                        n=bamggeom->RequiredEdgesSize[0];
    … …  
    261261                //SubDomain
    262262                if(bamggeom->SubDomains){
    263                         if(verbose>5) _printf_("      processing SubDomains" << "\n");
     263                        if(verbose>5) _printf_("      processing SubDomains\n");
    264264                        if (bamggeom->SubDomainsSize[1]!=4) _error_("SubDomains should have 4 columns");
    265265                        nbsubdomains=bamggeom->SubDomainsSize[0];
    … …  
    293293
    294294                /*Vertices*/
    295                 if(verbose>5) _printf_("      writing Vertices" << "\n");
     295                if(verbose>5) _printf_("      writing Vertices\n");
    296296                bamggeom->VerticesSize[0]=nbv;
    297297                bamggeom->VerticesSize[1]=3;
    … …  
    309309
    310310                /*Edges*/
    311                 if(verbose>5) _printf_("      writing Edges" << "\n");
     311                if(verbose>5) _printf_("      writing Edges\n");
    312312                bamggeom->EdgesSize[0]=nbe;
    313313                bamggeom->EdgesSize[1]=3;
    … …  
    327327
    328328                /*RequiredEdges*/
    329                 if(verbose>5) _printf_("      writing " << nbreq << " RequiredEdges" << "\n");
     329                if(verbose>5) _printf_("      writing " << nbreq << " RequiredEdges\n");
    330330                bamggeom->RequiredEdgesSize[0]=nbreq;
    331331                bamggeom->RequiredEdgesSize[1]=1;
    … …  
    344344
    345345                /*RequiredVertices*/
    346                 if(verbose>5) _printf_("      writing " << nbreqv << " RequiredVertices" << "\n");
     346                if(verbose>5) _printf_("      writing " << nbreqv << " RequiredVertices\n");
    347347                bamggeom->RequiredVerticesSize[0]=nbreqv;
    348348                bamggeom->RequiredVerticesSize[1]=1;
    … …  
    359359
    360360                /*SubDomains*/
    361                 if(verbose>5) _printf_("      writing SubDomains" << "\n");
     361                if(verbose>5) _printf_("      writing SubDomains\n");
    362362                bamggeom->SubDomainsSize[0]=nbsubdomains;
    363363                bamggeom->SubDomainsSize[1]=4;
    … …  
    373373
    374374                /*TangentAtEdges*/
    375                 if(verbose>5) _printf_("      writing TangentAtEdges" << "\n");
     375                if(verbose>5) _printf_("      writing TangentAtEdges\n");
    376376                bamggeom->TangentAtEdgesSize[0]=nbtan;
    377377                bamggeom->TangentAtEdgesSize[1]=4;
    … …  
    402402        void Geometry::Echo(void){
    403403
    404                 _printf_("Geometry:" << "\n");
     404                _printf_("Geometry:\n");
    405405                _printf_("   nbv  (number of vertices) : " << nbv << "\n");
    406406                _printf_("   nbe  (number of edges)    : " << nbe << "\n");
    … …  
    412412                _printf_("   subdomains: " << subdomains << "\n");
    413413                _printf_("   curves: " << curves << "\n");
    414                 _printf_("   pmin (x,y): (" << pmin.x << " " << pmin.y << ")" << "\n");
    415                 _printf_("   pmax (x,y): (" << pmax.x << " " << pmax.y << ")" << "\n");
     414                _printf_("   pmin (x,y): (" << pmin.x << " " << pmin.y << ")\n");
     415                _printf_("   pmax (x,y): (" << pmax.x << " " << pmax.y << ")\n");
    416416                _printf_("   coefIcoor: " << coefIcoor << "\n");
    417417                _printf_("   MaxCornerAngle: " << MaxCornerAngle << "\n");
    … …  
    526526                                _printf_("reference numbers: " << v->ReferenceNumber << " " << vertices[i].ReferenceNumber << "\n");
    527527                                _printf_("Id: " << i+1 << "\n");
    528                                 _printf_("Coords: ["<<v->r.x<<" "<<v->r.y<<"] ["<<vertices[i].r.x<<" "<<vertices[i].r.y<<"]" << "\n");
     528                                _printf_("Coords: ["<<v->r.x<<" "<<v->r.y<<"] ["<<vertices[i].r.x<<" "<<vertices[i].r.y<<"]\n");
    529529
    530530                                delete [] next_p;
    … …  
    869869                        if (bge<=0) {
    870870                                if(NbTry) {
    871                                         _printf_("Fatal Error: on the class Mesh before call Geometry::ProjectOnCurve" << "\n");
    872                                         _printf_("That bug might come from:" << "\n");
    873                                         _printf_(" 1)  a mesh edge  containing more than " << mxe/2 << " geometrical edges" << "\n");
    874                                         _printf_(" 2)  code bug : be sure that we call   Mesh::SetVertexFieldOn() before" << "\n");
    875                                         _printf_("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)" << "\n");
     871                                        _printf_("Fatal Error: on the class Mesh before call Geometry::ProjectOnCurve\n");
     872                                        _printf_("That bug might come from:\n");
     873                                        _printf_(" 1)  a mesh edge  containing more than " << mxe/2 << " geometrical edges\n");
     874                                        _printf_(" 2)  code bug : be sure that we call   Mesh::SetVertexFieldOn() before\n");
     875                                        _printf_("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)\n");
    876876                                        _error_("see above");
    877877                                }
    … …  
    886886                while (eg1 != (GeomEdge*) vg1  &&  (*eg1)(direction1) != (GeomVertex*) vg1) {
    887887                        if(tge>=mxe ) {
    888                                 _printf_("WARNING: on the class Mesh before call Geometry::ProjectOnCurve is having issues (isn't it Eric?)" << "\n");
     888                                _printf_("WARNING: on the class Mesh before call Geometry::ProjectOnCurve is having issues (isn't it Eric?)\n");
    889889                                NbTry++;
    890890                                if (NbTry<2) goto retry;
    891                                 _printf_("Fatal Error: on the class Mesh before call Geometry::ProjectOnCurve" << "\n");
    892                                 _printf_("That bug might come from:" << "\n");
    893                                 _printf_(" 1)  a mesh edge  contening more than " << mxe/2 << " geometrical edges" << "\n");
    894                                 _printf_(" 2)  code bug : be sure that we call   Mesh::SetVertexFieldOn() before" << "\n");
    895                                 _printf_("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)" << "\n");
     891                                _printf_("Fatal Error: on the class Mesh before call Geometry::ProjectOnCurve\n");
     892                                _printf_("That bug might come from:\n");
     893                                _printf_(" 1)  a mesh edge  contening more than " << mxe/2 << " geometrical edges\n");
     894                                _printf_(" 2)  code bug : be sure that we call   Mesh::SetVertexFieldOn() before\n");
     895                                _printf_("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)\n");
    896896                                _error_("see above");
    897897                        }

  • issm/trunk-jpl/src/c/bamg/Mesh.cpp

    r15100 r15104  
    3131                if(bamggeom->Edges==NULL) {
    3232                        /*Recreate geometry if needed*/
    33                         _printf_("WARNING: mesh present but no geometry found. Reconstructing..." << "\n");
     33                        _printf_("WARNING: mesh present but no geometry found. Reconstructing...\n");
    3434                        BuildGeometryFromMesh(bamgopts);
    3535                        Gh.PostRead();
    … …  
    281281
    282282                //Vertices
    283                 if (verbose) _printf_("Reading vertices (" << nbv << ")" << "\n");
     283                if (verbose) _printf_("Reading vertices (" << nbv << ")\n");
    284284                vertices=xNew<BamgVertex>(nbv);
    285285                orderedvertices=xNew<BamgVertex*>(nbv);
    … …  
    295295
    296296                //Triangles
    297                 if (verbose) _printf_("Reading triangles (" << nbt << ")" << "\n");
     297                if (verbose) _printf_("Reading triangles (" << nbt << ")\n");
    298298                triangles =new Triangle[maxnbt]; //we cannot allocate only nbt triangles since
    299299                nodeflags=xNew<bool>(nbv);
    … …  
    311311
    312312                /*Recreate geometry: */
    313                 if (verbose) _printf_("Building Geometry" << "\n");
     313                if (verbose) _printf_("Building Geometry\n");
    314314                BuildGeometryFromMesh();
    315                 if (verbose) _printf_("Completing geometry" << "\n");
     315                if (verbose) _printf_("Completing geometry\n");
    316316                Gh.PostRead();
    317317
    … …  
    320320                for(i=0;i<nbv;i++){
    321321                        if(!nodeflags[i]){
    322                                 _printf_("Vertex " << i+1 << " does not belong to any element" << "\n");
     322                                _printf_("Vertex " << i+1 << " does not belong to any element\n");
    323323                                isorphan=true;
    324324                        }
    … …  
    347347                //Vertices
    348348                if(bamgmesh->Vertices){
    349                         if(verbose>5) _printf_("      processing Vertices" << "\n");
     349                        if(verbose>5) _printf_("      processing Vertices\n");
    350350
    351351                        vertices=xNew<BamgVertex>(nbv);
    … …  
    368368                //Triangles
    369369                if(bamgmesh->Triangles){
    370                         if(verbose>5) _printf_("      processing Triangles" << "\n");
     370                        if(verbose>5) _printf_("      processing Triangles\n");
    371371                        triangles =new Triangle[maxnbt]; //we cannot allocate only nbt triangles since
    372372                        //other triangles will be added for each edge
    … …  
    386386                //Quadrilaterals
    387387                if(bamgmesh->Quadrilaterals){
    388                         if(verbose>5) _printf_("      processing Quadrilaterals" << "\n");
     388                        if(verbose>5) _printf_("      processing Quadrilaterals\n");
    389389                        long i1,i2,i3,i4;
    390390                        triangles =new Triangle[nbt];
    … …  
    408408                //VerticesOnGeomEdge
    409409                if(bamgmesh->VerticesOnGeomEdge){
    410                         if(verbose>5) _printf_("      processing VerticesOnGeomEdge" << "\n");
     410                        if(verbose>5) _printf_("      processing VerticesOnGeomEdge\n");
    411411                        NbVerticesOnGeomEdge=bamgmesh->VerticesOnGeomEdgeSize[0];
    412412                        VerticesOnGeomEdge= new  VertexOnGeom[NbVerticesOnGeomEdge] ;
    … …  
    423423                //VerticesOnGeomVertex
    424424                if(bamgmesh->VerticesOnGeomVertexSize[0]){
    425                         if(verbose>5) _printf_("      processing VerticesOnGeomVertex" << "\n");
     425                        if(verbose>5) _printf_("      processing VerticesOnGeomVertex\n");
    426426                        NbVerticesOnGeomVertex=bamgmesh->VerticesOnGeomVertexSize[0];
    427427                        VerticesOnGeomVertex  = new  VertexOnGeom[NbVerticesOnGeomVertex] ;
    … …  
    439439                        double* len=NULL;
    440440
    441                         if(verbose>5) _printf_("      processing Edges" << "\n");
     441                        if(verbose>5) _printf_("      processing Edges\n");
    442442                        nbe=bamgmesh->EdgesSize[0];
    443443                        edges= new Edge[nbe];
    … …  
    499499                //EdgeOnGeomEdge
    500500                if(bamgmesh->EdgesOnGeomEdge){
    501                         if(verbose>5) _printf_("      processing EdgesOnGeomEdge" << "\n");
     501                        if(verbose>5) _printf_("      processing EdgesOnGeomEdge\n");
    502502                        int i1,i2,i,j;
    503503                        i2=bamgmesh->EdgesOnGeomEdgeSize[0];
    … …  
    516516                if(bamgmesh->SubDomains){
    517517                        long i3,head,direction;
    518                         if(verbose>5) _printf_("      processing SubDomains" << "\n");
     518                        if(verbose>5) _printf_("      processing SubDomains\n");
    519519                        nbsubdomains=bamgmesh->SubDomainsSize[0];
    520520                        subdomains = new SubDomain [ nbsubdomains ];
    … …  
    585585
    586586                /*Vertices*/
    587                 if(verbose>5) _printf_("      writing Vertices" << "\n");
     587                if(verbose>5) _printf_("      writing Vertices\n");
    588588                bamgmesh->VerticesSize[0]=nbv;
    589589                bamgmesh->VerticesSize[1]=3;
    … …  
    598598
    599599                /*Edges*/
    600                 if(verbose>5) _printf_("      writing Edges" << "\n");
     600                if(verbose>5) _printf_("      writing Edges\n");
    601601                bamgmesh->EdgesSize[0]=nbe;
    602602                bamgmesh->EdgesSize[1]=3;
    … …  
    615615
    616616                /*Element edges*/
    617                 if(verbose>5) _printf_("      writing element edges" << "\n");
     617                if(verbose>5) _printf_("      writing element edges\n");
    618618                SetOfEdges4* edge4=new SetOfEdges4(nbt*3,nbv);
    619619                double* elemedge=NULL;
    … …  
    671671
    672672                /*IssmSegments*/
    673                 if(verbose>5) _printf_("      writing IssmSegments" << "\n");
     673                if(verbose>5) _printf_("      writing IssmSegments\n");
    674674                bamgmesh->IssmSegmentsSize[0]=NumIssmSegments;
    675675                bamgmesh->IssmSegmentsSize[1]=4;
    … …  
    714714
    715715                /*Triangles*/
    716                 if(verbose>5) _printf_("      writing Triangles" << "\n");
     716                if(verbose>5) _printf_("      writing Triangles\n");
    717717                k=nbInT-nbq*2;
    718718                num=0;
    … …  
    735735
    736736                /*Quadrilaterals*/
    737                 if(verbose>5) _printf_("      writing Quadrilaterals" << "\n");
     737                if(verbose>5) _printf_("      writing Quadrilaterals\n");
    738738                bamgmesh->QuadrilateralsSize[0]=nbq;
    739739                bamgmesh->QuadrilateralsSize[1]=5;
    … …  
    756756
    757757                /*SubDomains*/
    758                 if(verbose>5) _printf_("      writing SubDomains" << "\n");
     758                if(verbose>5) _printf_("      writing SubDomains\n");
    759759                bamgmesh->SubDomainsSize[0]=nbsubdomains;
    760760                bamgmesh->SubDomainsSize[1]=4;
    … …  
    770770
    771771                /*SubDomainsFromGeom*/
    772                 if(verbose>5) _printf_("      writing SubDomainsFromGeom" << "\n");
     772                if(verbose>5) _printf_("      writing SubDomainsFromGeom\n");
    773773                bamgmesh->SubDomainsFromGeomSize[0]=Gh.nbsubdomains;
    774774                bamgmesh->SubDomainsFromGeomSize[1]=4;
    … …  
    784784
    785785                /*VerticesOnGeomVertex*/
    786                 if(verbose>5) _printf_("      writing VerticesOnGeomVertex" << "\n");
     786                if(verbose>5) _printf_("      writing VerticesOnGeomVertex\n");
    787787                bamgmesh->VerticesOnGeomVertexSize[0]=NbVerticesOnGeomVertex;
    788788                bamgmesh->VerticesOnGeomVertexSize[1]=2;
    … …  
    798798
    799799                /*VertexOnGeomEdge*/
    800                 if(verbose>5) _printf_("      writing VerticesOnGeomEdge" << "\n");
     800                if(verbose>5) _printf_("      writing VerticesOnGeomEdge\n");
    801801                bamgmesh->VerticesOnGeomEdgeSize[0]=NbVerticesOnGeomEdge;
    802802                bamgmesh->VerticesOnGeomEdgeSize[1]=3;
    … …  
    815815
    816816                /*EdgesOnGeomEdge*/
    817                 if(verbose>5) _printf_("      writing EdgesOnGeomEdge" << "\n");
     817                if(verbose>5) _printf_("      writing EdgesOnGeomEdge\n");
    818818                k=0;
    819819                for (i=0;i<nbe;i++){
    … …  
    835835
    836836                /*Element Connectivity*/
    837                 if(verbose>5) _printf_("      writing Element connectivity" << "\n");
     837                if(verbose>5) _printf_("      writing Element connectivity\n");
    838838                bamgmesh->ElementConnectivitySize[0]=nbt-nbtout;
    839839                bamgmesh->ElementConnectivitySize[1]=3;
    … …  
    855855
    856856                /*ElementNodal Connectivity*/
    857                 if(verbose>5) _printf_("      writing Nodal element connectivity" << "\n");
     857                if(verbose>5) _printf_("      writing Nodal element connectivity\n");
    858858                bamgmesh->NodalElementConnectivitySize[0]=nbv;
    859859                bamgmesh->NodalElementConnectivitySize[1]=connectivitymax_1;
    … …  
    870870
    871871                /*Nodal Connectivity*/
    872                 if(verbose>5) _printf_("      writing Nodal connectivity" << "\n");
     872                if(verbose>5) _printf_("      writing Nodal connectivity\n");
    873873                //chaining algorithm (again...)
    874874                int* head_2=NULL;
    … …  
    924924
    925925                /*Cracked vertices*/
    926                 if(verbose>5) _printf_("      writing Cracked vertices" << "\n");
     926                if(verbose>5) _printf_("      writing Cracked vertices\n");
    927927                bamgmesh->CrackedVerticesSize[0]=NbCrackedVertices;
    928928                bamgmesh->CrackedVerticesSize[1]=2;
    … …  
    936936
    937937                /*Cracked vertices*/
    938                 if(verbose>5) _printf_("      writing Cracked vertices" << "\n");
     938                if(verbose>5) _printf_("      writing Cracked vertices\n");
    939939                bamgmesh->CrackedEdgesSize[0]=NbCrackedEdges;
    940940                bamgmesh->CrackedEdgesSize[1]=4;
    … …  
    964964                int  i,j;
    965965
    966                 if(bamgopts->verbose>3) _printf_("      processing metric" << "\n");
     966                if(bamgopts->verbose>3) _printf_("      processing metric\n");
    967967                double hmin = Max(bamgopts->hmin,MinimalHmin());
    968968                double hmax = Min(bamgopts->hmax,MaximalHmax());
    … …  
    13151315
    13161316                //display info
    1317                 if (verbose>1) _printf_("   construction of the geometry from the 2d mesh" << "\n");
     1317                if (verbose>1) _printf_("   construction of the geometry from the 2d mesh\n");
    13181318
    13191319                //check that the mesh is not empty
    … …  
    13751375                                //else (see 3 lines above), the edge has been called more than twice: return error
    13761376                                else {
    1377                                         _printf_("The edge (" << GetId(triangles[i][VerticesOfTriangularEdge[j][0]]) << "," << GetId(triangles[i][VerticesOfTriangularEdge[j][1]]) << ") belongs to more than 2 triangles (" << k << ")" << "\n");
     1377                                        _printf_("The edge (" << GetId(triangles[i][VerticesOfTriangularEdge[j][0]]) << "," << GetId(triangles[i][VerticesOfTriangularEdge[j][1]]) << ") belongs to more than 2 triangles (" << k << ")\n");
    13781378                                        _printf_("Edge " << j << " of triangle " << i << "\n");
    13791379                                        _printf_("Edge " << (-st[k]+2)%3 << " of triangle " << (-st[k]+2)/3 << "\n");
    … …  
    13901390                //display info
    13911391                if(verbose>5) {
    1392                         _printf_("         info on Mesh:" << "\n");
     1392                        _printf_("         info on Mesh:\n");
    13931393                        _printf_("            - number of vertices    = " << nbv << "\n");
    13941394                        _printf_("            - number of triangles   = " << nbt << "\n");
    … …  
    14521452                                                edges[add].ReferenceNumber=edgessave[i].ReferenceNumber;                     
    14531453                                                edges[add].GeomEdgeHook=edgessave[i].GeomEdgeHook; //  HACK to get required edges
    1454                                                 _printf_("oh no..." << "\n");
     1454                                                _printf_("oh no...\n");
    14551455                                        }
    14561456                                        else
    … …  
    17911791                //display infos
    17921792                if(verbose>1) {
    1793                         _printf_("   Construction of Metric: number of field: " << nbsol << " (nbt=" << nbt << ", nbv=" << nbv << ")" << "\n");
     1793                        _printf_("   Construction of Metric: number of field: " << nbsol << " (nbt=" << nbt << ", nbv=" << nbv << ")\n");
    17941794                }
    17951795
    … …  
    18671867                        //skip constant field
    18681868                        if (sdelta < 1.0e-10*Max(absmax,1e-20)){
    1869                                 _printf_("      Solution " << nusol << " is constant, skipping..." << "\n");
     1869                                _printf_("      Solution " << nusol << " is constant, skipping...\n");
    18701870                                continue;
    18711871                        }
    … …  
    19931993                //display infos
    19941994                if(verbose>1) {
    1995                         _printf_("   Construction of Metric: number of field: " << nbsol << " (nbt=" << nbt << ", nbv=" << nbv << ")" << "\n");
     1995                        _printf_("   Construction of Metric: number of field: " << nbsol << " (nbt=" << nbt << ", nbv=" << nbv << ")\n");
    19961996                }
    19971997
    … …  
    20802080                        //skip constant field
    20812081                        if (sdelta < 1.0e-10*Max(absmax,1e-20) ){
    2082                                 if (verbose>2) _printf_("      Solution " << nusol << " is constant, skipping..." << "\n");
     2082                                if (verbose>2) _printf_("      Solution " << nusol << " is constant, skipping...\n");
    20832083                                continue;
    20842084                        }
    … …  
    23752375                                        }
    23762376                                }
    2377                                 //_printf_(element_renu[GetId(ta.t)] << " -> " << GetId((*ta.t)[0])+1 << " " << GetId((*ta.t)[1])+1 << " " << GetId((*ta.t)[2])+1 << ", edge [" << i1 << "->" << j1 << " " << i2 << "->" << j2 << "]" << "\n");
     2377                                //_printf_(element_renu[GetId(ta.t)] << " -> " << GetId((*ta.t)[0])+1 << " " << GetId((*ta.t)[1])+1 << " " << GetId((*ta.t)[2])+1 << ", edge [" << i1 << "->" << j1 << " " << i2 << "->" << j2 << "]\n");
    23782378                                ta = Next(ta).Adj();
    23792379                                if (count++>50) _error_("Maximum number of iteration exceeded");
    … …  
    23992399                int i;
    24002400
    2401                 _printf_("Mesh Echo:" << "\n");
     2401                _printf_("Mesh Echo:\n");
    24022402                _printf_("   nbv = " << nbv << "\n");
    24032403                _printf_("   nbt = " << nbt << "\n");
    24042404                _printf_("   nbe = " << nbe << "\n");
    24052405                _printf_("   nbq = " << nbq << "\n");
    2406                 _printf_("   index:" << "\n");
     2406                _printf_("   index:\n");
    24072407                for (i=0;i<nbt;i++){
    24082408                        _printf_("   " << setw(4) << i+1 << ": ["
    … …  
    24112411                                                << setw(4) << (((BamgVertex *)triangles[i](0))?GetId(triangles[i][2])+1:0) << "]");
    24122412                }
    2413                 _printf_("   coordinates:" << "\n");
     2413                _printf_("   coordinates:\n");
    24142414                for (i=0;i<nbv;i++){
    2415                         _printf_("   " << setw(4) << i+1 << ": [" << vertices[i].r.x << " " << vertices[i].r.y << "]" << "\n");
     2415                        _printf_("   " << setw(4) << i+1 << ": [" << vertices[i].r.x << " " << vertices[i].r.y << "]\n");
    24162416                }
    24172417
    … …  
    24682468
    24692469                if (verbose >2){
    2470                         if (OutSide) _printf_("   Find all external sub-domain" << "\n");
    2471                         else _printf_("   Find all internal sub-domain" << "\n");
     2470                        if (OutSide) _printf_("   Find all external sub-domain\n");
     2471                        else _printf_("   Find all internal sub-domain\n");
    24722472                  }
    24732473                short * HeapArete = new short[nbt];
    … …  
    26922692
    26932693                                if (inew < nbsubdomains) {
    2694                                         if (verbose>5) _printf_("WARNING: " << nbsubdomains-inew << " SubDomains are being removed" << "\n");
     2694                                        if (verbose>5) _printf_("WARNING: " << nbsubdomains-inew << " SubDomains are being removed\n");
    26952695                                        nbsubdomains=inew;}
    26962696
    … …  
    27992799
    28002800                //Display info
    2801                 if (verbose>2) _printf_("   Insert initial " << nbv << " vertices" << "\n");
     2801                if (verbose>2) _printf_("   Insert initial " << nbv << " vertices\n");
    28022802
    28032803                //Compute integer coordinates for the existing vertices
    … …  
    28902890                /*Now, add the vertices One by One*/
    28912891                long NbSwap=0;
    2892                 if (verbose>3) _printf_("   Begining of insertion process..." << "\n");
     2892                if (verbose>3) _printf_("   Begining of insertion process...\n");
    28932893
    28942894                for (int icount=2; icount<nbv; icount++) {
    … …  
    29522952
    29532953                //display info if required
    2954                 if (verbose>5) _printf_("      Try to Insert " << nbvnew << " new points" << "\n");
     2954                if (verbose>5) _printf_("      Try to Insert " << nbvnew << " new points\n");
    29552955
    29562956                //return if no new points
    … …  
    30553055                        if (!e[i]){
    30563056                                kk++;
    3057                                 if(kk<10) _printf_("BUG: the geometrical edge " << i << " is on no edge curve" << "\n");
     3057                                if(kk<10) _printf_("BUG: the geometrical edge " << i << " is on no edge curve\n");
    30583058                        }
    30593059                }
    … …  
    30723072
    30733073                if (costheta >1) {
    3074                         if (verbose>5) _printf_("   do nothing: costheta > 1" << "\n");
     3074                        if (verbose>5) _printf_("   do nothing: costheta > 1\n");
    30753075                }
    30763076
    … …  
    32233223                /*First, insert old points if requested*/
    32243224                if (KeepVertices && (&Bh != this) && (nbv+Bh.nbv< maxnbv)){
    3225                         if (verbose>5) _printf_("         Inserting initial mesh points" << "\n");
     3225                        if (verbose>5) _printf_("         Inserting initial mesh points\n");
    32263226                        for (i=0;i<Bh.nbv;i++){
    32273227                                BamgVertex &bv=Bh[i];
    … …  
    32423242                // Big loop (most time consuming)
    32433243                int iter=0;
    3244                 if (verbose>5) _printf_("         Big loop" << "\n");
     3244                if (verbose>5) _printf_("         Big loop\n");
    32453245                do {
    32463246                        /*Update variables*/
    … …  
    34733473        // find extrema coordinates of vertices pmin,pmax
    34743474        long i;
    3475         if(verbose>2) _printf_("      Reconstruct mesh of " << nbv << " vertices" << "\n");
     3475        if(verbose>2) _printf_("      Reconstruct mesh of " << nbv << " vertices\n");
    34763476
    34773477        //initialize orderedvertices
    … …  
    35303530        //Display info if required
    35313531        if(verbose>5) {
    3532                 _printf_("         info of Mesh:" << "\n");
    3533                 _printf_("            - number of vertices    = " << nbv << " " << "\n");
    3534                 _printf_("            - number of triangles   = " << nbt << " " << "\n");
    3535                 _printf_("            - number of given edges = " << nbe << " " << "\n");
     3532                _printf_("         info of Mesh:\n");
     3533                _printf_("            - number of vertices    = " << nbv << " \n");
     3534                _printf_("            - number of triangles   = " << nbt << " \n");
     3535                _printf_("            - number of given edges = " << nbe << " \n");
    35363536                _printf_("            - number of all edges   = " << edge4->nb() << "\n");
    35373537                _printf_("            - Euler number 1 - nb of holes = " << nbt-edge4->nb()+nbv << "\n");
    … …  
    35553555                                }
    35563556                                else if (k==10){
    3557                                         _printf_("Other lost boundary edges not shown..." << "\n");
     3557                                        _printf_("Other lost boundary edges not shown...\n");
    35583558                                }
    35593559                        }
    … …  
    37433743                                        /*Check that the 2 vertices are on geometry AND required*/
    37443744                                        if(!edges[i][j].GeomEdgeHook->IsRequiredVertex()){
    3745                                                 _printf_("ReconstructExistingMesh error message: problem with the edge number " << i+1 << ": [" << GetId(edges[i][0])+1 << " " << GetId(edges[i][1])+1 << "]" << "\n");
     3745                                                _printf_("ReconstructExistingMesh error message: problem with the edge number " << i+1 << ": [" << GetId(edges[i][0])+1 << " " << GetId(edges[i][1])+1 << "]\n");
    37463746                                                _printf_("This edge is on geometrical edge number " << Gh.GetId(edges[i].GeomEdgeHook)+1 << "\n");
    37473747                                                if (edges[i][j].GeomEdgeHook->OnGeomVertex())
    … …  
    37523752                                                 _printf_("Its pointer is " << edges[i][j].GeomEdgeHook << "\n");
    37533753
    3754                                                 _printf_("This edge is on geometry and has no adjacent edge (open curve) and one of the tip is not required" << "\n");
     3754                                                _printf_("This edge is on geometry and has no adjacent edge (open curve) and one of the tip is not required\n");
    37553755                                                _error_("See above (might be cryptic...)");
    37563756                                        }
    … …  
    38403840                long it,ie,i;
    38413841
    3842                 _printf_("renumbering triangles" << "\n");
     3842                _printf_("renumbering triangles\n");
    38433843                for ( it=0;it<nbt;it++)
    38443844                 triangles[it].Renumbering(vertices,ve,renu);
    38453845
    3846                 _printf_("renumbering edges" << "\n");
     3846                _printf_("renumbering edges\n");
    38473847                for ( ie=0;ie<nbe;ie++)
    38483848                 edges[ie].Renumbering(vertices,ve,renu);
    38493849
    3850                 _printf_("renumbering vertices on geom" << "\n");
     3850                _printf_("renumbering vertices on geom\n");
    38513851                for (i=0;i< NbVerticesOnGeomVertex;i++)
    38523852                  {
    … …  
    38563856                  }
    38573857
    3858                 _printf_("renumbering vertices on edge" << "\n");
     3858                _printf_("renumbering vertices on edge\n");
    38593859                for (i=0;i< NbVerticesOnGeomEdge;i++)
    38603860                  {
    … …  
    38643864                  }
    38653865
    3866                 _printf_("renumbering vertices on Bth vertex" << "\n");
     3866                _printf_("renumbering vertices on Bth vertex\n");
    38673867                for (i=0;i< NbVertexOnBThVertex;i++)
    38683868                  {
    … …  
    39543954                                number_of_errors++;
    39553955                                if (number_of_errors<20){
    3956                                         _printf_("Area of Triangle " << i+1 << " < 0 (det=" << triangles[i].det << ")" << "\n");
     3956                                        _printf_("Area of Triangle " << i+1 << " < 0 (det=" << triangles[i].det << ")\n");
    39573957                                }
    39583958                        }
    … …  
    40224022        gammamn=sqrt(gammamn);
    40234023        gammamx=sqrt(gammamx);   
    4024         _printf_("   Adaptmesh info:" << "\n");
     4024        _printf_("   Adaptmesh info:\n");
    40254025        _printf_("      number of triangles = " << nt << "\n");
    40264026        _printf_("      hmin = " << hmin << ", hmax=" << hmax << "\n");
    … …  
    40624062                } 
    40634063        _printf_(" --- Histogram of the unit mesh,  nb of edges = " << nbedges << "\n");
    4064         _printf_("      length of edge in   | %% of edge  | Nb of edges " << "\n");
    4065         _printf_("      --------------------+-------------+-------------" << "\n");
     4064        _printf_("      length of edge in   | %% of edge  | Nb of edges \n");
     4065        _printf_("      --------------------+-------------+-------------\n");
    40664066        for   (i=0;i<=kmax;i++){
    40674067                if (i==0) _printf_( "      " << setw(10) << 0.);
    … …  
    40694069                if (i==kmax) _printf_("          +inf   ");
    40704070                else      _printf_( "      " << setw(10) << exp(lmin+(i+1)/delta));
    4071                 _printf_("|  " << setw(10) << (long((10000. * histo[i])/ nbedges)/100.) << " |" << "\n");
     4071                _printf_("|  " << setw(10) << (long((10000. * histo[i])/ nbedges)/100.) << " |\n");
    40724072                _printf_("  " << histo[i] << "\n");
    40734073        }
    4074         _printf_("      --------------------+-------------+-------------" << "\n");
     4074        _printf_("      --------------------+-------------+-------------\n");
    40754075}
    40764076/*}}}*/
    … …  
    47484748                        Triangle *tcvi = TriangleFindFromCoord(vi.i,det3);
    47494749                        if (tcvi && !tcvi->link) {
    4750                                 _printf_("problem inserting point in SplitInternalEdgeWithBorderVertices (tcvj && !tcvj->link)" << "\n");
     4750                                _printf_("problem inserting point in SplitInternalEdgeWithBorderVertices (tcvj && !tcvj->link)\n");
    47514751                        }
    47524752
    … …  
    47654765                }
    47664766        }
    4767         if (NbSplitEdge>nbv-nbvold) _printf_("WARNING: not enough vertices  to split all internal edges, we lost " << NbSplitEdge - ( nbv-nbvold) << " edges..." << "\n");
     4767        if (NbSplitEdge>nbv-nbvold) _printf_("WARNING: not enough vertices  to split all internal edges, we lost " << NbSplitEdge - ( nbv-nbvold) << " edges...\n");
    47684768        if (verbose>2) _printf_("SplitInternalEdgeWithBorderVertices: Number of splited edge " << NbSplitEdge << "\n");
    47694769
    … …  
    49364936
    49374937        //Vertices
    4938         if(verbose) _printf_("Reading vertices (" << nbv << ")" << "\n");
     4938        if(verbose) _printf_("Reading vertices (" << nbv << ")\n");
    49394939        for (i=0;i<nbv;i++){
    49404940                vertices[i].r.x=x[i];
    … …  
    52495249                //Insert points inside existing triangles
    52505250                if (verbose>4) _printf_("      -- current number of vertices = " << nbv << "\n");
    5251                 if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation..." << "\n");
    5252                 if (verbose>3) _printf_("         Inserting boundary points" << "\n");
     5251                if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation...\n");
     5252                if (verbose>3) _printf_("         Inserting boundary points\n");
    52535253                Insert();
    52545254
    52555255                //Force the boundary
    5256                 if (verbose>3) _printf_("         Forcing boundaries" << "\n");
     5256                if (verbose>3) _printf_("         Forcing boundaries\n");
    52575257                ForceBoundary();
    52585258
    52595259                //Extract SubDomains
    5260                 if (verbose>3) _printf_("         Extracting subdomains" << "\n");
     5260                if (verbose>3) _printf_("         Extracting subdomains\n");
    52615261                FindSubDomain();
    52625262
    5263                 if (verbose>3) _printf_("      Inserting internal points" << "\n");
     5263                if (verbose>3) _printf_("      Inserting internal points\n");
    52645264                NewPoints(*this,bamgopts,0) ;
    52655265                if (verbose>4) _printf_("      -- current number of vertices = " << nbv << "\n");
    … …  
    53675367
    53685368                                        //_printf_("Dealing with curve number " << nc << "\n");
    5369                                         //_printf_("edge on geometry is same as GhCurve? " << (ei.GeomEdgeHook==Gh.curves[nc].FirstEdge || ei.GeomEdgeHook==Gh.curves[nc].LastEdge)?"yes":"no" << "\n");
     5369                                        //_printf_("edge on geometry is same as GhCurve? " << (ei.GeomEdgeHook==Gh.curves[nc].FirstEdge || ei.GeomEdgeHook==Gh.curves[nc].LastEdge)?"yes":"no\n");
    53705370                                        //if(ei.GeomEdgeHook==Gh.curves[nc].FirstEdge || ei.GeomEdgeHook==Gh.curves[nc].LastEdge){
    5371                                         //      _printf_("Do we have the right extremity? curve first vertex -> " << ((GeomVertex *)*ei[je].GeomEdgeHook==&(*Gh.curves[nc].FirstEdge)[Gh.curves[nc].FirstVertexIndex])?"yes":"no" << "\n");
    5372                                         //      _printf_("Do we have the right extremity? curve last  vertex -> " << ((GeomVertex *)*ei[je].GeomEdgeHook==&(*Gh.curves[nc].LastEdge)[Gh.curves[nc].LastVertexIndex])?"yes":"no" << "\n");
     5371                                        //      _printf_("Do we have the right extremity? curve first vertex -> " << ((GeomVertex *)*ei[je].GeomEdgeHook==&(*Gh.curves[nc].FirstEdge)[Gh.curves[nc].FirstVertexIndex])?"yes":"no\n");
     5372                                        //      _printf_("Do we have the right extremity? curve last  vertex -> " << ((GeomVertex *)*ei[je].GeomEdgeHook==&(*Gh.curves[nc].LastEdge)[Gh.curves[nc].LastVertexIndex])?"yes":"no\n");
    53735373                                        //}
    53745374                                        //BUG FIX from original bamg
    … …  
    55785578                //Insert points inside existing triangles
    55795579                if (verbose>4) _printf_("      -- current number of vertices = " << nbv << "\n");
    5580                 if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation..." << "\n");
    5581                 if (verbose>3) _printf_("         Inserting boundary points" << "\n");
     5580                if (verbose>3) _printf_("      Creating initial Constrained Delaunay Triangulation...\n");
     5581                if (verbose>3) _printf_("         Inserting boundary points\n");
    55825582                Insert();
    55835583
    55845584                //Force the boundary
    5585                 if (verbose>3) _printf_("         Forcing boundaries" << "\n");
     5585                if (verbose>3) _printf_("         Forcing boundaries\n");
    55865586                ForceBoundary();
    55875587
    55885588                //Extract SubDomains
    5589                 if (verbose>3) _printf_("         Extracting subdomains" << "\n");
     5589                if (verbose>3) _printf_("         Extracting subdomains\n");
    55905590                FindSubDomain();
    55915591
    5592                 if (verbose>3) _printf_("      Inserting internal points" << "\n");
     5592                if (verbose>3) _printf_("      Inserting internal points\n");
    55935593                NewPoints(BTh,bamgopts,KeepVertices) ;
    55945594                if (verbose>4) _printf_("      -- current number of vertices = " << nbv << "\n");

  • issm/trunk-jpl/src/c/bamg/Metric.cpp

    r15100 r15104  
    7171        void Metric::Echo(void){
    7272
    73                 _printf_("Metric:" << "\n");
    74                 _printf_("   [a11 a21 a22]: [" << a11 << " " << a21 << " " << a22 << "]" << "\n");
     73                _printf_("Metric:\n");
     74                _printf_("   [a11 a21 a22]: [" << a11 << " " << a21 << " " << a22 << "]\n");
    7575
    7676                return;

  • issm/trunk-jpl/src/c/bamg/Triangle.cpp

    r15100 r15104  
    109109                int i;
    110110
    111                 _printf_("Triangle:" << "\n");
    112                 _printf_("   vertices pointer towards three vertices" << "\n");
     111                _printf_("Triangle:\n");
     112                _printf_("   vertices pointer towards three vertices\n");
    113113                _printf_("      vertices[0] vertices[1] vertices[2] = " << vertices[0] << " " << vertices[1] << " " << vertices[2] << "\n");
    114                 _printf_("   adj pointer towards three adjacent triangles" << "\n");
     114                _printf_("   adj pointer towards three adjacent triangles\n");
    115115                _printf_("      adj[0] adj[1] adj[2] = " << adj[0] << " " << adj[1] << " " << adj[2] << "\n");
    116116                _printf_("   det (integer triangle determinant) = " << det << "\n");
    … …  
    122122                }
    123123
    124                 _printf_("\nThree vertices:" << "\n");
     124                _printf_("\nThree vertices:\n");
    125125                for(i=0;i<3;i++){
    126126                        if (vertices[i]){
    … …  
    128128                        }
    129129                        else{
    130                                 _printf_("   vertex " << i+1 << " does not exist" << "\n");
     130                                _printf_("   vertex " << i+1 << " does not exist\n");
    131131                        }
    132132                }

  • issm/trunk-jpl/src/c/classes/Constraints/SpcDynamic.cpp

    r15100 r15104  
    4242void SpcDynamic::Echo(void){
    4343
    44         _printf_("SpcDynamic:" << "\n");
     44        _printf_("SpcDynamic:\n");
    4545        _printf_("   sid: " << sid << "\n");
    4646        _printf_("   nodeid: " << nodeid << "\n");

  • issm/trunk-jpl/src/c/classes/Constraints/SpcStatic.cpp

    r15100 r15104  
    4141void SpcStatic::Echo(void){
    4242
    43         _printf_("SpcStatic:" << "\n");
     43        _printf_("SpcStatic:\n");
    4444        _printf_("   sid: " << sid << "\n");
    4545        _printf_("   nodeid: " << nodeid << "\n");
    … …  
    5353void SpcStatic::DeepEcho(void){
    5454
    55         _printf_("SpcStatic:" << "\n");
     55        _printf_("SpcStatic:\n");
    5656        _printf_("   sid: " << sid << "\n");
    5757        _printf_("   nodeid: " << nodeid << "\n");

  • issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.cpp

    r15100 r15104  
    5656
    5757        int i;
    58         _printf_("SpcTransient:" << "\n");
     58        _printf_("SpcTransient:\n");
    5959        _printf_("   sid: " << sid << "\n");
    6060        _printf_("   nodeid: " << nodeid << "\n");
    … …  
    6262        _printf_("   nsteps: " << nsteps << "\n");
    6363        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    64         _printf_("   steps|times|values" << "\n");
     64        _printf_("   steps|times|values\n");
    6565        for(i=0;i<nsteps;i++){
    6666                _printf_(i << "-" << times[i] << ":" << values[i] << "\n");

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

    r15100 r15104  
    6161                        _printf_("  closed: " << (closed?"true":"false") << "\n");
    6262                        if(nods){
    63                                 _printf_("   x , y:" << "\n");
     63                                _printf_("   x , y:\n");
    6464                                for(int i=0;i<nods;i++){
    6565                                        _printf_(i << ": " << x[i] << " | " << y[i] << "\n");

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

    r15100 r15104  
    4040void DependentObject::Echo(void){
    4141
    42         _printf_("DependentObject:" << "\n");
     42        _printf_("DependentObject:\n");
    4343        _printf_("   name: " << EnumToStringx(this->name) << "\n");
    4444        if(this->type==0)
    45                 _printf_("   type: scalar" << "\n");
     45                _printf_("   type: scalar\n");
    4646        else if(this->type==1)
    47                 _printf_("   type: vertex" << "\n");
     47                _printf_("   type: vertex\n");
    4848        else
    4949                _error_(" unknown type: " << this->type);

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

    r15100 r15104  
    163163void DofIndexing::Echo(void){
    164164
    165         _printf_("DofIndexing:" << "\n");
     165        _printf_("DofIndexing:\n");
    166166        _printf_("   gsize: " << gsize << "\n");
    167167        _printf_("   clone: " << clone << "\n");
    … …  
    173173        int i;
    174174
    175         _printf_("DofIndexing:" << "\n");
     175        _printf_("DofIndexing:\n");
    176176        _printf_("   gsize: " << gsize << "\n");
    177177        _printf_("   fsize: " << fsize << "\n");
    … …  
    179179        _printf_("   clone: " << clone << "\n");
    180180
    181         _printf_("   set membership: f,s sets " << "\n");
     181        _printf_("   set membership: f,s sets \n");
    182182        for(i=0;i<gsize;i++){
    183183                _printf_("      dof " << i << ": " <<(f_set[i]?"true":"false")<< " " <<(s_set[i]?"true":"false") << "\n");
    … …  
    188188                if(this->s_set[i])_printf_(" " << svalues[i] << " |");
    189189        }
    190         _printf_("" << "\n");
     190        _printf_("\n");
    191191
    192192        if(doftype){
    … …  
    195195                        _printf_(" " << doftype[i] << " |");
    196196                }
    197                 _printf_("" << "\n");
    198         }
    199         else _printf_("   doftype: NULL" << "\n");
     197                _printf_("\n");
     198        }
     199        else _printf_("   doftype: NULL\n");
    200200
    201201        _printf_("   g_doflist (" << this->gsize << "): |");
    … …  
    203203                _printf_(" " << gdoflist[i] << " |");
    204204        }
    205         _printf_("" << "\n");
     205        _printf_("\n");
    206206
    207207        _printf_("   f_doflist (" << this->fsize << "): |");
    … …  
    209209                _printf_(" " << fdoflist[i] << " |");
    210210        }
    211         _printf_("" << "\n");
     211        _printf_("\n");
    212212
    213213        _printf_("   s_doflist (" << this->ssize << "): |");
    … …  
    215215                _printf_(" " << sdoflist[i] << " |");
    216216        }
    217         _printf_("" << "\n");
     217        _printf_("\n");
    218218}               
    219219/*}}}*/

  • issm/trunk-jpl/src/c/classes/ElementResults/BoolElementResult.cpp

    r15100 r15104  
    4242void BoolElementResult::DeepEcho(void){
    4343
    44         _printf_("BoolElementResult:" << "\n");
    45         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     44        _printf_("BoolElementResult:\n");
     45        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4646        _printf_("   value: "<<(this->value?"true":"false") << "\n");
    4747        _printf_("   step: " << this->step << "\n");

  • issm/trunk-jpl/src/c/classes/ElementResults/DoubleElementResult.cpp

    r15100 r15104  
    4242void DoubleElementResult::DeepEcho(void){
    4343
    44         _printf_("DoubleElementResult:" << "\n");
    45         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     44        _printf_("DoubleElementResult:\n");
     45        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4646        _printf_("   value: " << this->value << "\n");
    4747        _printf_("   step: " << this->step << "\n");

  • issm/trunk-jpl/src/c/classes/ElementResults/PentaP1ElementResult.cpp

    r15100 r15104  
    4444void PentaP1ElementResult::DeepEcho(void){
    4545
    46         _printf_("PentaP1ElementResult:" << "\n");
    47         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    48         _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << " " << this->values[3] << " " << this->values[4] << " " << this->values[5] << "]" << "\n");
     46        _printf_("PentaP1ElementResult:\n");
     47        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     48        _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << " " << this->values[3] << " " << this->values[4] << " " << this->values[5] << "]\n");
    4949        _printf_("   step: " << this->step << "\n");
    5050        _printf_("   time: " << this->time << "\n");

  • issm/trunk-jpl/src/c/classes/ElementResults/TriaP1ElementResult.cpp

    r15100 r15104  
    4444void TriaP1ElementResult::DeepEcho(void){
    4545
    46         _printf_("TriaP1ElementResult:" << "\n");
    47         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    48         _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << "]" << "\n");
     46        _printf_("TriaP1ElementResult:\n");
     47        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     48        _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << "]\n");
    4949        _printf_("   step: " << this->step << "\n");
    5050        _printf_("   time: " << this->time << "\n");

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

    r15100 r15104  
    668668void Penta::DeepEcho(void){
    669669
    670         _printf_("Penta:" << "\n");
     670        _printf_("Penta:\n");
    671671        _printf_("   id: " << id << "\n");
    672672        nodes[0]->DeepEcho();
    … …  
    679679        matpar->DeepEcho();
    680680        _printf_("   neighbor ids: " << verticalneighbors[0]->Id() << "-" << verticalneighbors[1]->Id() << "\n");
    681         _printf_("   parameters" << "\n");
     681        _printf_("   parameters\n");
    682682        parameters->DeepEcho();
    683         _printf_("   inputs" << "\n");
     683        _printf_("   inputs\n");
    684684        inputs->DeepEcho();
    685         _printf_("   results" << "\n");
     685        _printf_("   results\n");
    686686        results->DeepEcho();
    687         _printf_("neighboor sids: " << "\n");
     687        _printf_("neighboor sids: \n");
    688688        _printf_(" " << horizontalneighborsids[0] << " " << horizontalneighborsids[1] << " " << horizontalneighborsids[2] << "\n");
    689689}
    … …  
    13351335        }
    13361336        else{
    1337                 _printf_("Interpolation " << EnumToStringx(interp) << " not supported" << "\n");
     1337                _printf_("Interpolation " << EnumToStringx(interp) << " not supported\n");
    13381338        }
    13391339}
    … …  
    45604560                /*Check solution*/
    45614561                if(xIsNan<IssmDouble>(values[i])) _error_("NaN found in solution vector");
    4562                 //if(values[i]<0)      _printf_("temperature < 0°K found in solution vector" << "\n");
    4563                 //if(values[i]>275)    _printf_("temperature > 275°K found in solution vector (Paterson's rheology associated is negative)" << "\n");
     4562                //if(values[i]<0)      _printf_("temperature < 0°K found in solution vector\n");
     4563                //if(values[i]>275)    _printf_("temperature > 275°K found in solution vector (Paterson's rheology associated is negative)\n");
    45644564        }
    45654565

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

    r15100 r15104  
    935935void Tria::DeepEcho(void){
    936936
    937         _printf_("Tria:" << "\n");
     937        _printf_("Tria:\n");
    938938        _printf_("   id: " << id << "\n");
    939939        if(nodes){
    … …  
    942942                nodes[2]->DeepEcho();
    943943        }
    944         else _printf_("nodes = NULL" << "\n");
     944        else _printf_("nodes = NULL\n");
    945945
    946946        if (material) material->DeepEcho();
    947         else _printf_("material = NULL" << "\n");
     947        else _printf_("material = NULL\n");
    948948
    949949        if (matpar) matpar->DeepEcho();
    950         else _printf_("matpar = NULL" << "\n");
    951 
    952         _printf_("   parameters" << "\n");
     950        else _printf_("matpar = NULL\n");
     951
     952        _printf_("   parameters\n");
    953953        if (parameters) parameters->DeepEcho();
    954         else _printf_("parameters = NULL" << "\n");
    955 
    956         _printf_("   inputs" << "\n");
     954        else _printf_("parameters = NULL\n");
     955
     956        _printf_("   inputs\n");
    957957        if (inputs) inputs->DeepEcho();
    958         else _printf_("inputs=NULL" << "\n");
     958        else _printf_("inputs=NULL\n");
    959959
    960960        if (results) results->DeepEcho();
    961         else _printf_("results=NULL" << "\n");
    962 
    963         _printf_("neighboor sids: " << "\n");
     961        else _printf_("results=NULL\n");
     962
     963        _printf_("neighboor sids: \n");
    964964        _printf_(" " << horizontalneighborsids[0] << " " << horizontalneighborsids[1] << " " << horizontalneighborsids[2] << "\n");
    965965
    … …  
    10471047/*FUNCTION Tria::Echo{{{*/
    10481048void Tria::Echo(void){
    1049         _printf_("Tria:" << "\n");
     1049        _printf_("Tria:\n");
    10501050        _printf_("   id: " << id << "\n");
    10511051        if(nodes){
    … …  
    10541054                nodes[2]->Echo();
    10551055        }
    1056         else _printf_("nodes = NULL" << "\n");
     1056        else _printf_("nodes = NULL\n");
    10571057
    10581058        if (material) material->Echo();
    1059         else _printf_("material = NULL" << "\n");
     1059        else _printf_("material = NULL\n");
    10601060
    10611061        if (matpar) matpar->Echo();
    1062         else _printf_("matpar = NULL" << "\n");
    1063 
    1064         _printf_("   parameters" << "\n");
     1062        else _printf_("matpar = NULL\n");
     1063
     1064        _printf_("   parameters\n");
    10651065        if (parameters) parameters->Echo();
    1066         else _printf_("parameters = NULL" << "\n");
    1067 
    1068         _printf_("   inputs" << "\n");
     1066        else _printf_("parameters = NULL\n");
     1067
     1068        _printf_("   inputs\n");
    10691069        if (inputs) inputs->Echo();
    1070         else _printf_("inputs=NULL" << "\n");
     1070        else _printf_("inputs=NULL\n");
    10711071
    10721072        if (results) results->Echo();
    1073         else _printf_("results=NULL" << "\n");
    1074 
    1075         _printf_("neighboor sids: " << "\n");
     1073        else _printf_("results=NULL\n");
     1074
     1075        _printf_("neighboor sids: \n");
    10761076        _printf_(" " << horizontalneighborsids[0] << " " << horizontalneighborsids[1] << " " << horizontalneighborsids[2] << "\n");
    10771077}
    … …  
    15591559        }
    15601560        else{
    1561                 _printf_("Interpolation " << EnumToStringx(interp) << " not supported" << "\n");
     1561                _printf_("Interpolation " << EnumToStringx(interp) << " not supported\n");
    15621562        }
    15631563}

  • issm/trunk-jpl/src/c/classes/ExternalResults/GenericExternalResult.h

    r15100 r15104  
    3434                void GenericEcho(void){/*{{{*/
    3535                        _printf_("   id: " << this->id << "\n");
    36                         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     36                        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    3737                        _printf_("   step: " << this->step << "\n");
    3838                        _printf_("   time: " << this->time << "\n");
    … …  
    143143template <> inline void GenericExternalResult<bool>::DeepEcho(void){ /*{{{*/
    144144
    145         _printf_("GenericExternalResult<bool>:" << "\n");
     145        _printf_("GenericExternalResult<bool>:\n");
    146146        this->GenericEcho();
    147147        _printf_("   value: " <<(this->value?"true":"false") << "\n");
    … …  
    155155template <> inline void GenericExternalResult<int>::DeepEcho(void){ /*{{{*/
    156156
    157         _printf_("GenericExternalResult<int>:" << "\n");
     157        _printf_("GenericExternalResult<int>:\n");
    158158        this->GenericEcho();
    159159        _printf_("   value: " << this->value << "\n");
    … …  
    167167template <> inline void GenericExternalResult<double>::DeepEcho(void){ /*{{{*/
    168168
    169         _printf_("GenericExternalResult<double>:" << "\n");
     169        _printf_("GenericExternalResult<double>:\n");
    170170        this->GenericEcho();
    171171        _printf_("   value: " << this->value << "\n");
    … …  
    192192template <> inline void GenericExternalResult<char*>::DeepEcho(void){ /*{{{*/
    193193
    194         _printf_("GenericExternalResult<char*>:" << "\n");
     194        _printf_("GenericExternalResult<char*>:\n");
    195195        this->GenericEcho();
    196196        _printf_("   value: " << this->value << "\n");
    … …  
    252252template <> inline void GenericExternalResult<IssmPDouble*>::Echo(void){ /*{{{*/
    253253
    254         _printf_("GenericExternalResult<IssmPDouble*>:" << "\n");
     254        _printf_("GenericExternalResult<IssmPDouble*>:\n");
    255255        this->GenericEcho();
    256256        _printf_("   matrix size: " << this->M << "-" << this->N << "\n");
    … …  
    261261        int i,j;
    262262
    263         _printf_("GenericExternalResult<IssmPDouble*>:" << "\n");
     263        _printf_("GenericExternalResult<IssmPDouble*>:\n");
    264264        this->GenericEcho();
    265265
    … …  
    270270                        _printf_( " " << setw(11) << setprecision (5) << this->value[i*this->N+j]);
    271271                } 
    272                 _printf_(" ]" << "\n");
     272                _printf_(" ]\n");
    273273        } 
    274274

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

    r15100 r15104  
    118118        this->parameters->FindParam(&lockfilename,LockFileNameEnum);
    119119        if (waitonlock>0){
    120                 _printf0_("write lock file:" << "\n");
     120                _printf0_("write lock file:\n");
    121121                WriteLockFile(lockfilename);
    122122        }
    … …  
    137137        /*Before we delete the profiler, report statistics for this run: */
    138138        profiler->Tag(Finish);  //final tagging
    139         _printf0_("" << "\n");
     139        _printf0_("\n");
    140140        _printf0_("   "<<setw(40)<<left<<"FemModel initialization elapsed time:"<<profiler->DeltaTime(StartInit,FinishInit) << "\n");
    141141        _printf0_("   "<<setw(40)<<left<<"Core solution elapsed time:"<<profiler->DeltaTime(StartCore,FinishCore) << "\n");
    142         _printf0_("" << "\n");
     142        _printf0_("\n");
    143143        _printf0_("   Total elapsed time:"
    144144                        <<profiler->DeltaTimeModHour(Start,Finish)<<" hrs "
    … …  
    146146                        <<profiler->DeltaTimeModSec(Start,Finish)<<" sec"
    147147                        );
    148         _printf0_("" << "\n");
     148        _printf0_("\n");
    149149
    150150        /*Now delete: */
    … …  
    153153        /*Finalize PETSC for this model: */
    154154        #ifdef _HAVE_PETSC_
    155         _printf0_("closing Petsc" << "\n");
     155        _printf0_("closing Petsc\n");
    156156        PetscFinalize();
    157157        #endif
    … …  
    164164void FemModel::Echo(void){
    165165
    166         _printf_("FemModel echo: " << "\n");
     166        _printf_("FemModel echo: \n");
    167167        _printf_("   number of fem models: " << nummodels << "\n");
    168         _printf_("   analysis_type_list: " << "\n");
     168        _printf_("   analysis_type_list: \n");
    169169        for(int i=0;i<nummodels;i++)_printf_("     " << i << ": " << EnumToStringx(analysis_type_list[i]) << "\n");
    170         _printf_("   current analysis_type: " << "\n");
     170        _printf_("   current analysis_type: \n");
    171171        _printf_("     " << analysis_counter << ": " << EnumToStringx(analysis_type_list[analysis_counter]) << "\n");
    172172
    … …  
    208208        for(i=0;i<nummodels;i++){
    209209
    210                 if(VerboseMProcessor()) _printf0_("   Processing finite element model of analysis " << EnumToStringx(analysis_type_list[i]) << ":" << "\n");
     210                if(VerboseMProcessor()) _printf0_("   Processing finite element model of analysis " << EnumToStringx(analysis_type_list[i]) << ":\n");
    211211                analysis_type=analysis_type_list[i];
    212212                this->SetCurrentConfiguration(analysis_type);
    213213
    214214                if(i==0){
    215                         if(VerboseMProcessor()) _printf0_("      creating vertex degrees of freedom" << "\n");
     215                        if(VerboseMProcessor()) _printf0_("      creating vertex degrees of freedom\n");
    216216                        VerticesDofx(vertices,parameters); //only call once, we only have one set of vertices
    217217                }
    218218
    219                 if(VerboseMProcessor()) _printf0_("      resolving node constraints" << "\n");
     219                if(VerboseMProcessor()) _printf0_("      resolving node constraints\n");
    220220                SpcNodesx(nodes,constraints,parameters,analysis_type);
    221221
    222                 if(VerboseMProcessor()) _printf0_("      creating nodal degrees of freedom" << "\n");
     222                if(VerboseMProcessor()) _printf0_("      creating nodal degrees of freedom\n");
    223223                NodesDofx(nodes,parameters,analysis_type);
    224224
    225                 if(VerboseMProcessor()) _printf0_("      configuring element and loads" << "\n");
     225                if(VerboseMProcessor()) _printf0_("      configuring element and loads\n");
    226226                ConfigureObjectsx(elements, loads, nodes, vertices, materials,parameters);
    227227        }
    … …  
    247247void FemModel::OutputResults(void){
    248248
    249         _printf0_("write results to disk:" << "\n");
     249        _printf0_("write results to disk:\n");
    250250
    251251        /*Just call the OutputResultsx module: */
    … …  
    317317        void (*solutioncore)(FemModel*)=NULL; //core solution function pointer
    318318
    319         _printf0_("call computational core:" << "\n");
     319        _printf0_("call computational core:\n");
    320320
    321321        /*Retrieve solution_type from parameters: */
    … …  
    343343                solution_memory=profiler->Memory(FinishCore);
    344344
    345                 _printf0_("Solution elapsed time  : " << solution_time << "  Seconds" << "\n");
    346                 _printf0_("Solution elapsed flops : " << solution_flops << "  Flops" << "\n");
    347                 _printf0_("Solution memory used   : " << solution_memory << "  Bytes" << "\n");
     345                _printf0_("Solution elapsed time  : " << solution_time << "  Seconds\n");
     346                _printf0_("Solution elapsed flops : " << solution_flops << "  Flops\n");
     347                _printf0_("Solution memory used   : " << solution_memory << "  Bytes\n");
    348348
    349349                /*Add to results: */
    … …  
    829829
    830830        /*Display message*/
    831         if(VerboseModule()) _printf0_("   Generating matrices" << "\n");
     831        if(VerboseModule()) _printf0_("   Generating matrices\n");
    832832
    833833        /*retrive parameters: */
    … …  
    15511551                for(i=0;i<d_numresponses-1;i++)_printf_(d_responses[i] << "|");
    15521552                _printf_(d_responses[d_numresponses-1]);
    1553                 _printf_("" << "\n");
     1553                _printf_("\n");
    15541554        }
    15551555        /*}}}*/

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

    r15100 r15104  
    6666        int i;
    6767        if (num){
    68                 _printf_("   Hook: " << "\n");
     68                _printf_("   Hook: \n");
    6969                _printf_("      num=" << this->num << "\n");
    7070                _printf_("      ids: ");
    7171                for (i=0;i<this->num;i++) _printf_(this->ids[i] << " ");
    72                 _printf_("" << "\n");
     72                _printf_("\n");
    7373                _printf_("      offsets: ");
    7474                for (i=0;i<this->num;i++) _printf_(this->offsets[i] << " ");
    75                 _printf_("" << "\n");
     75                _printf_("\n");
    7676        }
    7777        else{
    78                 _printf_("   Hook: num=0 " << "\n");
     78                _printf_("   Hook: num=0 \n");
    7979        }
    8080}
    … …  
    8585        int i;
    8686        if (num){
    87                 _printf_("   Hook: " << "\n");
     87                _printf_("   Hook: \n");
    8888                _printf_("      num=" << this->num << "\n");
    8989                _printf_("      ids: ");
    9090                for (i=0;i<this->num;i++) _printf_(this->ids[i] << " ");
    91                 _printf_("" << "\n");
     91                _printf_("\n");
    9292                _printf_("      offsets: ");
    9393                for (i=0;i<this->num;i++) _printf_(this->offsets[i] << " ");
    94                 _printf_("" << "\n");
    95                 if (!objects) _printf_("      warning: object not hooked yet" << "\n");
     94                _printf_("\n");
     95                if (!objects) _printf_("      warning: object not hooked yet\n");
    9696                else{
    9797                        _printf_("      objects:\n   ");
    … …  
    9999                                _printf_("         object " << i << "\n");
    100100                                if(objects[i]) objects[i]->DeepEcho();
    101                                 else           _printf_("            no object hooked yet (not configured)" << "\n");
     101                                else           _printf_("            no object hooked yet (not configured)\n");
    102102                        }
    103103                }
    104104        }
    105105        else{
    106                 _printf_("   Hook: num=0 " << "\n");
     106                _printf_("   Hook: num=0 \n");
    107107        }
    108108}

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

    r15100 r15104  
    4141void IndependentObject::Echo(void){
    4242
    43         _printf_("IndependentObject:" << "\n");
     43        _printf_("IndependentObject:\n");
    4444        _printf_("   name: " << EnumToStringx(this->name) << "\n");
    4545        if(this->type==0)
    46                 _printf_("   type: scalar" << "\n");
     46                _printf_("   type: scalar\n");
    4747        else if(this->type==1)
    48                 _printf_("   type: vertex" << "\n");
     48                _printf_("   type: vertex\n");
    4949        else
    5050                _error_(" unknown type: " << this->type);

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

    r15100 r15104  
    4040void BoolInput::DeepEcho(void){
    4141
    42         _printf_("BoolInput:" << "\n");
    43         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     42        _printf_("BoolInput:\n");
     43        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4444        _printf_("   value: " <<(value?"true":"false") << "\n");
    4545}

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

    r15100 r15104  
    6767void ControlInput::DeepEcho(void){
    6868
    69         _printf_("ControlInput:" << "\n");
    70         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    71         _printf_("---values: " << "\n");     if (values)      values->Echo();
    72         _printf_("---savedvalues: " << "\n");if (savedvalues) savedvalues->Echo();
    73         _printf_("---minvalues: " << "\n");  if (minvalues)   minvalues->Echo();
    74         _printf_("---maxvalues: " << "\n");  if (maxvalues)   maxvalues->Echo();
    75         _printf_("---gradient: " << "\n");   if (gradient)    gradient->Echo();
     69        _printf_("ControlInput:\n");
     70        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     71        _printf_("---values: \n");     if (values)      values->Echo();
     72        _printf_("---savedvalues: \n");if (savedvalues) savedvalues->Echo();
     73        _printf_("---minvalues: \n");  if (minvalues)   minvalues->Echo();
     74        _printf_("---maxvalues: \n");  if (maxvalues)   maxvalues->Echo();
     75        _printf_("---gradient: \n");   if (gradient)    gradient->Echo();
    7676}
    7777/*}}}*/

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

    r15100 r15104  
    4242void DatasetInput::DeepEcho(void){
    4343
    44         _printf_("DatasetInput:" << "\n");
    45         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    46         _printf_("---inputs: " << "\n"); inputs->Echo();
     44        _printf_("DatasetInput:\n");
     45        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     46        _printf_("---inputs: \n"); inputs->Echo();
    4747}
    4848/*}}}*/

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

    r15100 r15104  
    4040void DoubleInput::DeepEcho(void){
    4141
    42         _printf_("DoubleInput:" << "\n");
    43         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     42        _printf_("DoubleInput:\n");
     43        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4444        _printf_("   value: " << this->value << "\n");
    4545}

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

    r15100 r15104  
    3535void IntInput::DeepEcho(void){
    3636
    37         _printf_("IntInput:" << "\n");
    38         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     37        _printf_("IntInput:\n");
     38        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    3939        _printf_("   value: " << (int)this->value << "\n");
    4040}

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

    r15100 r15104  
    5151void PentaP1Input::DeepEcho(void){
    5252
    53         _printf_("PentaP1Input:" << "\n");
    54         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    55         _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << " " << this->values[3] << " " << this->values[4] << " " << this->values[5] << "]" << "\n");
     53        _printf_("PentaP1Input:\n");
     54        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     55        _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << " " << this->values[3] << " " << this->values[4] << " " << this->values[5] << "]\n");
    5656}
    5757/*}}}*/

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

    r15100 r15104  
    6161        int i;
    6262
    63         _printf_("TransientInput:" << "\n");
    64         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     63        _printf_("TransientInput:\n");
     64        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    6565        _printf_("   numtimesteps: " << this->numtimesteps << "\n");
    66         _printf_("---inputs: " << "\n");
     66        _printf_("---inputs: \n");
    6767        for(i=0;i<this->numtimesteps;i++){
    68                 _printf_("   time: " << this->timesteps[i] << "  " << "\n");
     68                _printf_("   time: " << this->timesteps[i] << "  \n");
    6969                ((Input*)this->inputs->GetObjectByOffset(i))->Echo();
    7070        }

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

    r15100 r15104  
    5151void TriaP1Input::DeepEcho(void){
    5252
    53         _printf_("TriaP1Input:" << "\n");
    54         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
    55         _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << "]" << "\n");
     53        _printf_("TriaP1Input:\n");
     54        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
     55        _printf_("   values: [" << this->values[0] << " " << this->values[1] << " " << this->values[2] << "]\n");
    5656}
    5757/*}}}*/

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

    r15100 r15104  
    8484                for(int i=0;i<MaximumNumberOfEnums;i++){
    8585                        if(this->data[i]){
    86                                 _printf0_("Info: previous pointer of " << EnumToStringx(i) << " has not been freed (DeleteData has not been called)" << "\n");
     86                                _printf0_("Info: previous pointer of " << EnumToStringx(i) << " has not been freed (DeleteData has not been called)\n");
    8787                        }
    8888                }
    … …  
    125125                else{
    126126                        if(record_enum!=MaximumNumberOfEnums){
    127                                 _printf0_("" << "\n");
    128                                 _printf0_("=========================================================================" << "\n");
    129                                 _printf0_(" Enums in marshalled file are not compatible with compiled code          " << "\n");
    130                                 _printf0_("                                                                         " << "\n");
    131                                 _printf0_("   * If you are running ISSM on a remote cluster:                        " << "\n");
    132                                 _printf0_("     make sure that you are using the same version of ISSM on your local " << "\n");
    133                                 _printf0_("     machine and remote cluster (you might need to run svn update)       " << "\n");
    134                                 _printf0_("   * If you are running ISSM on your local machine:                      " << "\n");
    135                                 _printf0_("     make sure that all the code is compiled (modules and executables)   " << "\n");
    136                                 _printf0_("   * If you are a developer and just added a new Enum:                   " << "\n");
    137                                 _printf0_("     you might need to run ./Synchronize.sh in src/c/EnumDefinitions     " << "\n");
    138                                 _printf0_("     and recompile                                                       " << "\n");
    139                                 _printf0_("=========================================================================" << "\n");
    140                                 _printf0_("" << "\n");
     127                                _printf0_("\n");
     128                                _printf0_("=========================================================================\n");
     129                                _printf0_(" Enums in marshalled file are not compatible with compiled code          \n");
     130                                _printf0_("                                                                         \n");
     131                                _printf0_("   * If you are running ISSM on a remote cluster:                        \n");
     132                                _printf0_("     make sure that you are using the same version of ISSM on your local \n");
     133                                _printf0_("     machine and remote cluster (you might need to run svn update)       \n");
     134                                _printf0_("   * If you are running ISSM on your local machine:                      \n");
     135                                _printf0_("     make sure that all the code is compiled (modules and executables)   \n");
     136                                _printf0_("   * If you are a developer and just added a new Enum:                   \n");
     137                                _printf0_("     you might need to run ./Synchronize.sh in src/c/EnumDefinitions     \n");
     138                                _printf0_("     and recompile                                                       \n");
     139                                _printf0_("=========================================================================\n");
     140                                _printf0_("\n");
    141141                                _error_("Enums not consistent (See error message above)");
    142142                        }

  • issm/trunk-jpl/src/c/classes/Loads/Friction.cpp

    r15100 r15104  
    4343/*FUNCTION Friction::Echo {{{*/
    4444void Friction::Echo(void){
    45         _printf_("Friction:" << "\n");
     45        _printf_("Friction:\n");
    4646        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    4747        _printf_("   element_type: " << this->element_type << "\n");

  • issm/trunk-jpl/src/c/classes/Loads/Icefront.cpp

    r15100 r15104  
    134134/*FUNCTION Icefront::Echo {{{*/
    135135void Icefront::Echo(void){
    136         _printf_("Icefront:" << "\n");
     136        _printf_("Icefront:\n");
    137137        _printf_("   id: " << id << "\n");
    138138        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    150150void Icefront::DeepEcho(void){
    151151
    152         _printf_("Icefront:" << "\n");
     152        _printf_("Icefront:\n");
    153153        _printf_("   id: " << id << "\n");
    154154        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");

  • issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp

    r15100 r15104  
    155155/*FUNCTION Numericalflux::Echo {{{*/
    156156void Numericalflux::Echo(void){
    157         _printf_("Numericalflux:" << "\n");
     157        _printf_("Numericalflux:\n");
    158158        _printf_("   id: " << id << "\n");
    159159        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    168168void Numericalflux::DeepEcho(void){
    169169
    170         _printf_("Numericalflux:" << "\n");
     170        _printf_("Numericalflux:\n");
    171171        _printf_("   id: " << id << "\n");
    172172        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    174174        hvertices->DeepEcho();
    175175        helement->DeepEcho();
    176         _printf_("   parameters" << "\n");
     176        _printf_("   parameters\n");
    177177        if(parameters)
    178178         parameters->DeepEcho();
    179179        else
    180          _printf_("      NULL" << "\n");
    181         _printf_("   inputs" << "\n");
     180         _printf_("      NULL\n");
     181        _printf_("   inputs\n");
    182182        inputs->DeepEcho();
    183183

  • issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp

    r15100 r15104  
    103103void Pengrid::DeepEcho(void){
    104104
    105         _printf_("Pengrid:" << "\n");
     105        _printf_("Pengrid:\n");
    106106        _printf_("   id: " << id << "\n");
    107107        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    111111        _printf_("   active " << this->active << "\n");
    112112        _printf_("   zigzag_counter " << this->zigzag_counter << "\n");
    113         _printf_("   parameters" << "\n");
     113        _printf_("   parameters\n");
    114114        parameters->DeepEcho();
    115         _printf_("   inputs" << "\n");
     115        _printf_("   inputs\n");
    116116        inputs->DeepEcho();
    117117}

  • issm/trunk-jpl/src/c/classes/Loads/Penpair.cpp

    r15100 r15104  
    5151void Penpair::Echo(void){
    5252
    53         _printf_("Penpair:" << "\n");
     53        _printf_("Penpair:\n");
    5454        _printf_("   id: " << id << "\n");
    5555        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    6262void Penpair::DeepEcho(void){
    6363
    64         _printf_("Penpair:" << "\n");
     64        _printf_("Penpair:\n");
    6565        _printf_("   id: " << id << "\n");
    6666        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");

  • issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp

    r15100 r15104  
    138138        input=(Input*)this->inputs->GetInput(FractionIncrementEnum); input->GetInputValue(&fractionincrement);
    139139
    140         _printf_("Riftfront:" << "\n");
     140        _printf_("Riftfront:\n");
    141141        _printf_("   id: " << id << "\n");
    142142        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    146146        _printf_("   parameters: " << parameters << "\n");
    147147        _printf_("   inputs: " << inputs << "\n");
    148         _printf_("   internal parameters: " << "\n");
     148        _printf_("   internal parameters: \n");
    149149        _printf_("   normal: " << normal[0] << "|" << normal[1] << "\n");
    150150        _printf_("   length: " << length << "\n");
    … …  
    166166void Riftfront::DeepEcho(void){
    167167
    168         _printf_("Riftfront:" << "\n");
     168        _printf_("Riftfront:\n");
    169169        _printf_("   id: " << id << "\n");
    170170        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    … …  
    172172        helements->DeepEcho();
    173173        hmatpar->DeepEcho();
    174         _printf_("   parameters" << "\n");
     174        _printf_("   parameters\n");
    175175        if(parameters)parameters->DeepEcho();
    176         _printf_("   inputs" << "\n");
     176        _printf_("   inputs\n");
    177177        if(inputs)inputs->DeepEcho();
    178178}
    … …  
    702702        this->active=activate;
    703703
    704         //if ((penetration>0) && (this->active==1))_printf_("Riftfront " << Id() << " wants to be released" << "\n");
     704        //if ((penetration>0) && (this->active==1))_printf_("Riftfront " << Id() << " wants to be released\n");
    705705
    706706        /*assign output pointer: */

  • issm/trunk-jpl/src/c/classes/Materials/Matdamageice.cpp

    r15100 r15104  
    5555void Matdamageice::Echo(void){
    5656
    57         _printf_("Matdamageice:" << "\n");
     57        _printf_("Matdamageice:\n");
    5858        _printf_("   mid: " << mid << "\n");
    59         _printf_("   inputs:" << "\n");
     59        _printf_("   inputs:\n");
    6060        inputs->Echo();
    61         _printf_("   element:" << "\n");
     61        _printf_("   element:\n");
    6262        helement->Echo();
    6363}
    … …  
    6666void Matdamageice::DeepEcho(void){
    6767
    68         _printf_("Matdamageice:" << "\n");
     68        _printf_("Matdamageice:\n");
    6969        _printf_("   mid: " << mid << "\n");
    70         _printf_("   inputs:" << "\n");
     70        _printf_("   inputs:\n");
    7171        inputs->DeepEcho();
    72         _printf_("   element:" << "\n");
     72        _printf_("   element:\n");
    7373        helement->Echo();
    7474}               

  • issm/trunk-jpl/src/c/classes/Materials/Matice.cpp

    r15100 r15104  
    6969void Matice::Echo(void){
    7070
    71         _printf_("Matice:" << "\n");
     71        _printf_("Matice:\n");
    7272        _printf_("   mid: " << mid << "\n");
    73         _printf_("   inputs:" << "\n");
     73        _printf_("   inputs:\n");
    7474        inputs->Echo();
    75         _printf_("   element:" << "\n");
     75        _printf_("   element:\n");
    7676        helement->Echo();
    7777}
    … …  
    8080void Matice::DeepEcho(void){
    8181
    82         _printf_("Matice:" << "\n");
     82        _printf_("Matice:\n");
    8383        _printf_("   mid: " << mid << "\n");
    84         _printf_("   inputs:" << "\n");
     84        _printf_("   inputs:\n");
    8585        inputs->DeepEcho();
    86         _printf_("   element:" << "\n");
     86        _printf_("   element:\n");
    8787        helement->Echo();
    8888}               

  • issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp

    r15100 r15104  
    9595void Matpar::Echo(void){
    9696
    97         _printf_("Matpar:" << "\n");
     97        _printf_("Matpar:\n");
    9898        _printf_("   mid: " << mid << "\n");
    9999        _printf_("   rho_ice: " << rho_ice << "\n");

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

    r15100 r15104  
    155155void Node::Echo(void){
    156156
    157         _printf_("Node:" << "\n");
     157        _printf_("Node:\n");
    158158        _printf_("   id: " << id << "\n");
    159159        _printf_("   sid: " << sid << "\n");
    160160        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    161161        indexing.Echo();
    162         _printf_("   hvertex:     not displayed" << "\n");
     162        _printf_("   hvertex:     not displayed\n");
    163163        _printf_("   inputs:      " << inputs << "\n");
    164164
    … …  
    168168void Node::DeepEcho(void){
    169169
    170         _printf_("Node:" << "\n");
     170        _printf_("Node:\n");
    171171        _printf_("   id: " << id << "\n");
    172172        _printf_("   sid: " << sid << "\n");
    173173        _printf_("   analysis_type: " << EnumToStringx(analysis_type) << "\n");
    174174        indexing.DeepEcho();
    175         _printf_("Vertex:" << "\n");
     175        _printf_("Vertex:\n");
    176176        hvertex->DeepEcho();
    177         _printf_("   inputs" << "\n");
     177        _printf_("   inputs\n");
    178178
    179179}

  • issm/trunk-jpl/src/c/classes/Options/GenericOption.h

    r15100 r15104  
    6464                        bool  flag=true;
    6565
    66                         if(flag) _printf0_(indent << "         name: \"" << name << "\"" << "\n");
     66                        if(flag) _printf0_(indent << "         name: \"" << name << "\"\n");
    6767                        if(flag) _printf0_(indent << "         numel: " << numel << "\n");
    6868                        if(flag) _printf0_(indent << "         ndims: " << ndims << "\n");
    … …  
    7171                                if(flag) _printf0_(indent << "          size: " << cstr << "\n");
    7272                        }
    73                         else if(flag) _printf0_(indent << "          size: [empty]" << "\n");
     73                        else if(flag) _printf0_(indent << "          size: [empty]\n");
    7474                        _printf_(indent << "         value: " << value << "\n");;
    7575                } /*}}}*/

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

    r15100 r15104  
    6464                        }
    6565                        else{
    66                                 if(GetOption(name)) _printf_("WARNING: option "<<name<<" found but fetched format not consistent, defaulting..." << "\n");
     66                                if(GetOption(name)) _printf_("WARNING: option "<<name<<" found but fetched format not consistent, defaulting...\n");
    6767                                *pvalue=default_value;
    6868                        }

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

    r15100 r15104  
    4343void BoolParam::DeepEcho(void){
    4444
    45         _printf_("BoolParam:" << "\n");
    46         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     45        _printf_("BoolParam:\n");
     46        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4747        _printf_("   value: " <<(this->value?"true":"false") << "\n");
    4848}

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

    r15100 r15104  
    4444void DataSetParam::DeepEcho(void){
    4545
    46         _printf_("DataSetParam:" << "\n");
    47         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     46        _printf_("DataSetParam:\n");
     47        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4848        this->value->Echo();
    4949}

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

    r15100 r15104  
    8282void DoubleMatArrayParam::Echo(void){
    8383
    84         _printf_("DoubleMatArrayParam:" << "\n");
    85         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     84        _printf_("DoubleMatArrayParam:\n");
     85        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    8686        _printf_("   array size: " << this->M << "\n");
    8787        _printf_("   array pointer: " << this->array << "\n");
    … …  
    9696        IssmDouble* matrix=NULL;
    9797
    98         _printf_("DoubleMatArrayParam:" << "\n");
    99         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     98        _printf_("DoubleMatArrayParam:\n");
     99        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    100100        _printf_("   array size: " << this->M << "\n");
    101101        for(i=0;i<M;i++){
    102                 _printf_("   array " << i << " (" << mdim_array[i] << "x" << ndim_array[i] << "):" << "\n");
     102                _printf_("   array " << i << " (" << mdim_array[i] << "x" << ndim_array[i] << "):\n");
    103103                matrix=array[i];
    104104                m=mdim_array[i];
    … …  
    108108                        _printf_("   ");
    109109                        for(k=0;k<n;k++)_printf_(*(matrix+n*j+k) << " ");
    110                         _printf_("" << "\n");
     110                        _printf_("\n");
    111111                }
    112112        }

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

    r15100 r15104  
    4343void DoubleMatParam::Echo(void){
    4444
    45         _printf_("DoubleMatParam:" << "\n");
    46         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     45        _printf_("DoubleMatParam:\n");
     46        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4747        _printf_("   matrix size: " << this->M << "x" << this->N << "\n");
    4848
    … …  
    5454        int i,j;
    5555
    56         _printf_("DoubleMatParam:" << "\n");
    57         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     56        _printf_("DoubleMatParam:\n");
     57        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5858        _printf_("   matrix size: " << this->M << "x" << this->N << "\n");
    5959        for(i=0;i<this->M;i++){

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

    r15100 r15104  
    4040void DoubleParam::DeepEcho(void){
    4141
    42         _printf_("DoubleParam:" << "\n");
    43         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     42        _printf_("DoubleParam:\n");
     43        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4444        _printf_("   value: " << this->value << "\n");
    4545}

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

    r15100 r15104  
    4242void DoubleVecParam::Echo(void){
    4343
    44         _printf_("DoubleVecParam:" << "\n");
    45         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     44        _printf_("DoubleVecParam:\n");
     45        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4646        _printf_("   vector size: " << this->M << "\n");
    4747
    … …  
    5353        int i;
    5454
    55         _printf_("DoubleVecParam:" << "\n");
    56         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     55        _printf_("DoubleVecParam:\n");
     56        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5757        _printf_("   vector size: " << this->M << "\n");
    5858        for(i=0;i<this->M;i++){

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

    r15100 r15104  
    4343void FileParam::DeepEcho(void){
    4444
    45         _printf_("FileParam:" << "\n");
    46         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     45        _printf_("FileParam:\n");
     46        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4747        _printf_("   value: " << this->value << "\n");
    4848}

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

    r15100 r15104  
    3939                // that any structured P must provide the friend << operator
    4040                void  DeepEcho() {
    41                   _printf_("GenericParam:" << "\n");
    42                   _printf_("   enum:  " << myEnumVal << " (" << EnumToStringx(myEnumVal) << ")" << "\n");
     41                  _printf_("GenericParam:\n");
     42                  _printf_("   enum:  " << myEnumVal << " (" << EnumToStringx(myEnumVal) << ")\n");
    4343                  _printf_("   value: " << myP << "\n");;
    4444                }

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

    r15100 r15104  
    4343void IntMatParam::Echo(void){
    4444
    45         _printf_("IntMatParam:" << "\n");
    46         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     45        _printf_("IntMatParam:\n");
     46        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4747        _printf_("   matrix size: " << this->M << "x" << this->N << "\n");
    4848
    … …  
    5454        int i,j;
    5555
    56         _printf_("IntMatParam:" << "\n");
    57         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     56        _printf_("IntMatParam:\n");
     57        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5858        _printf_("   matrix size: " << this->M << "x" << this->N << "\n");
    5959        for(i=0;i<this->M;i++){

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

    r15100 r15104  
    4343void IntParam::DeepEcho(void){
    4444
    45         _printf_("IntParam:" << "\n");
    46         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     45        _printf_("IntParam:\n");
     46        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4747        _printf_("   value: " << this->value << "\n");
    4848}

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

    r15100 r15104  
    5858void IntVecParam::Echo(void){
    5959
    60         _printf_("IntVecParam:" << "\n");
    61         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     60        _printf_("IntVecParam:\n");
     61        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    6262        _printf_("   vector size: " << this->M << "\n");
    6363
    … …  
    6969        int i;
    7070
    71         _printf_("IntVecParam:" << "\n");
    72         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     71        _printf_("IntVecParam:\n");
     72        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    7373        _printf_("   vector size: " << this->M << "\n");
    7474        for(i=0;i<this->M;i++){

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

    r15100 r15104  
    4242void MatrixParam::Echo(void){
    4343
    44         _printf_("MatrixParam:" << "\n");
    45         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     44        _printf_("MatrixParam:\n");
     45        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4646
    4747}
    … …  
    5050void MatrixParam::DeepEcho(void){
    5151
    52         _printf_("MatrixParam:" << "\n");
    53         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     52        _printf_("MatrixParam:\n");
     53        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5454        this->value->Echo();
    5555}

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

    r15100 r15104  
    7070        char* string=NULL;
    7171
    72         _printf_("StringArrayParam:" << "\n");
    73         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     72        _printf_("StringArrayParam:\n");
     73        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    7474        for(i=0;i<this->numstrings;i++){
    7575                string=this->value[i];

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

    r15100 r15104  
    4444/*FUNCTION StringParam::DeepEcho{{{*/
    4545void StringParam::DeepEcho(void){
    46         _printf_("StringParam:" << "\n");
    47         _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     46        _printf_("StringParam:\n");
     47        _printf_("   enum:  " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4848        _printf_("   value: " << this->value << "\n");
    4949}

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

    r15100 r15104  
    4747void TransientParam::Echo(void){
    4848
    49         _printf_("TransientParam:" << "\n");
    50         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     49        _printf_("TransientParam:\n");
     50        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5151        _printf_("   size: " << this->N << "\n");
    5252
    … …  
    5656void TransientParam::DeepEcho(void){
    5757
    58         _printf_("TransientParam:" << "\n");
    59         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     58        _printf_("TransientParam:\n");
     59        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    6060        _printf_("   size: " << this->N << "\n");
    6161        for(int i=0;i<this->N;i++){

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

    r15100 r15104  
    4444void VectorParam::Echo(void){
    4545
    46         _printf_("VectorParam:" << "\n");
    47         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     46        _printf_("VectorParam:\n");
     47        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    4848
    4949}
    … …  
    5252void VectorParam::DeepEcho(void){
    5353
    54         _printf_("VectorParam:" << "\n");
    55         _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")" << "\n");
     54        _printf_("VectorParam:\n");
     55        _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
    5656        value->Echo();
    5757}

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

    r15100 r15104  
    3636void Profiler::Echo(void){
    3737
    38         _printf_("Profiler:" << "\n");
    39         _printf_("   time tags: " << "\n");
     38        _printf_("Profiler:\n");
     39        _printf_("   time tags: \n");
    4040        this->time->Echo();
    4141
    … …  
    4545void Profiler::DeepEcho(void){
    4646
    47         _printf_("Profiler:" << "\n");
    48         _printf_("   time tags: " << "\n");
     47        _printf_("Profiler:\n");
     48        _printf_("   time tags: \n");
    4949        this->time->DeepEcho();
    5050

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

    r15100 r15104  
    5252                void Echo(void){
    5353
    54                         _printf_("Segment:" << "\n");
     54                        _printf_("Segment:\n");
    5555                        _printf_("   eid: " << eid << "\n");
    5656                        _printf_("   node 1: " << this->x1 << "|" << this->y1 << "\n");

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

    r15100 r15104  
    6262void Vertex::Echo(void){
    6363
    64         _printf_("Vertex:" << "\n");
     64        _printf_("Vertex:\n");
    6565        _printf_("   id: " << id << "\n");
    6666        _printf_("   sid: " << sid << "\n");

  • issm/trunk-jpl/src/c/classes/gauss/GaussPenta.cpp

    r15100 r15104  
    251251void GaussPenta::Echo(void){
    252252
    253         _printf_("GaussPenta:" << "\n");
     253        _printf_("GaussPenta:\n");
    254254        _printf_("   numgauss: " << numgauss << "\n");
    255255
    … …  
    257257         _printf_("   weights = [");
    258258         for(int i=0;i<numgauss;i++) _printf_(" " << weights[i] << "\n");
    259          _printf_("]" << "\n");
    260         }
    261         else _printf_("weights = NULL" << "\n");
     259         _printf_("]\n");
     260        }
     261        else _printf_("weights = NULL\n");
    262262        if (coords1){
    263263         _printf_("   coords1 = [");
    264264         for(int i=0;i<numgauss;i++) _printf_(" " << coords1[i] << "\n");
    265          _printf_("]" << "\n");
    266         }
    267         else _printf_("coords1 = NULL" << "\n");
     265         _printf_("]\n");
     266        }
     267        else _printf_("coords1 = NULL\n");
    268268        if (coords2){
    269269         _printf_("   coords2 = [");
    270270         for(int i=0;i<numgauss;i++) _printf_(" " << coords2[i] << "\n");
    271          _printf_("]" << "\n");
    272         }
    273         else _printf_("coords2 = NULL" << "\n");
     271         _printf_("]\n");
     272        }
     273        else _printf_("coords2 = NULL\n");
    274274        if (coords3){
    275275         _printf_("   coords3 = [");
    276276         for(int i=0;i<numgauss;i++) _printf_(" " << coords3[i] << "\n");
    277          _printf_("]" << "\n");
    278         }
    279         else _printf_("coords3 = NULL" << "\n");
     277         _printf_("]\n");
     278        }
     279        else _printf_("coords3 = NULL\n");
    280280        if (coords4){
    281281                _printf_("   coords4 = [");
    282282                for(int i=0;i<numgauss;i++) _printf_(" " << coords4[i] << "\n");
    283                 _printf_("]" << "\n");
    284         }
    285         else _printf_("coords4 = NULL" << "\n");
     283                _printf_("]\n");
     284        }
     285        else _printf_("coords4 = NULL\n");
    286286
    287287        _printf_("   weight = " << weight << "\n");

  • issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp

    r15100 r15104  
    263263void GaussTria::Echo(void){
    264264
    265         _printf_("GaussTria:" << "\n");
     265        _printf_("GaussTria:\n");
    266266        _printf_("   numgauss: " << numgauss << "\n");
    267267
    … …  
    269269         _printf_("   weights = [");
    270270         for(int i=0;i<numgauss;i++) _printf_(" " << weights[i] << "\n");
    271          _printf_("]" << "\n");
    272         }
    273         else _printf_("weights = NULL" << "\n");
     271         _printf_("]\n");
     272        }
     273        else _printf_("weights = NULL\n");
    274274        if (coords1){
    275275         _printf_("   coords1 = [");
    276276         for(int i=0;i<numgauss;i++) _printf_(" " << coords1[i] << "\n");
    277          _printf_("]" << "\n");
    278         }
    279         else _printf_("coords1 = NULL" << "\n");
     277         _printf_("]\n");
     278        }
     279        else _printf_("coords1 = NULL\n");
    280280        if (coords2){
    281281         _printf_("   coords2 = [");
    282282         for(int i=0;i<numgauss;i++) _printf_(" " << coords2[i] << "\n");
    283          _printf_("]" << "\n");
    284         }
    285         else _printf_("coords2 = NULL" << "\n");
     283         _printf_("]\n");
     284        }
     285        else _printf_("coords2 = NULL\n");
    286286        if (coords3){
    287287         _printf_("   coords3 = [");
    288288         for(int i=0;i<numgauss;i++) _printf_(" " << coords3[i] << "\n");
    289          _printf_("]" << "\n");
    290         }
    291         else _printf_("coords3 = NULL" << "\n");
     289         _printf_("]\n");
     290        }
     291        else _printf_("coords3 = NULL\n");
    292292
    293293        _printf_("   weight = " << weight << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/ExponentialVariogram.cpp

    r15100 r15104  
    4747/*FUNCTION ExponentialVariogram::Echo {{{*/
    4848void ExponentialVariogram::Echo(void){
    49         _printf_("ExponentialVariogram" << "\n");
     49        _printf_("ExponentialVariogram\n");
    5050        _printf_("   nugget: " << this->nugget << "\n");
    5151        _printf_("   sill  : " << this->sill << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/GaussianVariogram.cpp

    r15100 r15104  
    4747/*FUNCTION GaussianVariogram::Echo {{{*/
    4848void GaussianVariogram::Echo(void){
    49         _printf_("GaussianVariogram" << "\n");
     49        _printf_("GaussianVariogram\n");
    5050        _printf_("   nugget: " << this->nugget << "\n");
    5151        _printf_("   sill  : " << this->sill << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/Observation.cpp

    r15100 r15104  
    3535void Observation::Echo(void){
    3636
    37         _printf_("Observation" << "\n");
     37        _printf_("Observation\n");
    3838        _printf_("   index : " << this->index << "\n");
    3939        _printf_("   x     : " << this->x << "\n");
    4040        _printf_("   y     : " << this->y << "\n");
    41         _printf_("   xi    : " << "\n"); printbinary(this->xi); _printf_("\n");
    42         _printf_("   yi    : " << "\n"); printbinary(this->yi); _printf_("\n");
     41        _printf_("   xi    : \n"); printbinary(this->xi); _printf_("\n");
     42        _printf_("   yi    : \n"); printbinary(this->yi); _printf_("\n");
    4343        _printf_("   weight: " << this->weight << "\n");
    4444        _printf_("   value : " << this->value << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/Observations.cpp

    r15100 r15104  
    107107                }
    108108        }
    109         _printf0_("done" << "\n");
     109        _printf0_("done\n");
    110110        _printf0_("Initial number of observations: " << n << "\n");
    111111        _printf0_("  Final number of observations: " << this->quadtree->NbObs << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/PowerVariogram.cpp

    r15100 r15104  
    4848/*FUNCTION PowerVariogram::Echo {{{*/
    4949void PowerVariogram::Echo(void){
    50         _printf_("PowerVariogram" << "\n");
     50        _printf_("PowerVariogram\n");
    5151        _printf_("   nugget: " << this->nugget << "\n");
    5252        _printf_("   slope : " << this->slope << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/Quadtree.cpp

    r15100 r15104  
    297297void  Quadtree::Echo(void){
    298298
    299         _printf_("Quadtree:" << "\n");
     299        _printf_("Quadtree:\n");
    300300        _printf_("   MaxDepth      = " << this->MaxDepth << "\n");
    301301        _printf_("   NbQuadtreeBox = " << this->NbQuadtreeBox << "\n");
    … …  
    307307void  Quadtree::DeepEcho(void){
    308308
    309         _printf_("Quadtree:" << "\n");
     309        _printf_("Quadtree:\n");
    310310        _printf_("   MaxDepth      = " << this->MaxDepth << "\n");
    311311        _printf_("   NbQuadtreeBox = " << this->NbQuadtreeBox << "\n");
    … …  
    512512void  Quadtree::QuadtreeBox::Echo(void){
    513513
    514         _printf_("QuadtreeBox:" << "\n");
     514        _printf_("QuadtreeBox:\n");
    515515        _printf_("   nbitems = " << this->nbitems << "\n");
    516516        _printf_("   xcenter = " << this->xcenter << "\n");

  • issm/trunk-jpl/src/c/classes/kriging/SphericalVariogram.cpp

    r15100 r15104  
    4747/*FUNCTION SphericalVariogram::Echo {{{*/
    4848void SphericalVariogram::Echo(void){
    49         _printf_("SphericalVariogram" << "\n");
     49        _printf_("SphericalVariogram\n");
    5050        _printf_("   nugget: " << this->nugget << "\n");
    5151        _printf_("   sill  : " << this->sill << "\n");

  • issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp

    r15100 r15104  
    377377
    378378        int i,j;
    379         _printf_("Element Matrix echo:" << "\n");
     379        _printf_("Element Matrix echo:\n");
    380380        _printf_("   nrows: " << this->nrows << "\n");
    381381        _printf_("   ncols: " << this->ncols << "\n");
    382382        _printf_("   dofsymmetrical: " << (dofsymmetrical?"true":"false") << "\n");
    383383
    384         _printf_("   values:" << "\n");
     384        _printf_("   values:\n");
    385385        for(i=0;i<nrows;i++){
    386386                _printf_(setw(4) << right << i << ": ");
    387387                for(j=0;j<ncols;j++) _printf_( " " << setw(11) << setprecision (5) << right << values[i*ncols+j]);
    388                 _printf_("" << "\n");
     388                _printf_("\n");
    389389        }
    390390
    391391        _printf_("   gglobaldoflist (" << gglobaldoflist << "): ");
    392         if(gglobaldoflist) for(i=0;i<nrows;i++) _printf_(" " << gglobaldoflist[i] << "\n"); _printf_("" << "\n");
     392        if(gglobaldoflist) for(i=0;i<nrows;i++) _printf_(" " << gglobaldoflist[i] << "\n"); _printf_("\n");
    393393
    394394        _printf_("   row_fsize: " << row_fsize << "\n");
    395395        _printf_("   row_flocaldoflist  (" << row_flocaldoflist << "): ");
    396         if(row_flocaldoflist) for(i=0;i<row_fsize;i++) _printf_(" " << row_flocaldoflist[i] << "\n"); _printf_(" " << "\n");
     396        if(row_flocaldoflist) for(i=0;i<row_fsize;i++) _printf_(" " << row_flocaldoflist[i] << "\n"); _printf_(" \n");
    397397        _printf_("   row_fglobaldoflist  (" << row_fglobaldoflist << "): ");
    398         if(row_fglobaldoflist) for(i=0;i<row_fsize;i++) _printf_(" " << row_fglobaldoflist[i] << "\n"); _printf_(" " << "\n");
     398        if(row_fglobaldoflist) for(i=0;i<row_fsize;i++) _printf_(" " << row_fglobaldoflist[i] << "\n"); _printf_(" \n");
    399399
    400400        _printf_("   row_ssize: " << row_ssize << "\n");
    401401        _printf_("   row_slocaldoflist  (" << row_slocaldoflist << "): ");
    402         if(row_slocaldoflist) for(i=0;i<row_ssize;i++) _printf_(" " << row_slocaldoflist[i] << "\n"); _printf_(" " << "\n");
     402        if(row_slocaldoflist) for(i=0;i<row_ssize;i++) _printf_(" " << row_slocaldoflist[i] << "\n"); _printf_(" \n");
    403403        _printf_("   row_sglobaldoflist  (" << row_sglobaldoflist << "): ");
    404         if(row_sglobaldoflist) for(i=0;i<row_ssize;i++) _printf_(" " << row_sglobaldoflist[i] << "\n"); _printf_(" " << "\n");
     404        if(row_sglobaldoflist) for(i=0;i<row_ssize;i++) _printf_(" " << row_sglobaldoflist[i] << "\n"); _printf_(" \n");
    405405
    406406        if(!dofsymmetrical){
    407407                _printf_("   col_fsize: " << col_fsize << "\n");
    408408                _printf_("   col_flocaldoflist  (" << col_flocaldoflist << "): ");
    409                 if(col_flocaldoflist) for(i=0;i<col_fsize;i++) _printf_(" " << col_flocaldoflist[i] << "\n"); _printf_(" " << "\n");
     409                if(col_flocaldoflist) for(i=0;i<col_fsize;i++) _printf_(" " << col_flocaldoflist[i] << "\n"); _printf_(" \n");
    410410                _printf_("   col_fglobaldoflist  (" << col_fglobaldoflist << "): ");
    411                 if(col_fglobaldoflist) for(i=0;i<col_fsize;i++) _printf_(" " << col_fglobaldoflist[i] << "\n"); _printf_(" " << "\n");
     411                if(col_fglobaldoflist) for(i=0;i<col_fsize;i++) _printf_(" " << col_fglobaldoflist[i] << "\n"); _printf_(" \n");
    412412
    413413                _printf_("   col_ssize: " << col_ssize << "\n");
    414414                _printf_("   col_slocaldoflist  (" << col_slocaldoflist << "): ");
    415                 if(col_slocaldoflist) for(i=0;i<col_ssize;i++) _printf_(" " << col_slocaldoflist[i] << "\n"); _printf_(" " << "\n");
     415                if(col_slocaldoflist) for(i=0;i<col_ssize;i++) _printf_(" " << col_slocaldoflist[i] << "\n"); _printf_(" \n");
    416416                _printf_("   col_sglobaldoflist  (" << col_sglobaldoflist << "): ");
    417                 if(col_sglobaldoflist) for(i=0;i<col_ssize;i++) _printf_(" " << col_sglobaldoflist[i] << "\n"); _printf_(" " << "\n");
     417                if(col_sglobaldoflist) for(i=0;i<col_ssize;i++) _printf_(" " << col_sglobaldoflist[i] << "\n"); _printf_(" \n");
    418418        }
    419419}

  • issm/trunk-jpl/src/c/classes/matrix/ElementVector.cpp

    r15100 r15104  
    219219        int i;
    220220
    221         _printf_("Element Vector echo:" << "\n");
     221        _printf_("Element Vector echo:\n");
    222222        _printf_("   nrows: " << nrows << "\n");
    223         _printf_("   values:" << "\n");
     223        _printf_("   values:\n");
    224224        for(i=0;i<nrows;i++){
    225225                _printf_(setw(4) << right << i << ": " << setw(10) << values[i] << "\n");
    … …  
    228228        _printf_("   gglobaldoflist (" << gglobaldoflist << "): ");
    229229        if(gglobaldoflist) for(i=0;i<nrows;i++) _printf_(" " << gglobaldoflist[i] );
    230         _printf_(" " << "\n");
     230        _printf_(" \n");
    231231
    232232        _printf_("   fsize: " << fsize << "\n");
    233233        _printf_("   flocaldoflist  (" << flocaldoflist << "): ");
    234234        if(flocaldoflist) for(i=0;i<fsize;i++) _printf_(" " << flocaldoflist[i] );
    235         _printf_(" " << "\n");
     235        _printf_(" \n");
    236236        _printf_("   fglobaldoflist (" << fglobaldoflist << "): ");
    237237        if(fglobaldoflist) for(i=0;i<fsize;i++) _printf_(" " << fglobaldoflist[i] );
    238         _printf_(" " << "\n");
     238        _printf_(" \n");
    239239}
    240240/*}}}*/

  • issm/trunk-jpl/src/c/datastructures/DataSet.cpp

    r15100 r15104  
    130130        if(this==NULL)_error_("trying to echo a NULL dataset");
    131131
    132         _printf0_("DataSet echo: " << objects.size() << " objects" << "\n");
     132        _printf0_("DataSet echo: " << objects.size() << " objects\n");
    133133
    134134        for ( object=objects.begin() ; object < objects.end(); object++ ){
    … …  
    147147        if(this==NULL)_error_("trying to echo a NULL dataset");
    148148
    149         _printf0_("DataSet echo: " << objects.size() << " objects" << "\n");
     149        _printf0_("DataSet echo: " << objects.size() << " objects\n");
    150150
    151151        for ( object=objects.begin() ; object < objects.end(); object++ ){

  • issm/trunk-jpl/src/c/kml/KMLFileReadUtils.cpp

    r15100 r15104  
    133133
    134134//      if      (itag)
    135 //              _printf0_("tag buffer (length=" << ibuf << "):" << "\n");
     135//              _printf0_("tag buffer (length=" << ibuf << "):\n");
    136136//      else if (ifield)
    137 //              _printf0_("field buffer (length=" << ibuf << "):" << "\n");
     137//              _printf0_("field buffer (length=" << ibuf << "):\n");
    138138//      _printf0_(buffer << "\n");
    139139
    … …  
    199199                }
    200200
    201 //      _printf0_("comment buffer (length=" << ibuf << "):" << "\n");
     201//      _printf0_("comment buffer (length=" << ibuf << "):\n");
    202202//      _printf0_(buffer << "\n");
    203203
    … …  
    265265
    266266        ktokn=strtok(ktagi,"< >");
    267 //      _printf0_("KMLFileTagName -- initial token=\"" << ktokn << "\"." << "\n");
     267//      _printf0_("KMLFileTagName -- initial token=\"" << ktokn << "\".\n");
    268268
    269269        if (!pname) {
    … …  
    275275
    276276        if (maxlen && (maxlen < strlen(ktokn))) {
    277                 _printf0_("KMLFileTagName -- string field too short for " << ktag << "." << "\n");
    278                 _printf0_("KMLFileTagName -- \"" << ktokn << "\" truncated to " << maxlen << " characters." << "\n");
     277                _printf0_("KMLFileTagName -- string field too short for " << ktag << ".\n");
     278                _printf0_("KMLFileTagName -- \"" << ktokn << "\" truncated to " << maxlen << " characters.\n");
    279279                strncpy(pname,ktokn,maxlen);
    280280        }
    … …  
    311311        /*  return first non blank and move past subsequent blank  */
    312312        ktokn=strtok(ktagi," ");
    313 //      _printf0_("KMLFileTagAttrib -- initial token=\"" << ktokn << "\"." << "\n");
     313//      _printf0_("KMLFileTagAttrib -- initial token=\"" << ktokn << "\".\n");
    314314
    315315        /*  return next non " =?/>" and move past subsequent " =?/>"  */
    … …  
    318318                /*  return next non quote and move past subsequent quote  */
    319319                ktokv=strtok(NULL,quote);
    320 //              _printf0_("KMLFileTagAttrib -- attribute " << ktokn << "=\"" << ktokv << "\"." << "\n");
     320//              _printf0_("KMLFileTagAttrib -- attribute " << ktokn << "=\"" << ktokv << "\".\n");
    321321
    322322/*  add the attribute to the dataset  */
    … …  
    334334                (!strncmp(&ktag[0],"<"        ,1) && !strncmp(&ktag[strlen(ktag)-2],"/>",2)))
    335335                isolo=1;
    336 //      _printf0_("KMLFileTagAttrib -- isolo=" << isolo << "." << "\n");
     336//      _printf0_("KMLFileTagAttrib -- isolo=" << isolo << ".\n");
    337337
    338338        return(isolo);
    … …  
    368368                        xDelete<char>(kstr);
    369369
    370 //      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pival << "." << "\n");
     370//      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pival << ".\n");
    371371
    372372        return(0);
    … …  
    402402                        xDelete<char>(kstr);
    403403
    404 //      _printf0_("KMLFileTokenParse -- " << ktag << "=" << (*pbval ? "true" : "false") << "." << "\n");
     404//      _printf0_("KMLFileTokenParse -- " << ktag << "=" << (*pbval ? "true" : "false") << ".\n");
    405405
    406406        return(0);
    … …  
    439439
    440440        if (maxlen && (maxlen < strlen(kstr))) {
    441                 _printf0_("KMLFileTokenParse -- string field too short for " << ktag << "." << "\n");
    442                 _printf0_("KMLFileTokenParse -- \"" << kstr << "\" truncated to " << maxlen << " characters." << "\n");
     441                _printf0_("KMLFileTokenParse -- string field too short for " << ktag << ".\n");
     442                _printf0_("KMLFileTokenParse -- \"" << kstr << "\" truncated to " << maxlen << " characters.\n");
    443443                strncpy(pstr,kstr,maxlen);
    444444        }
    … …  
    463463                        xDelete<char>(kstr);
    464464
    465 //      _printf0_("KMLFileTokenParse -- " << ktag << "=\"" << pstr << "\"." << "\n");
     465//      _printf0_("KMLFileTokenParse -- " << ktag << "=\"" << pstr << "\".\n");
    466466
    467467        return(pstr);
    … …  
    497497                        xDelete<char>(kstr);
    498498
    499 //      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pfval << "." << "\n");
     499//      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pfval << ".\n");
    500500
    501501        return(0);
    … …  
    531531                        xDelete<char>(kstr);
    532532
    533 //      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pdval << "." << "\n");
     533//      _printf0_("KMLFileTokenParse -- " << ktag << "=" << *pdval << ".\n");
    534534
    535535        return(0);
    … …  
    595595                        xDelete<char>(kstr);
    596596
    597 //      _printf0_("KMLFileTokenParse -- " << ktag << "=..." << "\n");
     597//      _printf0_("KMLFileTokenParse -- " << ktag << "=...\n");
    598598//      for (j=0; j<=i; j++)
    599 //              _printf0_("   [" << j << "]: " << (*pdval)[j] << "g" << "\n");
     599//              _printf0_("   [" << j << "]: " << (*pdval)[j] << "g\n");
    600600
    601601        return(0);
    … …  
    651651
    652652        if ((i+1) % n)
    653                 _printf0_("KMLFileTokenParse -- Double [m x n] field for " << ktag << " does not have multiple of n values." << "\n");
     653                _printf0_("KMLFileTokenParse -- Double [m x n] field for " << ktag << " does not have multiple of n values.\n");
    654654
    655655/*  get additional token and compare to closing tag  */
    … …  
    665665                        xDelete<char>(kstr);
    666666
    667 //      _printf0_("KMLFileTokenParse -- " << ktag << "=..." << "\n");
     667//      _printf0_("KMLFileTokenParse -- " << ktag << "=...\n");
    668668//      for (j=0; j<=i; j++)
    669 //              _printf0_("   [" << j << "]: " << (*pdval)[j] << "g" << "\n");
     669//              _printf0_("   [" << j << "]: " << (*pdval)[j] << "g\n");
    670670
    671671        return(0);
    … …  
    681681        opening tag, must find corresponding closing tag  */
    682682
    683         _printf0_("KMLFileTagSkip -- input tag " << ktag << "." << "\n");
     683        _printf0_("KMLFileTagSkip -- input tag " << ktag << ".\n");
    684684
    685685/*  if next token is a closing tag, compare to input  */
    … …  
    690690                                 (kstr[1] == '/') &&
    691691                                 (!strncmp(&(kstr[2]),&(ktag[1]),(strcspn(ktag," >")-1)/sizeof(char)))) {
    692                         _printf0_("KMLFileTagSkip -- closing tag " << kstr << "." << "\n");
     692                        _printf0_("KMLFileTagSkip -- closing tag " << kstr << ".\n");
    693693                        xDelete<char>(kstr);
    694694                        return(0);
    … …  
    699699                else if ((kstr[0] == '<') &&
    700700                                 (kstr[1] != '/')) {
    701                         _printf0_("KMLFileTagSkip -- opening tag " << kstr << "." << "\n");
     701                        _printf0_("KMLFileTagSkip -- opening tag " << kstr << ".\n");
    702702                        KMLFileTagSkip(kstr,
    703703                                                   fid);

  • issm/trunk-jpl/src/c/kml/KML_Attribute.cpp

    r15100 r15104  
    4343        for (i=0;i<10-strlen(name);i++)
    4444                if(flag) _printf0_(" ");
    45         if(flag) _printf0_(name << ": \"" << value << "\"" << "\n");
     45        if(flag) _printf0_(name << ": \"" << value << "\"\n");
    4646
    4747        return;
    … …  
    6767        for (i=0;i<10-strlen(name);i++)
    6868                if(flag) _printf0_(" ");
    69         if(flag) _printf0_(name << ": \"" << value << "\"" << "\n");
     69        if(flag) _printf0_(name << ": \"" << value << "\"\n");
    7070
    7171        return;

  • issm/trunk-jpl/src/c/kml/KML_Container.cpp

    r15100 r15104  
    4747        KML_Feature::Echo();
    4848
    49         if(flag) _printf0_("       feature: (size=" << feature->Size() << ")" << "\n");
     49        if(flag) _printf0_("       feature: (size=" << feature->Size() << ")\n");
    5050
    5151        return;
    … …  
    7878        if (feature->Size())
    7979                for (i=0; i<feature->Size(); i++) {
    80                         if(flag) _printf0_(indent << "       feature: -------- begin [" << i << "] --------" << "\n");
     80                        if(flag) _printf0_(indent << "       feature: -------- begin [" << i << "] --------\n");
    8181                        ((KML_Feature *)feature->GetObjectByOffset(i))->DeepEcho(indent2);
    82                         if(flag) _printf0_(indent << "       feature: --------  end  [" << i << "] --------" << "\n");
     82                        if(flag) _printf0_(indent << "       feature: --------  end  [" << i << "] --------\n");
    8383                }
    8484        else
    85                 if(flag) _printf0_(indent << "       feature: [empty]" << "\n");
     85                if(flag) _printf0_(indent << "       feature: [empty]\n");
    8686
    8787        return;

  • issm/trunk-jpl/src/c/kml/KML_Document.cpp

    r15100 r15104  
    3838        bool  flag=true;
    3939
    40         if(flag) _printf0_("KML_Document:" << "\n");
     40        if(flag) _printf0_("KML_Document:\n");
    4141        KML_Container::Echo();
    4242
    … …  
    5959        bool  flag=true;
    6060
    61         if(flag) _printf0_(indent << "KML_Document:" << "\n");
     61        if(flag) _printf0_(indent << "KML_Document:\n");
    6262        KML_Container::DeepEcho(indent);
    6363

  • issm/trunk-jpl/src/c/kml/KML_Feature.cpp

    r15100 r15104  
    6060        KML_Object::Echo();
    6161
    62         if(flag) _printf0_("          name: \"" << name << "\"" << "\n");
     62        if(flag) _printf0_("          name: \"" << name << "\"\n");
    6363        if(flag) _printf0_("    visibility: " << (visibility ? "true" : "false") << "\n");
    6464        if(flag) _printf0_("          open: " << (open ? "true" : "false") << "\n");
    65         if(flag) _printf0_("       snippet: \"" << snippet << "\"" << "\n");
    66         if(flag) _printf0_("      descript: \"" << descript << "\"" << "\n");
    67         if(flag) _printf0_("      styleurl: \"" << styleurl << "\"" << "\n");
    68         if(flag) _printf0_("         style: (size=" << style->Size() << ")" << "\n");
     65        if(flag) _printf0_("       snippet: \"" << snippet << "\"\n");
     66        if(flag) _printf0_("      descript: \"" << descript << "\"\n");
     67        if(flag) _printf0_("      styleurl: \"" << styleurl << "\"\n");
     68        if(flag) _printf0_("         style: (size=" << style->Size() << ")\n");
    6969
    7070        return;
    … …  
    9090        KML_Object::DeepEcho(indent);
    9191
    92         if(flag) _printf0_(indent << "          name: \"" << name << "\"" << "\n");
     92        if(flag) _printf0_(indent << "          name: \"" << name << "\"\n");
    9393        if(flag) _printf0_(indent << "    visibility: " << (visibility ? "true" : "false") << "\n");
    9494        if(flag) _printf0_(indent << "          open: " << (open ? "true" : "false") << "\n");
    95         if(flag) _printf0_(indent << "       snippet: \"" << snippet << "\"" << "\n");
    96         if(flag) _printf0_(indent << "      descript: \"" << descript << "\"" << "\n");
    97         if(flag) _printf0_(indent << "      styleurl: \"" << styleurl << "\"" << "\n");
     95        if(flag) _printf0_(indent << "       snippet: \"" << snippet << "\"\n");
     96        if(flag) _printf0_(indent << "      descript: \"" << descript << "\"\n");
     97        if(flag) _printf0_(indent << "      styleurl: \"" << styleurl << "\"\n");
    9898
    9999/*  loop over any styles for the feature  */
    … …  
    104104        if (style->Size())
    105105                for (i=0; i<style->Size(); i++) {
    106                         if(flag) _printf0_(indent << "         style: -------- begin [" << i << "] --------" << "\n");
     106                        if(flag) _printf0_(indent << "         style: -------- begin [" << i << "] --------\n");
    107107                        ((KML_Style *)style->GetObjectByOffset(i))->DeepEcho(indent2);
    108                         if(flag) _printf0_(indent << "         style: --------  end  [" << i << "] --------" << "\n");
     108                        if(flag) _printf0_(indent << "         style: --------  end  [" << i << "] --------\n");
    109109                }
    110110        else
    111                 if(flag) _printf0_(indent << "         style: [empty]" << "\n");
     111                if(flag) _printf0_(indent << "         style: [empty]\n");
    112112
    113113        return;

  • issm/trunk-jpl/src/c/kml/KML_File.cpp

    r15100 r15104  
    3939        bool  flag=true;
    4040
    41         if(flag) _printf0_("KML_File:" << "\n");
     41        if(flag) _printf0_("KML_File:\n");
    4242        KML_Object::Echo();
    4343
    … …  
    6060        bool  flag=true;
    6161
    62         if(flag) _printf0_(indent << "KML_File:" << "\n");
     62        if(flag) _printf0_(indent << "KML_File:\n");
    6363        KML_Object::DeepEcho(indent);
    6464

  • issm/trunk-jpl/src/c/kml/KML_Folder.cpp

    r15100 r15104  
    3838        bool  flag=true;
    3939
    40         if(flag) _printf0_("KML_Folder:" << "\n");
     40        if(flag) _printf0_("KML_Folder:\n");
    4141        KML_Container::Echo();
    4242
    … …  
    5959        bool  flag=true;
    6060
    61         if(flag) _printf0_(indent << "KML_Folder:" << "\n");
     61        if(flag) _printf0_(indent << "KML_Folder:\n");
    6262        KML_Container::DeepEcho(indent);
    6363

  • issm/trunk-jpl/src/c/kml/KML_GroundOverlay.cpp

    r15100 r15104  
    4444void  KML_GroundOverlay::Echo(){
    4545
    46         _printf_("KML_GroundOverlay:" << "\n");
     46        _printf_("KML_GroundOverlay:\n");
    4747        KML_Overlay::Echo();
    4848
    … …  
    6767        char  indent2[81];
    6868
    69         _printf_(indent << "KML_GroundOverlay:" << "\n");
     69        _printf_(indent << "KML_GroundOverlay:\n");
    7070        KML_Overlay::DeepEcho(indent);
    7171

  • issm/trunk-jpl/src/c/kml/KML_Icon.cpp

    r15100 r15104  
    4545        bool  flag=true;
    4646
    47         if(flag) _printf0_("KML_Icon:" << "\n");
     47        if(flag) _printf0_("KML_Icon:\n");
    4848        KML_Object::Echo();
    4949
    50         if(flag) _printf0_("          href: \"" << href << "\"" << "\n");
    51         if(flag) _printf0_("       refmode: \"" << refmode << "\"" << "\n");
     50        if(flag) _printf0_("          href: \"" << href << "\"\n");
     51        if(flag) _printf0_("       refmode: \"" << refmode << "\"\n");
    5252        if(flag) _printf0_("        refint: " << refint << "\n");
    53         if(flag) _printf0_("      vrefmode: \"" << vrefmode << "\"" << "\n");
     53        if(flag) _printf0_("      vrefmode: \"" << vrefmode << "\"\n");
    5454        if(flag) _printf0_("      vreftime: " << vreftime << "\n");
    5555        if(flag) _printf0_("      vboundsc: " << vboundsc << "\n");
    56         if(flag) _printf0_("       vformat: \"" << vformat << "\"" << "\n");
    57         if(flag) _printf0_("        hquery: \"" << hquery << "\"" << "\n");
     56        if(flag) _printf0_("       vformat: \"" << vformat << "\"\n");
     57        if(flag) _printf0_("        hquery: \"" << hquery << "\"\n");
    5858
    5959        return;
    … …  
    7575        bool  flag=true;
    7676
    77         if(flag) _printf0_(indent << "KML_Icon:" << "\n");
     77        if(flag) _printf0_(indent << "KML_Icon:\n");
    7878        KML_Object::DeepEcho(indent);
    7979
    80         if(flag) _printf0_(indent << "          href: \"" << href << "\"" << "\n");
    81         if(flag) _printf0_(indent << "       refmode: \"" << refmode << "\"" << "\n");
     80        if(flag) _printf0_(indent << "          href: \"" << href << "\"\n");
     81        if(flag) _printf0_(indent << "       refmode: \"" << refmode << "\"\n");
    8282        if(flag) _printf0_(indent << "        refint: " << refint << "\n");
    83         if(flag) _printf0_(indent << "      vrefmode: \"" << vrefmode << "\"" << "\n");
     83        if(flag) _printf0_(indent << "      vrefmode: \"" << vrefmode << "\"\n");
    8484        if(flag) _printf0_(indent << "      vreftime: " << vreftime << "\n");
    8585        if(flag) _printf0_(indent << "      vboundsc: " << vboundsc << "\n");
    86         if(flag) _printf0_(indent << "       vformat: \"" << vformat << "\"" << "\n");
    87         if(flag) _printf0_(indent << "        hquery: \"" << hquery << "\"" << "\n");
     86        if(flag) _printf0_(indent << "       vformat: \"" << vformat << "\"\n");
     87        if(flag) _printf0_(indent << "        hquery: \"" << hquery << "\"\n");
    8888
    8989        return;

  • issm/trunk-jpl/src/c/kml/KML_LatLonBox.cpp

    r15100 r15104  
    4040void  KML_LatLonBox::Echo(){
    4141
    42         _printf_("KML_LatLonBox:" << "\n");
     42        _printf_("KML_LatLonBox:\n");
    4343        KML_Object::Echo();
    4444
    … …  
    6363void  KML_LatLonBox::DeepEcho(const char* indent){
    6464
    65         _printf_(indent << "KML_LatLonBox:" << "\n");
     65        _printf_(indent << "KML_LatLonBox:\n");
    6666        KML_Object::DeepEcho(indent);
    6767

  • issm/trunk-jpl/src/c/kml/KML_LineString.cpp

    r15100 r15104  
    4646        bool  flag=true;
    4747
    48         if(flag) _printf0_("KML_LineString:" << "\n");
     48        if(flag) _printf0_("KML_LineString:\n");
    4949        KML_Geometry::Echo();
    5050
    5151        if(flag) _printf0_("       extrude: " << (extrude ? "true" : "false") << "\n");
    5252        if(flag) _printf0_("    tessellate: " << (tessellate ? "true" : "false") << "\n");
    53         if(flag) _printf0_("       altmode: \"" << altmode << "\"" << "\n");
    54         if(flag) _printf0_("        coords: (ncoord=" << ncoord << ")" << "\n");
     53        if(flag) _printf0_("       altmode: \"" << altmode << "\"\n");
     54        if(flag) _printf0_("        coords: (ncoord=" << ncoord << ")\n");
    5555
    5656        return;
    … …  
    7373        bool  flag=true;
    7474
    75         if(flag) _printf0_(indent << "KML_LineString:" << "\n");
     75        if(flag) _printf0_(indent << "KML_LineString:\n");
    7676        KML_Geometry::DeepEcho(indent);
    7777
    7878        if(flag) _printf0_(indent << "       extrude: " << (extrude ? "true" : "false") << "\n");
    7979        if(flag) _printf0_(indent << "    tessellate: " << (tessellate ? "true" : "false") << "\n");
    80         if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"" << "\n");
    81         if(flag) _printf0_(indent << "        coords: (ncoord=" << ncoord << ")" << "\n");
     80        if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"\n");
     81        if(flag) _printf0_(indent << "        coords: (ncoord=" << ncoord << ")\n");
    8282        for (i=0; i<ncoord; i++)
    83                 if(flag) _printf0_(indent << "                (" << coords[3*i+0] << "," << coords[3*i+1] << "," << coords[3*i+2] << ")" << "\n");
     83                if(flag) _printf0_(indent << "                (" << coords[3*i+0] << "," << coords[3*i+1] << "," << coords[3*i+2] << ")\n");
    8484
    8585        return;

  • issm/trunk-jpl/src/c/kml/KML_LineStyle.cpp

    r15100 r15104  
    3838        bool  flag=true;
    3939
    40         if(flag) _printf0_("KML_LineStyle:" << "\n");
     40        if(flag) _printf0_("KML_LineStyle:\n");
    4141        KML_ColorStyle::Echo();
    4242
    … …  
    6161        bool  flag=true;
    6262
    63         if(flag) _printf0_(indent << "KML_LineStyle:" << "\n");
     63        if(flag) _printf0_(indent << "KML_LineStyle:\n");
    6464        KML_ColorStyle::DeepEcho(indent);
    6565

  • issm/trunk-jpl/src/c/kml/KML_LinearRing.cpp

    r15100 r15104  
    4646        bool  flag=true;
    4747
    48         if(flag) _printf0_("KML_LinearRing:" << "\n");
     48        if(flag) _printf0_("KML_LinearRing:\n");
    4949        KML_Geometry::Echo();
    5050
    5151        if(flag) _printf0_("       extrude: " << (extrude ? "true" : "false") << "\n");
    5252        if(flag) _printf0_("    tessellate: " << (tessellate ? "true" : "false") << "\n");
    53         if(flag) _printf0_("       altmode: \"" << altmode << "\"" << "\n");
    54         if(flag) _printf0_("        coords: (ncoord=" << ncoord << ")" << "\n");
     53        if(flag) _printf0_("       altmode: \"" << altmode << "\"\n");
     54        if(flag) _printf0_("        coords: (ncoord=" << ncoord << ")\n");
    5555
    5656        return;
    … …  
    7373        bool  flag=true;
    7474
    75         if(flag) _printf0_(indent << "KML_LinearRing:" << "\n");
     75        if(flag) _printf0_(indent << "KML_LinearRing:\n");
    7676        KML_Geometry::DeepEcho(indent);
    7777
    7878        if(flag) _printf0_(indent << "       extrude: " << (extrude ? "true" : "false") << "\n");
    7979        if(flag) _printf0_(indent << "    tessellate: " << (tessellate ? "true" : "false") << "\n");
    80         if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"" << "\n");
    81         if(flag) _printf0_(indent << "        coords: (ncoord=" << ncoord << ")" << "\n");
     80        if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"\n");
     81        if(flag) _printf0_(indent << "        coords: (ncoord=" << ncoord << ")\n");
    8282        for (i=0; i<ncoord; i++)
    83                 if(flag)_printf_(indent << "                (" <<coords[3*i+0] << "," <<coords[3*i+1] << "," <<coords[3*i+2] << ")\n" << "\n");
     83                if(flag)_printf_(indent << "                (" <<coords[3*i+0] << "," <<coords[3*i+1] << "," <<coords[3*i+2] << ")\n\n");
    8484
    8585        return;

  • issm/trunk-jpl/src/c/kml/KML_MultiGeometry.cpp

    r15100 r15104  
    4646        bool  flag=true;
    4747
    48         if(flag) _printf0_("KML_Multigeometry:" << "\n");
     48        if(flag) _printf0_("KML_Multigeometry:\n");
    4949        KML_Geometry::Echo();
    5050
    51         if(flag) _printf0_("      geometry: (size=" << geometry->Size() << ")" << "\n");
     51        if(flag) _printf0_("      geometry: (size=" << geometry->Size() << ")\n");
    5252
    5353        return;
    … …  
    7171        bool  flag=true;
    7272
    73         if(flag) _printf0_(indent << "KML_Multigeometry:" << "\n");
     73        if(flag) _printf0_(indent << "KML_Multigeometry:\n");
    7474        KML_Geometry::DeepEcho(indent);
    7575
    … …  
    8181        if (geometry->Size())
    8282                for (i=0; i<geometry->Size(); i++) {
    83                         if(flag) _printf0_(indent << "      geometry: -------- begin [" << i << "] --------" << "\n");
     83                        if(flag) _printf0_(indent << "      geometry: -------- begin [" << i << "] --------\n");
    8484                        ((KML_Geometry *)geometry->GetObjectByOffset(i))->DeepEcho(indent2);
    85                         if(flag) _printf0_(indent << "      geometry: --------  end  [" << i << "] --------" << "\n");
     85                        if(flag) _printf0_(indent << "      geometry: --------  end  [" << i << "] --------\n");
    8686                }
    8787        else
    88                 if(flag) _printf0_(indent << "      geometry: [empty]" << "\n");
     88                if(flag) _printf0_(indent << "      geometry: [empty]\n");
    8989
    9090        return;

  • issm/trunk-jpl/src/c/kml/KML_Object.cpp

    r15100 r15104  
    6666        bool  flag=true;
    6767
    68         if(flag) _printf0_("        attrib: (size=" << attrib->Size() << ")" << "\n");
    69         if(flag) _printf0_("        commnt: (size=" << commnt->Size() << ")" << "\n");
    70         if(flag) _printf0_("        kmlobj: (size=" << kmlobj->Size() << ")" << "\n");
     68        if(flag) _printf0_("        attrib: (size=" << attrib->Size() << ")\n");
     69        if(flag) _printf0_("        commnt: (size=" << commnt->Size() << ")\n");
     70        if(flag) _printf0_("        kmlobj: (size=" << kmlobj->Size() << ")\n");
    7171
    7272        return;
    … …  
    9797                }
    9898        else
    99                 if(flag) _printf0_(indent << "        attrib: [empty]" << "\n");
     99                if(flag) _printf0_(indent << "        attrib: [empty]\n");
    100100
    101101/*  loop over the comments for the object  */
    … …  
    106106                }
    107107        else
    108                 if(flag) _printf0_(indent << "        commnt: [empty]" << "\n");
     108                if(flag) _printf0_(indent << "        commnt: [empty]\n");
    109109
    110110/*  loop over the unknown objects for the object  */
    … …  
    115115        if (kmlobj->Size())
    116116                for (i=0; i<kmlobj->Size(); i++) {
    117             if(flag) _printf0_(indent << "        kmlobj: -------- begin [" << i << "] --------" << "\n");
     117            if(flag) _printf0_(indent << "        kmlobj: -------- begin [" << i << "] --------\n");
    118118                        ((KML_Unknown *)kmlobj->GetObjectByOffset(i))->DeepEcho(indent2);
    119             if(flag) _printf0_(indent << "        kmlobj: --------  end  [" << i << "] --------" << "\n");
     119            if(flag) _printf0_(indent << "        kmlobj: --------  end  [" << i << "] --------\n");
    120120                }
    121121        else
    122                 if(flag) _printf0_(indent << "        kmlobj: [empty]" << "\n");
     122                if(flag) _printf0_(indent << "        kmlobj: [empty]\n");
    123123
    124124        return;
    … …  
    264264
    265265        else if (!strncmp(kstr,"<",1)) {
    266                 _printf0_("KML_Object::Read -- Unrecognized opening tag " << kstr << "." << "\n");
     266                _printf0_("KML_Object::Read -- Unrecognized opening tag " << kstr << ".\n");
    267267//              KMLFileTagSkip(kstr,
    268268//                                         fid);

  • issm/trunk-jpl/src/c/kml/KML_Overlay.cpp

    r15100 r15104  
    4545
    4646        KML_Feature::Echo();
    47         _printf0_("         color: \"" << color << "\"" << "\n");
     47        _printf0_("         color: \"" << color << "\"\n");
    4848        _printf0_("       draword: " << draword << "\n");
    4949        _printf0_("          icon: " << icon << "\n");

  • issm/trunk-jpl/src/c/kml/KML_Placemark.cpp

    r15100 r15104  
    4747        bool  flag=true;
    4848
    49         if(flag) _printf0_("KML_Placemark:" << "\n");
     49        if(flag) _printf0_("KML_Placemark:\n");
    5050        KML_Feature::Echo();
    5151
    52         if(flag) _printf0_("      geometry: (size=" << geometry->Size() << ")" << "\n");
     52        if(flag) _printf0_("      geometry: (size=" << geometry->Size() << ")\n");
    5353
    5454        return;
    … …  
    7272        bool  flag=true;
    7373
    74         if(flag) _printf0_(indent << "KML_Placemark:" << "\n");
     74        if(flag) _printf0_(indent << "KML_Placemark:\n");
    7575        KML_Feature::DeepEcho(indent);
    7676
    … …  
    8282        if (geometry->Size())
    8383                for (i=0; i<geometry->Size(); i++) {
    84                         if(flag) _printf0_(indent << "      geometry: -------- begin [" << i << "] --------" << "\n");
     84                        if(flag) _printf0_(indent << "      geometry: -------- begin [" << i << "] --------\n");
    8585                        ((KML_Geometry *)geometry->GetObjectByOffset(i))->DeepEcho(indent2);
    86                         if(flag) _printf0_(indent << "      geometry: --------  end  [" << i << "] --------" << "\n");
     86                        if(flag) _printf0_(indent << "      geometry: --------  end  [" << i << "] --------\n");
    8787                }
    8888        else
    89                 if(flag) _printf0_(indent << "      geometry: [empty]" << "\n");
     89                if(flag) _printf0_(indent << "      geometry: [empty]\n");
    9090
    9191        return;

  • issm/trunk-jpl/src/c/kml/KML_Point.cpp

    r15100 r15104  
    4343        bool  flag=true;
    4444
    45         if(flag) _printf0_("KML_Point:" << "\n");
     45        if(flag) _printf0_("KML_Point:\n");
    4646        KML_Geometry::Echo();
    4747
    4848        if(flag) _printf0_("       extrude: " << (extrude ? "true" : "false") << "\n");
    49         if(flag) _printf0_("       altmode: \"" << altmode << "\"" << "\n");
    50         if(flag) _printf0_("        coords: (" << coords[0] << "," << coords[1] << "," << coords[2] << ")" << "\n");
     49        if(flag) _printf0_("       altmode: \"" << altmode << "\"\n");
     50        if(flag) _printf0_("        coords: (" << coords[0] << "," << coords[1] << "," << coords[2] << ")\n");
    5151
    5252        return;
    … …  
    6868        bool  flag=true;
    6969
    70         if(flag) _printf0_(indent << "KML_Point:" << "\n");
     70        if(flag) _printf0_(indent << "KML_Point:\n");
    7171        KML_Geometry::DeepEcho(indent);
    7272
    7373        if(flag) _printf0_(indent << "       extrude: " << (extrude ? "true" : "false") << "\n");
    74         if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"" << "\n");
    75         if(flag) _printf0_(indent << "        coords: (" << coords[0] << "," << coords[1] << "," << coords[2] << ")" << "\n");
     74        if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"\n");
     75        if(flag) _printf0_(indent << "        coords: (" << coords[0] << "," << coords[1] << "," << coords[2] << ")\n");
    7676
    7777        return;

  • issm/trunk-jpl/src/c/kml/KML_PolyStyle.cpp

    r15100 r15104  
    4141        bool  flag=true;
    4242
    43         if(flag) _printf0_("KML_PolyStyle:" << "\n");
     43        if(flag) _printf0_("KML_PolyStyle:\n");
    4444        KML_ColorStyle::Echo();
    4545
    … …  
    6565        bool  flag=true;
    6666
    67         if(flag) _printf0_(indent << "KML_PolyStyle:" << "\n");
     67        if(flag) _printf0_(indent << "KML_PolyStyle:\n");
    6868        KML_ColorStyle::DeepEcho(indent);
    6969

  • issm/trunk-jpl/src/c/kml/KML_Polygon.cpp

    r15100 r15104  
    5252        bool  flag=true;
    5353
    54         if(flag) _printf0_("KML_Polygon:" << "\n");
     54        if(flag) _printf0_("KML_Polygon:\n");
    5555        KML_Geometry::Echo();
    5656
    5757        if(flag) _printf0_("       extrude: " << (extrude ? "true" : "false") << "\n");
    5858        if(flag) _printf0_("    tessellate: " << (tessellate ? "true" : "false") << "\n");
    59         if(flag) _printf0_("       altmode: \"" << altmode << "\"" << "\n");
    60         if(flag) _printf0_("         outer: (size=" << outer->Size() << ")" << "\n");
    61         if(flag) _printf0_("         inner: (size=" << inner->Size() << ")" << "\n");
     59        if(flag) _printf0_("       altmode: \"" << altmode << "\"\n");
     60        if(flag) _printf0_("         outer: (size=" << outer->Size() << ")\n");
     61        if(flag) _printf0_("         inner: (size=" << inner->Size() << ")\n");
    6262
    6363        return;
    … …  
    8181        bool  flag=true;
    8282
    83         if(flag) _printf0_(indent << "KML_Polygon:" << "\n");
     83        if(flag) _printf0_(indent << "KML_Polygon:\n");
    8484        KML_Geometry::DeepEcho(indent);
    8585
    8686        if(flag) _printf0_(indent << "       extrude: " << (extrude ? "true" : "false") << "\n");
    8787        if(flag) _printf0_(indent << "    tessellate: " << (tessellate ? "true" : "false") << "\n");
    88         if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"" << "\n");
     88        if(flag) _printf0_(indent << "       altmode: \"" << altmode << "\"\n");
    8989
    9090        memcpy(indent2,indent,(strlen(indent)+1)*sizeof(char));
    … …  
    9393        if (outer->Size())
    9494                for (i=0; i<outer->Size(); i++) {
    95                         if(flag) _printf0_(indent << "         outer: -------- begin [" << i << "] --------" << "\n");
     95                        if(flag) _printf0_(indent << "         outer: -------- begin [" << i << "] --------\n");
    9696                        ((KML_LinearRing *)outer->GetObjectByOffset(i))->DeepEcho(indent2);
    97                         if(flag) _printf0_(indent << "         outer: --------  end  [" << i << "] --------" << "\n");
     97                        if(flag) _printf0_(indent << "         outer: --------  end  [" << i << "] --------\n");
    9898                }
    9999        else
    100                 if(flag) _printf0_(indent << "         outer: [empty]" << "\n");
     100                if(flag) _printf0_(indent << "         outer: [empty]\n");
    101101
    102102        if (inner->Size())
    103103                for (i=0; i<inner->Size(); i++) {
    104                         if(flag) _printf0_(indent << "         inner: -------- begin [" << i << "] --------" << "\n");
     104                        if(flag) _printf0_(indent << "         inner: -------- begin [" << i << "] --------\n");
    105105                        ((KML_LinearRing *)inner->GetObjectByOffset(i))->DeepEcho(indent2);
    106                         if(flag) _printf0_(indent << "         inner: --------  end  [" << i << "] --------" << "\n");
     106                        if(flag) _printf0_(indent << "         inner: --------  end  [" << i << "] --------\n");
    107107                }
    108108        else
    109                 if(flag) _printf0_(indent << "         inner: [empty]" << "\n");
     109                if(flag) _printf0_(indent << "         inner: [empty]\n");
    110110
    111111        return;

  • issm/trunk-jpl/src/c/kml/KML_Style.cpp

    r15100 r15104  
    6868        bool  flag=true;
    6969
    70         if(flag) _printf0_("KML_Style:" << "\n");
     70        if(flag) _printf0_("KML_Style:\n");
    7171        KML_StyleSelector::Echo();
    7272
    … …  
    9797        bool  flag=true;
    9898
    99         if(flag) _printf0_(indent << "KML_Style:" << "\n");
     99        if(flag) _printf0_(indent << "KML_Style:\n");
    100100        KML_StyleSelector::DeepEcho(indent);
    101101

  • issm/trunk-jpl/src/c/kml/KML_Unknown.cpp

    r15100 r15104  
    4040        bool  flag=true;
    4141
    42         if(flag) _printf0_("KML_Unknown " << name << ":" << "\n");
     42        if(flag) _printf0_("KML_Unknown " << name << ":\n");
    4343        KML_Object::Echo();
    4444
    4545        if (value     )
    46                 if(flag) _printf0_("         value: \"" << value << "\"" << "\n");
     46                if(flag) _printf0_("         value: \"" << value << "\"\n");
    4747    else
    48         if(flag) _printf0_("         value: [none]" << "\n");
     48        if(flag) _printf0_("         value: [none]\n");
    4949
    5050        return;
    … …  
    6969        bool         flag=true;
    7070
    71         if(flag) _printf0_(indent << "KML_Unknown " << name << ":" << "\n");
     71        if(flag) _printf0_(indent << "KML_Unknown " << name << ":\n");
    7272        KML_Object::DeepEcho(indent);
    7373
    … …  
    8181                while (vtoken=strtok(NULL,nl))
    8282                        if(flag) _printf0_("\n" << indent << "                 " << vtoken);
    83                 if(flag) _printf0_("\"" << "\n");
     83                if(flag) _printf0_("\"\n");
    8484
    8585                xDelete<char>(valuei);
    8686        }
    8787    else
    88         if(flag) _printf0_(indent << "         value: [none]" << "\n");
     88        if(flag) _printf0_(indent << "         value: [none]\n");
    8989
    9090        return;
    … …  
    136136        name=KMLFileTagName(NULL,
    137137                                                kstr);
    138 //      _printf0_("KML_Unknown::Read -- opening name=" << name << "." << "\n");
     138//      _printf0_("KML_Unknown::Read -- opening name=" << name << ".\n");
    139139
    140140/*  get object attributes and check for solo tag  */
    … …  
    148148        while (kstri=KMLFileToken(fid,
    149149                                                          &ncom,&pcom)) {
    150 //              _printf0_("KML_Unknown::Read -- kstri=" << kstri << "." << "\n");
     150//              _printf0_("KML_Unknown::Read -- kstri=" << kstri << ".\n");
    151151                if      (!strncmp(&kstri[0],"</", 2) &&
    152152                                 !strncmp(&kstri[2],name,strlen(name))) {
    153 //                      _printf0_("KML_Unknown::Read -- closing name=" << name << "." << "\n");
     153//                      _printf0_("KML_Unknown::Read -- closing name=" << name << ".\n");
    154154                        xDelete<char>(kstri);
    155155                        break;

  • issm/trunk-jpl/src/c/main/kriging.cpp

    r15100 r15104  
    5151        pfclose(input_fid,binfilename);
    5252
    53         _printf0_("call computational core:" << "\n");
     53        _printf0_("call computational core:\n");
    5454        pKrigingx(&predictions,&error,x,y,data,nobs,x_interp,y_interp,ninterp,options);
    5555
    56         _printf0_("write results to disk:" << "\n");
     56        _printf0_("write results to disk:\n");
    5757        Results *results = new Results();
    5858        if(IssmComm::GetRank()==0){
    … …  
    6868
    6969        /*Close output and toolkits options file and write lock file if requested*/
    70         _printf0_("write lock file:" << "\n");
     70        _printf0_("write lock file:\n");
    7171        WriteLockFile(lockfilename);
    7272

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

    r15100 r15104  
    4040
    4141                //Step1: generate geometry Gh
    42                 if (verbosity>0) _printf_("Construction of a mesh from a given geometry" << "\n");
    43                 if (verbosity>1) _printf_("   Processing geometry..." << "\n");
     42                if (verbosity>0) _printf_("Construction of a mesh from a given geometry\n");
     43                if (verbosity>1) _printf_("   Processing geometry...\n");
    4444                Geometry Gh(bamggeom_in,bamgopts);
    4545
    … …  
    4949
    5050                //build metric using geometry
    51                 if (verbosity>1) _printf_("   Generating Metric..." << "\n");
     51                if (verbosity>1) _printf_("   Generating Metric...\n");
    5252                for(i=0;i<Gh.nbv;i++){
    5353                        Metric M=Gh[i];
    … …  
    5959
    6060                //generate mesh
    61                 if (verbosity>1) _printf_("   Generating Mesh..." << "\n");
     61                if (verbosity>1) _printf_("   Generating Mesh...\n");
    6262                Mesh Th(maxnbv,Gh,bamgopts);
    6363
    … …  
    7272
    7373                //Build output
    74                 if (verbosity>1) _printf_("   Write Mesh..." << "\n");
     74                if (verbosity>1) _printf_("   Write Mesh...\n");
    7575                Th.WriteMesh(bamgmesh_out,bamgopts);
    76                 if (verbosity>1) _printf_("   Write Geometry..." << "\n");
     76                if (verbosity>1) _printf_("   Write Geometry...\n");
    7777                Gh.WriteGeometry(bamggeom_out,bamgopts);
    7878
    … …  
    8686
    8787                // read background mesh
    88                 if (verbosity>0) _printf_("Anisotropic mesh adaptation" << "\n");
    89                 if (verbosity>1) _printf_("   Processing initial mesh and geometry..." << "\n");
     88                if (verbosity>0) _printf_("Anisotropic mesh adaptation\n");
     89                if (verbosity>1) _printf_("   Processing initial mesh and geometry...\n");
    9090                Mesh BTh(bamggeom_in,bamgmesh_in,bamgopts);
    9191
    … …  
    9999                //Generate initial metric
    100100                if (bamgopts->metric){
    101                         if (verbosity>1) _printf_("   Processing Metric..." << "\n");
     101                        if (verbosity>1) _printf_("   Processing Metric...\n");
    102102                        BTh.ReadMetric(bamgopts);
    103103                }
    104104                else {
    105                         if (verbosity>1) _printf_("   Generating initial Metric..." << "\n");
     105                        if (verbosity>1) _printf_("   Generating initial Metric...\n");
    106106                        Metric Mhmax(bamgopts->hmax);
    107107                        for (int iv=0;iv<BTh.nbv;iv++) BTh[iv].m = Mhmax;
    … …  
    110110                //use present fields to generate metric if present
    111111                if (bamgopts->field){
    112                         if (verbosity>1) _printf_("   Merge metric with field provided..." << "\n");
     112                        if (verbosity>1) _printf_("   Merge metric with field provided...\n");
    113113                        BTh.AddMetric(bamgopts);
    114114                }
    … …  
    116116                // change using hVertices if provided
    117117                if(bamgopts->hVertices && bamgopts->hVerticesSize[0]==BTh.nbv){
    118                         if (verbosity>1) _printf_("   Merging Metric with hVertices..." << "\n");
     118                        if (verbosity>1) _printf_("   Merging Metric with hVertices...\n");
    119119                        for (i=0;i<BTh.nbv;i++){
    120120                                if (!xIsNan<IssmPDouble>(bamgopts->hVertices[i])){
    … …  
    126126                // change using hminVertices if provided
    127127                if (bamgopts->hminVertices){
    128                         if (verbosity>1) _printf_("   Merging Metric with hminVertices..." << "\n");
     128                        if (verbosity>1) _printf_("   Merging Metric with hminVertices...\n");
    129129                        for (i=0;i<BTh.nbv;i++){
    130130                                if (!xIsNan<IssmPDouble>(bamgopts->hminVertices[i])){
    … …  
    139139                // change using hmaxVertices if provided
    140140                if (bamgopts->hmaxVertices){
    141                         if (verbosity>1) _printf_("   Merging Metric with hmaxVertices..." << "\n");
     141                        if (verbosity>1) _printf_("   Merging Metric with hmaxVertices...\n");
    142142                        for (i=0;i<BTh.nbv;i++){
    143143                                if (!xIsNan<IssmPDouble>(bamgopts->hmaxVertices[i])){
    … …  
    163163
    164164                //Build new mesh
    165                 if (verbosity>1) _printf_("   Generating Mesh..." << "\n");
     165                if (verbosity>1) _printf_("   Generating Mesh...\n");
    166166                Thr=&BTh,Thb=0;
    167167                Mesh & Th( *(0 ?  new Mesh(*Thr,&Thr->Gh,Thb,maxnbv) :  new Mesh(maxnbv,BTh,bamgopts,bamgopts->KeepVertices)));
    … …  
    192192
    193193                //Build output
    194                 if (verbosity>1) _printf_("   Write Mesh..." << "\n");
     194                if (verbosity>1) _printf_("   Write Mesh...\n");
    195195                Th.WriteMesh(bamgmesh_out,bamgopts);
    196                 if (verbosity>1) _printf_("   Write Geometry..." << "\n");
     196                if (verbosity>1) _printf_("   Write Geometry...\n");
    197197                Th.Gh.WriteGeometry(bamggeom_out,bamgopts);
    198                 if (verbosity>1) _printf_("   Write Metric..." << "\n");
     198                if (verbosity>1) _printf_("   Write Metric...\n");
    199199                BTh.WriteMetric(bamgopts);
    200200
    … …  
    206206
    207207        /*No error return*/
    208         if (verbosity>1) _printf_("   Exiting Bamg." << "\n");
     208        if (verbosity>1) _printf_("   Exiting Bamg.\n");
    209209        return noerr;
    210210

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

    r15100 r15104  
    6363
    6464        if (DEBUG_TRACE > 0) {
    65                 _printf_("<Entering main>" << "\n");
     65                _printf_("<Entering main>\n");
    6666        }
    6767
    6868        if (PRINT_HEADERS) {
    69                 _printf_("\n                    Chaco 2.0" << "\n");
    70                 _printf_("          Sandia National Laboratories\n" << "\n");
     69                _printf_("\n                    Chaco 2.0\n");
     70                _printf_("          Sandia National Laboratories\n\n");
    7171        }
    7272
    … …  
    166166
    167167        if (DEBUG_MEMORY > 0) {
    168                 _printf_("" << "\n");
     168                _printf_("\n");
    169169                smalloc_stats();
    170170        }
    … …  
    174174
    175175        if (DEBUG_TRACE > 1) {
    176                 _printf_("<Leaving main>" << "\n");
     176                _printf_("<Leaving main>\n");
    177177        }
    178178

  • issm/trunk-jpl/src/c/modules/Chacox/input_parse.cpp

    r15100 r15104  
    3737
    3838        if (DEBUG_TRACE > 0) {
    39                 _printf_("<Entering input_parse>" << "\n");
     39                _printf_("<Entering input_parse>\n");
    4040        }
    4141
    4242        if (PROMPT) {
    43                 _printf_("Parallel machine architecture:" << "\n");
    44                 _printf_("  (0) Hypercube" << "\n");
    45                 _printf_("  (1) One-dimensional mesh" << "\n");
    46                 _printf_("  (2) Two-dimensional mesh" << "\n");
    47                 _printf_("  (3) Three-dimensional mesh" << "\n");
     43                _printf_("Parallel machine architecture:\n");
     44                _printf_("  (0) Hypercube\n");
     45                _printf_("  (1) One-dimensional mesh\n");
     46                _printf_("  (2) Two-dimensional mesh\n");
     47                _printf_("  (3) Three-dimensional mesh\n");
    4848        }
    4949        *architecture = (int)options[OPT_ARCH];
    … …  
    7373        else {
    7474                if (PROMPT) {
    75                         _printf_("Global partitioning method:" << "\n");
    76                         _printf_("  (1) Multilevel-KL" << "\n");
    77                         _printf_("  (2) Spectral" << "\n");
    78                         _printf_("  (3) Inertial" << "\n");
    79                         _printf_("  (4) Linear" << "\n");
    80                         _printf_("  (5) Random" << "\n");
    81                         _printf_("  (6) Scattered" << "\n");
    82                         _printf_("  (7) Read-from-file" << "\n");
     75                        _printf_("Global partitioning method:\n");
     76                        _printf_("  (1) Multilevel-KL\n");
     77                        _printf_("  (2) Spectral\n");
     78                        _printf_("  (3) Inertial\n");
     79                        _printf_("  (4) Linear\n");
     80                        _printf_("  (5) Random\n");
     81                        _printf_("  (6) Scattered\n");
     82                        _printf_("  (7) Read-from-file\n");
    8383                }
    8484                *global_method = (int)options[OPT_GLOBAL];
    … …  
    101101        else if (*global_method == 2) {
    102102                if (PROMPT) {
    103                         _printf_("Eigensolver:" << "\n");
    104                         _printf_("  (1) Multilevel RQI/Symmlq" << "\n");
    105                         _printf_("  (2) Lanczos" << "\n");
     103                        _printf_("Eigensolver:\n");
     104                        _printf_("  (1) Multilevel RQI/Symmlq\n");
     105                        _printf_("  (2) Lanczos\n");
    106106                }
    107107                eigensolver = (int)options[OPT_RQI];
    … …  
    162162        else {
    163163                if (PROMPT) {
    164                         _printf_("Local refinement method:" << "\n");
    165                         _printf_("  (1) Kernighan-Lin" << "\n");
    166                         _printf_("  (2) None" << "\n");
     164                        _printf_("Local refinement method:\n");
     165                        _printf_("  (1) Kernighan-Lin\n");
     166                        _printf_("  (2) None\n");
    167167                }
    168168                *local_method = (int)options[OPT_LOCAL];
    … …  
    181181                *ndims_tot = nparts[0];
    182182                if (*ndims_tot < 1) {
    183                         _printf_(" Number of divisions must be at least 1" << "\n");
     183                        _printf_(" Number of divisions must be at least 1\n");
    184184                        printf("%s -- Number of divisions %d must be at least 1.\n",
    185185                                   __FUNCT__,nparts[0]);
    … …  
    220220        else if (*nprocs <= 7) {
    221221                if (PROMPT) {
    222                         _printf_("Partitioning dimension: " << "\n");
    223                         _printf_("  (1) Bisection" << "\n");
    224                         _printf_("  (2) Quadrisection" << "\n");
     222                        _printf_("Partitioning dimension: \n");
     223                        _printf_("  (1) Bisection\n");
     224                        _printf_("  (2) Quadrisection\n");
    225225                }
    226226                *ndims = (int)options[OPT_NDIMS];
    … …  
    232232        else {
    233233                if (PROMPT) {
    234                         _printf_("Partitioning dimension: " << "\n");
    235                         _printf_("  (1) Bisection" << "\n");
    236                         _printf_("  (2) Quadrisection" << "\n");
    237                         _printf_("  (3) Octasection" << "\n");
     234                        _printf_("Partitioning dimension: \n");
     235                        _printf_("  (1) Bisection\n");
     236                        _printf_("  (2) Quadrisection\n");
     237                        _printf_("  (3) Octasection\n");
    238238                }
    239239                *ndims = (int)options[OPT_NDIMS];

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

    r15100 r15104  
    2323        parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
    2424
    25         if(VerboseMProcessor()) _printf0_("      Configuring elements..." << "\n");
     25        if(VerboseMProcessor()) _printf0_("      Configuring elements...\n");
    2626        for (i=0;i<elements->Size();i++){
    2727                element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
    2828                element->Configure(elements,loads,nodes,vertices,materials,parameters);
    2929        }
    30         if(VerboseMProcessor()) _printf0_("      Configuring loads..." << "\n");
     30        if(VerboseMProcessor()) _printf0_("      Configuring loads...\n");
    3131        for (i=0;i<loads->Size();i++){
    3232                load=(Load*)loads->GetObjectByOffset(i);
    … …  
    3535                }
    3636        }
    37         if(VerboseMProcessor()) _printf0_("      Configuring nodes..." << "\n");
     37        if(VerboseMProcessor()) _printf0_("      Configuring nodes...\n");
    3838        for (i=0;i<nodes->Size();i++){
    3939                node=(Node*)nodes->GetObjectByOffset(i);
    … …  
    4343        }
    4444
    45         if(VerboseMProcessor()) _printf0_("      Configuring materials..." << "\n");
     45        if(VerboseMProcessor()) _printf0_("      Configuring materials...\n");
    4646        for (i=0;i<materials->Size();i++){
    4747                material=(Material*)materials->GetObjectByOffset(i);

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

    r15100 r15104  
    1717
    1818        /*Display message*/
    19         if(VerboseModule()) _printf0_("   Constraining penalties" << "\n");
     19        if(VerboseModule()) _printf0_("   Constraining penalties\n");
    2020
    2121        /*recover parameters: */

  • issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp

    r15100 r15104  
    5151        }
    5252        else if(num_unstable_constraints<=min_mechanical_constraints){
    53                 if(VerboseModule()) _printf0_("   freezing constraints" << "\n");
     53                if(VerboseModule()) _printf0_("   freezing constraints\n");
    5454                RiftFreezeConstraints(loads,configuration_type);
    5555        }

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

    r15100 r15104  
    1717        int slocalsize;
    1818
    19         if(VerboseModule()) _printf0_("   Create nodal constraints" << "\n");
     19        if(VerboseModule()) _printf0_("   Create nodal constraints\n");
    2020
    2121        /*figure out how many dofs we have: */

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

    r15100 r15104  
    5353        if (!ExpRead(&nprof,&pnvert,&pprofx,&pprofy,&closed,filexp))
    5454                _error_("Error reading exp file.");
    55         _printf0_("Exp2Kmlx -- Reading " << nprof << " exp profiles from file \"" << filexp << "\"." << "\n");
     55        _printf0_("Exp2Kmlx -- Reading " << nprof << " exp profiles from file \"" << filexp << "\".\n");
    5656//      for (i=0; i<nprof; i++)
    5757//              _printf_("i=" << i << "; nvert=" << pnvert[i] << ", closed=" << closed[i] << "\n");
    … …  
    118118
    119119        if (holes && nprof && (pnvert[0] <= 1 || pprofx[0][pnvert[0]-1] != pprofx[0][0] || pprofy[0][pnvert[0]-1] != pprofy[0][0])) {
    120                 _printf0_("Warning -- Outer profile is not closed, so \"holes\" option will be ignored." << "\n");
     120                _printf0_("Warning -- Outer profile is not closed, so \"holes\" option will be ignored.\n");
    121121                holes=false;
    122122        }
    … …  
    150150                for (i=1; i<nprof; i++) {
    151151                        if (pnvert[i] <= 1 || pprofx[i][pnvert[i]-1] != pprofx[i][0] || pprofy[i][pnvert[i]-1] != pprofy[i][0]) {
    152                                 _printf0_("Warning -- Inner profile " << i+1 << " is not closed with \"holes\" specified, so it will be ignored." << "\n");
     152                                _printf0_("Warning -- Inner profile " << i+1 << " is not closed with \"holes\" specified, so it will be ignored.\n");
    153153                                continue;
    154154                        }
    … …  
    272272/*  write kml file  */
    273273
    274         _printf0_("Exp2Kmlx -- Writing kml document to file \"" << filkml << "\"." << "\n");
     274        _printf0_("Exp2Kmlx -- Writing kml document to file \"" << filkml << "\".\n");
    275275        fid=fopen(filkml,"w");
    276276        fprintf(fid,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
    … …  
    287287        clock1=clock();
    288288        time1 =time(NULL);
    289         _printf_("Exp2Kmlx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     289        _printf_("Exp2Kmlx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0) << " elapsed seconds.\n\n\n");
    290290
    291291        return(iret);

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

    r15100 r15104  
    1919        Vector<IssmDouble>* solution=NULL;
    2020
    21         if(VerboseModule()) _printf0_("   Get solution from inputs" << "\n");
     21        if(VerboseModule()) _printf0_("   Get solution from inputs\n");
    2222
    2323        /*retrive parameters: */

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

    r15100 r15104  
    1717        Element            *element                          = NULL;
    1818
    19         if(VerboseModule()) _printf0_("   Migrating grounding line" << "\n");
     19        if(VerboseModule()) _printf0_("   Migrating grounding line\n");
    2020
    2121        /*retrieve parameters: */

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

    r15100 r15104  
    268268
    269269                        #ifdef _DEBUG2_
    270                                 //_printf_(temp << " " << elev << " " << range << " " << "\n");
     270                                //_printf_(temp << " " << elev << " " << range << " \n");
    271271                        #endif
    272272
    … …  
    338338        #ifdef _DEBUG2_
    339339                _printf_( "\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b" );
    340                 _printf_("Number of zeroes remaining:          0\n" << "\n");
     340                _printf_("Number of zeroes remaining:          0\n\n");
    341341                printf ( "\n");
    342342        #endif

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

    r15100 r15104  
    1010
    1111        int i;
    12         if(VerboseModule()) _printf0_("   Input updates from constant" << "\n");
     12        if(VerboseModule()) _printf0_("   Input updates from constant\n");
    1313
    1414        /*Elements and loads drive the update: */
    … …  
    3131
    3232        int i;
    33         if(VerboseModule()) _printf0_("   Input updates from constant" << "\n");
     33        if(VerboseModule()) _printf0_("   Input updates from constant\n");
    3434
    3535        /*Elements and loads drive the update: */
    … …  
    5252
    5353        int i;
    54         if(VerboseModule()) _printf0_("   Input updates from constant" << "\n");
     54        if(VerboseModule()) _printf0_("   Input updates from constant\n");
    5555
    5656        /*Elements and loads drive the update: */

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

    r15100 r15104  
    8787        /*launch the thread manager with InterpFromGridToMeshxt as a core: */
    8888        LaunchThread(InterpFromGridToMeshxt,(void*)&gate,num);
    89         _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     89        _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    9090
    9191        /*Assign output pointers:*/
    … …  
    176176                                        break;
    177177                                default:
    178                                         _printf_("Interpolation " << EnumToStringx(interpenum) << " not supported yet" << "\n");
     178                                        _printf_("Interpolation " << EnumToStringx(interpenum) << " not supported yet\n");
    179179                                        return NULL; /*WARNING: no error because it would blow up the multithreading!*/
    180180                        }

  • issm/trunk-jpl/src/c/modules/InterpFromMesh2dx/InterpFromMesh2dxt.cpp

    r15100 r15104  
    9999        }
    100100        if(debug && my_thread==0)
    101          _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     101         _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    102102        return NULL;
    103103}

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

    r15100 r15104  
    168168        }
    169169        if (debug)
    170          _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     170         _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    171171
    172172        /*Assign output pointers:*/

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

    r15100 r15104  
    142142                }
    143143        }
    144         //if(N_interp>=100) _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     144        //if(N_interp>=100) _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    145145
    146146        /*clean-up and return*/

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

    r15100 r15104  
    132132        }
    133133        if (debug)
    134          _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     134         _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    135135
    136136        /*Assign output pointers:*/

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

    r15100 r15104  
    4545
    4646        if (kxml) {
    47                 _printf0_("XML declaration:" << "\n");
     47                _printf0_("XML declaration:\n");
    4848                kxml->DeepEcho("  ");
    4949                delete kxml;
    5050        }
    5151        if (kdtd) {
    52                 _printf0_("DTD declaration (not yet implemented):" << "\n");
     52                _printf0_("DTD declaration (not yet implemented):\n");
    5353                kdtd->DeepEcho("  ");
    5454                delete kdtd;
    … …  
    5757        clock1=clock();
    5858        time1 =time(NULL);
    59         _printf_("KMLFileReadx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     59        _printf_("KMLFileReadx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0) << " elapsed seconds.\n\n\n");
    6060
    6161        return(kfil);

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

    r15100 r15104  
    8080
    8181        if (cmap) {
    82                 _printf0_("Writing " << mcmap << " Matlab colors as KML style templates." << "\n");
     82                _printf0_("Writing " << mcmap << " Matlab colors as KML style templates.\n");
    8383                ipt=0;
    8484                for (i=0; i<mcmap; i++) {
    … …  
    115115
    116116        if (!nodecon) {
    117                 _printf0_("Creating the node connectivity table." << "\n");
     117                _printf0_("Creating the node connectivity table.\n");
    118118                nncon=mxepg+1;
    119119                nodecon=xNewZeroInit<int>(mncon*nncon);
    … …  
    146146
    147147                else if (mdata == mncon) {
    148                         _printf0_("Averaging nodal data to element data." << "\n");
     148                        _printf0_("Averaging nodal data to element data.\n");
    149149                        edata=xNewZeroInit<double>(melem*ndata);
    150150                        edfree=true;
    … …  
    185185        clock0a=clock();
    186186        time0a =time(NULL);
    187         _printf_("  Constructed kml document -- " << ((double)(clock0a-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0a,time0) << " elapsed seconds.\n\n" << "\n");
     187        _printf_("  Constructed kml document -- " << ((double)(clock0a-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0a,time0) << " elapsed seconds.\n\n\n");
    188188
    189189/*  write kml file  */
    190190
    191         _printf0_("Writing kml document to file." << "\n");
     191        _printf0_("Writing kml document to file.\n");
    192192        fprintf(fid,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
    193193        fprintf(fid,"<kml xmlns=\"http://www.opengis.net/kml/2.2\">\n");
    … …  
    196196        clock0b=clock();
    197197        time0b =time(NULL);
    198         _printf_("  Wrote kml file -- " << ((double)(clock0b-clock0a))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0b,time0a) << " elapsed seconds.\n\n" << "\n");
    199 
    200         _printf0_("Deleting kml document." << "\n");
     198        _printf_("  Wrote kml file -- " << ((double)(clock0b-clock0a))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0b,time0a) << " elapsed seconds.\n\n\n");
     199
     200        _printf0_("Deleting kml document.\n");
    201201        delete kdoc;
    202202        clock0c=clock();
    203203        time0c =time(NULL);
    204         _printf_("  Deleted kml document -- " << ((double)(clock0c-clock0b))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0c,time0b) << " elapsed seconds.\n\n" << "\n");
     204        _printf_("  Deleted kml document -- " << ((double)(clock0c-clock0b))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time0c,time0b) << " elapsed seconds.\n\n\n");
    205205
    206206        clock1=clock();
    207207        time1 =time(NULL);
    208         _printf_("KMLMeshWritex Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     208        _printf_("KMLMeshWritex Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n\n");
    209209
    210210        return;
    … …  
    297297/*  write each element as a polygon placemark  */
    298298
    299         _printf0_("Writing " << melem << " tria elements as KML polygons." << "\n");
     299        _printf0_("Writing " << melem << " tria elements as KML polygons.\n");
    300300
    301301        for (i=0; i<melem; i++) {
    … …  
    349349
    350350//              if (!(int)fmod((double)(i+1),1000))
    351 //                      _printf0_("  " << (i+1) << " tria elements written." << "\n");
     351//                      _printf0_("  " << (i+1) << " tria elements written.\n");
    352352        }
    353         _printf0_("  " << melem << " tria elements written." << "\n");
     353        _printf0_("  " << melem << " tria elements written.\n");
    354354
    355355        return(kfold);

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

    r15100 r15104  
    7979/*  write kml file  */
    8080
    81         _printf0_("Writing kml document to file." << "\n");
     81        _printf0_("Writing kml document to file.\n");
    8282        fprintf(fid,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
    8383        kfile->Write(fid,indent);
    … …  
    8787        clock1=clock();
    8888        time1 =time(NULL);
    89         _printf_("KMLOverlayx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     89        _printf_("KMLOverlayx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n\n");
    9090
    9191        return;

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

    r15100 r15104  
    4646/*  open exp file  */
    4747
    48         _printf0_("Writing exp profiles to file." << "\n");
     48        _printf0_("Writing exp profiles to file.\n");
    4949        fido=fopen(filexp,"w");
    5050
    … …  
    6161        clock1=clock();
    6262        time1 =time(NULL);
    63         _printf_("Kml2Expx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     63        _printf_("Kml2Expx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n\n");
    6464
    6565        return(iret);

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

    r15100 r15104  
    8383                observations->ObservationList(&x,&y,&data,&nobs);
    8484
    85                 _printf_("Generation Delaunay Triangulation" << "\n");
     85                _printf_("Generation Delaunay Triangulation\n");
    8686                BamgTriangulatex(&index,&nel,x,y,nobs);
    8787
    88                 _printf_("Interpolating" << "\n");
     88                _printf_("Interpolating\n");
    8989                xDelete<double>(predictions);
    9090                InterpFromMeshToMesh2dx(&predictions,index,x,y,nobs,nel,data,nobs,1,x_interp,y_interp,n_interp,options);
    … …  
    110110                /*launch the thread manager with Krigingxt as a core: */
    111111                LaunchThread(NearestNeighbort,(void*)&gate,num);
    112                 _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     112                _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    113113                xDelete<int>(gate.numdone);
    114114        }
    … …  
    132132                /*launch the thread manager with Krigingxt as a core: */
    133133                LaunchThread(idwt,(void*)&gate,num);
    134                 _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     134                _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    135135                xDelete<int>(gate.numdone);
    136136        }
    … …  
    154154                /*launch the thread manager with Krigingxt as a core: */
    155155                LaunchThread(v4t,(void*)&gate,num);
    156                 _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     156                _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    157157                xDelete<int>(gate.numdone);
    158158        }
    … …  
    177177                /*launch the thread manager with Krigingxt as a core: */
    178178                LaunchThread(Krigingxt,(void*)&gate,num);
    179                 _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     179                _printf_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    180180                xDelete<int>(gate.numdone);
    181181        }

  • issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp

    r15100 r15104  
    6868                /*partition loop across threads: */
    6969                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    70                         _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%" << "\n");
     70                        _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n");
    7171                        observations->InterpolationKriging(&predictions[idx],&error[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,variogram);
    7272                }
    73                 _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%" << "\n");
     73                _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n");
    7474
    7575#ifdef _HAVE_MPI_
    … …  
    8686                /*partition loop across threads: */
    8787                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    88                         _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%" << "\n");
     88                        _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n");
    8989                        observations->InterpolationV4(&predictions[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata);
    9090                }
    91                 _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%" << "\n");
     91                _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n");
    9292
    9393#ifdef _HAVE_MPI_
    … …  
    101101                /*partition loop across threads: */
    102102                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    103                         _printf0_("      interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%" << "\n");
     103                        _printf0_("      interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n");
    104104                        observations->InterpolationNearestNeighbor(&predictions[idx],x_interp[idx],y_interp[idx],radius);
    105105                }
    106                 _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%" << "\n");
     106                _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n");
    107107
    108108#ifdef _HAVE_MPI_
    … …  
    118118                /*partition loop across threads: */
    119119                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    120                         _printf0_("      interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%" << "\n");
     120                        _printf0_("      interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n");
    121121                        observations->InterpolationIDW(&predictions[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,power);
    122122                }
    123                 _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%" << "\n");
     123                _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n");
    124124
    125125#ifdef _HAVE_MPI_
    … …  
    142142
    143143        MPI_Barrier(MPI_COMM_WORLD); finish=MPI_Wtime( );
    144         _printf_("\n   " << setw(34) << left << "Observation fitering elapsed time: " << finish_init-start_init << " seconds  \n" << "\n");
    145         _printf_("   " << setw(34) << left << "Kriging prediction elapsed time: " << finish_core-start_core << " seconds  \n" << "\n");
    146         _printf_("\n   " << "Total elapsed time " << int((finish-start)/3600) << " hrs " << int(int(finish-start)%3600/60) << " min " << int(finish-start)%60 << " sec\n\n" << "\n");
     144        _printf_("\n   " << setw(34) << left << "Observation fitering elapsed time: " << finish_init-start_init << " seconds  \n\n");
     145        _printf_("   " << setw(34) << left << "Kriging prediction elapsed time: " << finish_core-start_core << " seconds  \n\n");
     146        _printf_("\n   " << "Total elapsed time " << int((finish-start)/3600) << " hrs " << int(int(finish-start)%3600/60) << " min " << int(finish-start)%60 << " sec\n\n\n");
    147147        return 1;
    148148#else

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

    r15100 r15104  
    1717
    1818        /*Display message*/
    19         if(VerboseModule()) _printf0_("   Merging solution vector from fset to gset" << "\n");
     19        if(VerboseModule()) _printf0_("   Merging solution vector from fset to gset\n");
    2020
    2121        /*first, get gsize, fsize and ssize: */

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

    r15100 r15104  
    4141        SetVerbosityLevel(verbose);
    4242
    43         if(VerboseMProcessor()) _printf0_("   starting model processor " << "\n");
     43        if(VerboseMProcessor()) _printf0_("   starting model processor \n");
    4444
    4545        for(i=0;i<nummodels;i++){
    … …  
    6868                CreateDataSets(&elements,&nodes,&vertices,&materials,&constraints,&loads,&parameters,iomodel,rootpath,solution_type,analysis_type,nummodels,i);
    6969        }
    70         if(VerboseMProcessor()) _printf0_("   done with model processor " << "\n");
     70        if(VerboseMProcessor()) _printf0_("   done with model processor \n");
    7171
    7272        /*Free resources:*/

  • issm/trunk-jpl/src/c/modules/PointCloudFindNeighborsx/PointCloudFindNeighborsxt.cpp

    r15100 r15104  
    7070        }
    7171        if (my_thread==0)
    72          _printf_("\r      loop progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%" << "\n");
     72         _printf_("\r      loop progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%\n");
    7373
    7474        /*Free ressources:*/

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

    r15100 r15104  
    5454
    5555  if (itm >= NPDMAX){
    56     _printf_("increase NPDMAX in massBalance.cpp" << "\n");
     56    _printf_("increase NPDMAX in massBalance.cpp\n");
    5757    exit (1);
    5858      }
    … …  
    8080  itm = reCast<int,IssmDouble>((PDCUT+2.*siglimc)/DT + 1.5);
    8181  if (itm >= NPDCMAX){
    82     _printf_("'increase NPDCMAX in p35com'" << "\n");
     82    _printf_("'increase NPDCMAX in p35com'\n");
    8383    exit (1);
    8484      }

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

    r15100 r15104  
    2222        bool        oldalloc  = false;
    2323
    24         if(VerboseModule()) _printf0_("   Dirichlet lifting applied to load vector" << "\n");
     24        if(VerboseModule()) _printf0_("   Dirichlet lifting applied to load vector\n");
    2525
    2626        Kfs->GetSize(&global_m,&global_n);

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

    r15100 r15104  
    1717        bool        oldalloc  = false;
    1818
    19         if(VerboseModule()) _printf0_("   Reduce vector from g to f set" << "\n");
     19        if(VerboseModule()) _printf0_("   Reduce vector from g to f set\n");
    2020
    2121        /*first figure out fsize: */

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

    r15100 r15104  
    1414
    1515        /*Display message*/
    16         if(VerboseModule()) _printf0_("   Resetting penalties" << "\n");
     16        if(VerboseModule()) _printf0_("   Resetting penalties\n");
    1717
    1818        /*recover parameters: */

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

    r15100 r15104  
    295295
    296296                        if (pnpart[i] > 0)
    297                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n" << "\n");
     297                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n\n");
    298298                        if (pnvert[i] > 1)
    299                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnvert[i] << " > 1 vertices.\n" << "\n");
     299                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnvert[i] << " > 1 vertices.\n\n");
    300300
    301301                        kpoint=new KML_Point();
    … …  
    464464                                else {
    465465                                        if (!kpoly) {
    466                                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\", part " << j << ", expected to be outer loop (cw).\n" << "\n");
     466                                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\", part " << j << ", expected to be outer loop (cw).\n\n");
    467467                                                continue;
    468468                                        }
    … …  
    512512
    513513                        if (pnpart[i] > 0)
    514                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n" << "\n");
     514                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n\n");
    515515
    516516/*  create a multigeometry to hold all the points  */
    … …  
    555555
    556556                else if (pstype[i] == SHPT_MULTIPATCH) {
    557                         _printf_("Warning -- Shape " << i << " of type \"" <<SHPTypeName( pstype[i] ) << "\" will be ignored.\n" << "\n");
     557                        _printf_("Warning -- Shape " << i << " of type \"" <<SHPTypeName( pstype[i] ) << "\" will be ignored.\n\n");
    558558                        continue;
    559559                }
    … …  
    562562
    563563                else {
    564                         _printf_(,"Warning -- Shape " << i << " of type \"" <<SHPTypeName( pstype[i] ) << "\" will be ignored.\n" << "\n");
     564                        _printf_(,"Warning -- Shape " << i << " of type \"" <<SHPTypeName( pstype[i] ) << "\" will be ignored.\n\n");
    565565                }
    566566        }
    … …  
    575575/*  open exp file  */
    576576
    577         _printf0_("Writing exp profiles to file." << "\n");
     577        _printf0_("Writing exp profiles to file.\n");
    578578        fid=fopen(filexp,"w");
    579579
    … …  
    609609        clock1=clock();
    610610        time1 =time(NULL);
    611         _printf_("Shp2Expx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0)  << " elapsed seconds.\n\n" << "\n");
     611        _printf_("Shp2Expx Module -- " <<((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " <<difftime(time1,time0)  << " elapsed seconds.\n\n\n");
    612612
    613613        return(iret);

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

    r15100 r15104  
    283283
    284284                        if (pnpart[i] > 0)
    285                                 _printf_("Warning -- Shape "<< i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n" << "\n");
     285                                _printf_("Warning -- Shape "<< i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n\n");
    286286                        if (pnvert[i] > 1)
    287                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 1 vertices.\n" << "\n");
     287                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 1 vertices.\n\n");
    288288
    289289                        kpoint=new KML_Point();
    … …  
    452452                                else {
    453453                                        if (!kpoly) {
    454                                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\", part " << j << ", expected to be outer loop (cw).\n" << "\n");
     454                                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\", part " << j << ", expected to be outer loop (cw).\n\n");
    455455                                                continue;
    456456                                        }
    … …  
    500500
    501501                        if (pnpart[i] > 0)
    502                                 _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n" << "\n");
     502                                _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" should not have " << pnpart[i] << " > 0 parts.\n\n");
    503503
    504504/*  create a multigeometry to hold all the points  */
    … …  
    543543
    544544                else if (pstype[i] == SHPT_MULTIPATCH) {
    545                         _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" will be ignored.\n" << "\n");
     545                        _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" will be ignored.\n\n");
    546546                        continue;
    547547                }
    … …  
    550550
    551551                else {
    552                         _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" will be ignored.\n" << "\n");
     552                        _printf_("Warning -- Shape " << i << " of type \"" << SHPTypeName( pstype[i] ) << "\" will be ignored.\n\n");
    553553                }
    554554        }
    … …  
    563563/*  write kml file  */
    564564
    565         _printf0_("Writing kml document to file." << "\n");
     565        _printf0_("Writing kml document to file.\n");
    566566        fid=fopen(filkml,"w");
    567567        fprintf(fid,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
    … …  
    592592        clock1=clock();
    593593        time1 =time(NULL);
    594         _printf_("Shp2Kmlx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n" << "\n");
     594        _printf_("Shp2Kmlx Module -- " << ((double)(clock1-clock0))/CLOCKS_PER_SEC << " CPU seconds; " << difftime(time1,time0) << " elapsed seconds.\n\n\n");
    595595
    596596        return(iret);

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

    r15100 r15104  
    2020        Vector<IssmDouble> *uf=NULL;
    2121
    22         if(VerboseModule()) _printf0_("   Solving matrix system" << "\n");
     22        if(VerboseModule()) _printf0_("   Solving matrix system\n");
    2323
    2424        /*Initialize solver: */

  • issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp

    r15100 r15104  
    55
    66void printarray(IssmPDouble* array,int lines,int cols){
    7         _printf_("" << "\n");
     7        _printf_("\n");
    88        for(int i=0;i<lines;i++){ 
    99                _printf_("   [ ");
    1010                for(int j=0;j<cols;j++) _printf_( " " << setw(11) << setprecision (5) << array[i*cols+j]);
    11                 _printf_(" ]" << "\n");
     11                _printf_(" ]\n");
    1212        } 
    13         _printf_("" << "\n");
     13        _printf_("\n");
    1414}
    1515void printsparsity(IssmPDouble* array,int lines,int cols){
    1616        int count;
    17         _printf_("" << "\n");
     17        _printf_("\n");
    1818        for(int i=0;i<lines;i++){ 
    1919                _printf_("   [ ");
    … …  
    2828                _printf_(" ] "<<i<<" => "<<count << "\n");
    2929        } 
    30         _printf_("" << "\n");
     30        _printf_("\n");
    3131}
    3232void printarray(int* array,int lines,int cols){
    33         _printf_("" << "\n");
     33        _printf_("\n");
    3434        for(int i=0;i<lines;i++){ 
    3535                _printf_("   [ ");
    3636                for(int j=0;j<cols;j++) _printf_( " " << setw(11) << setprecision (5) << array[i*cols+j]);
    37                 _printf_(" ]" << "\n");
     37                _printf_(" ]\n");
    3838        } 
    39         _printf_("" << "\n");
     39        _printf_("\n");
    4040}
    4141void printarray(bool* array,int lines,int cols){
    42         _printf_("" << "\n");
     42        _printf_("\n");
    4343        for(int i=0;i<lines;i++){ 
    4444                _printf_("   [ ");
    … …  
    4747                        else _printf_( " 0");
    4848                }
    49                 _printf_(" ]" << "\n");
     49                _printf_(" ]\n");
    5050        } 
    51         _printf_("" << "\n");
     51        _printf_("\n");
    5252}
    5353void printbinary(int n){

  • issm/trunk-jpl/src/c/shared/Exceptions/Exceptions.cpp

    r15100 r15104  
    4747                if(num_procs==1){
    4848                        _printf_("\n??? Error in ==> " << file_name.c_str() << ":" << file_line << "\n");
    49                         _printf_(function_name.c_str() << " error message: " << what() << "\n" << "\n");
     49                        _printf_(function_name.c_str() << " error message: " << what() << "\n\n");
    5050                }
    5151                else{
    5252                        _printf_("\n[" << my_rank<< "] ??? Error using ==> " << file_name.c_str() << ":" << file_line << "\n");
    53                         _printf_(  "[" << my_rank << "] " << function_name.c_str() << " error message: " << what() << "\n" << "\n");
     53                        _printf_(  "[" << my_rank << "] " << function_name.c_str() << " error message: " << what() << "\n\n");
    5454                }
    5555        }

  • issm/trunk-jpl/src/c/shared/Exceptions/exceptions.h

    r15100 r15104  
    6464        }\
    6565        catch (exception& e) {\
    66                 _printf_("Standard exception: " << e.what() << "\n" << "\n");\
     66                _printf_("Standard exception: " << e.what() << "\n\n");\
    6767                return 0;\
    6868        }\
    6969        catch(...){\
    70                 _printf_("An unexpected error occurred \n" << "\n");\
     70                _printf_("An unexpected error occurred \n\n");\
    7171                return 0;\
    7272        }

  • issm/trunk-jpl/src/c/shared/LatLong/Ll2xyx.cpp

    r15100 r15104  
    114114                *pdelta= 45;
    115115                *pslat = 70;
    116                 if(flag) _printf0_("Info: creating coordinates in polar stereographic (Std Latitude: 70N Meridian: 45)." << "\n");
     116                if(flag) _printf0_("Info: creating coordinates in polar stereographic (Std Latitude: 70N Meridian: 45).\n");
    117117        }
    118118        else if (sgn == -1) {
    119119                *pdelta= 0;
    120120                *pslat = 71;
    121                 if(flag) _printf0_("Info: creating coordinates in polar stereographic (Std Latitude: 71S Meridian: 0)." << "\n");
     121                if(flag) _printf0_("Info: creating coordinates in polar stereographic (Std Latitude: 71S Meridian: 0).\n");
    122122        }
    123123        else _error_("Sign should be either +1 or -1.\n");

  • issm/trunk-jpl/src/c/shared/LatLong/Xy2llx.cpp

    r15100 r15104  
    122122                *pdelta= 45;
    123123                *pslat = 70;
    124                 if(flag) _printf0_("Warning: expecting coordinates in polar stereographic (Std Latitude: 70N Meridian: 45)." << "\n");
     124                if(flag) _printf0_("Warning: expecting coordinates in polar stereographic (Std Latitude: 70N Meridian: 45).\n");
    125125        }
    126126        else if (sgn == -1) {
    127127                *pdelta= 0;
    128128                *pslat = 71;
    129                 if(flag) _printf0_("Warning: expecting coordinates in polar stereographic (Std Latitude: 71S Meridian: 0)." << "\n");
     129                if(flag) _printf0_("Warning: expecting coordinates in polar stereographic (Std Latitude: 71S Meridian: 0).\n");
    130130        }
    131131        else _error_("Sign should be either +1 or -1.\n");

  • issm/trunk-jpl/src/c/shared/Matrix/MatrixUtils.cpp

    r15100 r15104  
    205205                pindx[ipiv]++;
    206206
    207                 //              _printf_("pivot for i=" << i << ": irow=" << irow << ", icol=" << icol  << ", pindx[" << ipiv << "]=" << pindx[ipiv] << "\n" << "\n");
     207                //              _printf_("pivot for i=" << i << ": irow=" << irow << ", icol=" << icol  << ", pindx[" << ipiv << "]=" << pindx[ipiv] << "\n\n");
    208208
    209209                /*  switch rows to put pivot element on diagonal, noting that the
    … …  
    211211
    212212                if (irow != icol) {
    213                         //                      _printf_("row switch for i=" << i << ": irow=" << irow << ", icol=" << icol << "\n" << "\n");
     213                        //                      _printf_("row switch for i=" << i << ": irow=" << irow << ", icol=" << icol << "\n\n");
    214214
    215215                        ipt=irow*ndim;
    … …  
    292292                        icol=pivrc2[j];
    293293
    294                         //                      _printf_("column switch back for j=" << j << ": irow=" << irow << ", icol=" << icol << "\n" << "\n");
     294                        //                      _printf_("column switch back for j=" << j << ": irow=" << irow << ", icol=" << icol << "\n\n");
    295295
    296296                        ipt=0;

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

    r15100 r15104  
    5050        if (xIsNan<IssmDouble>(fxmin)) _error_("Function evaluation returned NaN");
    5151        cout<<setprecision(5);
    52         if(VerboseControl()) _printf0_("" << "\n");
    53         if(VerboseControl()) _printf0_("       Iteration         x           f(x)       Tolerance         Procedure" << "\n");
    54         if(VerboseControl()) _printf0_("" << "\n");
    55         if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<xmin<<"  "<<setw(12)<<fxmin<<"           N/A         boundary" << "\n");
     52        if(VerboseControl()) _printf0_("\n");
     53        if(VerboseControl()) _printf0_("       Iteration         x           f(x)       Tolerance         Procedure\n");
     54        if(VerboseControl()) _printf0_("\n");
     55        if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<xmin<<"  "<<setw(12)<<fxmin<<"           N/A         boundary\n");
    5656        fxmax = (*f)(xmax,optargs);
    5757        if (xIsNan<IssmDouble>(fxmax)) _error_("Function evaluation returned NaN");
    58         if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<xmax<<"  "<<setw(12)<<fxmax<<"           N/A         boundary" << "\n");
     58        if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<xmax<<"  "<<setw(12)<<fxmax<<"           N/A         boundary\n");
    5959
    6060        /*test if jump option activated and xmin==0*/
    … …  
    9393        /*4: print result*/
    9494        if(VerboseControl())
    95          _printf0_("         "<<setw(5)<<iter<<"    "<<setw(12)<<xbest<<"  "<<setw(12)<<fxbest<<"  "<<setw(12)<<pow(pow(xbest-xm,2),0.5)<<"         initial" << "\n");
     95         _printf0_("         "<<setw(5)<<iter<<"    "<<setw(12)<<xbest<<"  "<<setw(12)<<fxbest<<"  "<<setw(12)<<pow(pow(xbest-xm,2),0.5)<<"         initial\n");
    9696        if (!xIsNan<IssmDouble>(cm_jump) && (xmin==0) && ((fxbest/fxmin)<cm_jump)){
    9797                if(VerboseControl()) _printf0_("      optimization terminated: current x satisfies criteria 'cm_jump'=" << cm_jump << "\n");
    … …  
    198198                }
    199199                else if (iter>=maxiter){
    200                         if(VerboseControl()) _printf0_("      exiting: Maximum number of iterations has been exceeded  ('maxiter'=" << maxiter << ")" << "\n");
     200                        if(VerboseControl()) _printf0_("      exiting: Maximum number of iterations has been exceeded  ('maxiter'=" << maxiter << ")\n");
    201201                        loop=false;
    202202                }

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

    r15100 r15104  
    12001200        //      _printf0_("GaussLegendreTria - ngaus=" << *pngaus << "\n");
    12011201        //      for (i=0; i<*pngaus; i++)
    1202         //              _printf0_("i=" << i << ": l1gaus=" << (*pl1 )[i] << ",l2gaus=" << (*pl2 )[i] << ",l3gaus=" << (*pl3 )[i]<< ",wgt=" << (*pwgt)[i]<< "\n" << "\n");
     1202        //              _printf0_("i=" << i << ": l1gaus=" << (*pl1 )[i] << ",l2gaus=" << (*pl2 )[i] << ",l3gaus=" << (*pl3 )[i]<< ",wgt=" << (*pwgt)[i]<< "\n\n");
    12031203
    12041204        return;

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

    r15100 r15104  
    4747        if (xIsNan<IssmDouble>(fx1)) _error_("Function evaluation returned NaN");
    4848        cout<<setprecision(5);
    49         if(VerboseControl()) _printf0_("" << "\n");
    50         if(VerboseControl()) _printf0_("       Iteration         x           f(x)       Tolerance" << "\n");
    51         if(VerboseControl()) _printf0_("" << "\n");
    52         if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<x1<<"  "<<setw(12)<<fx1<<"           N/A" << "\n");
     49        if(VerboseControl()) _printf0_("\n");
     50        if(VerboseControl()) _printf0_("       Iteration         x           f(x)       Tolerance\n");
     51        if(VerboseControl()) _printf0_("\n");
     52        if(VerboseControl()) _printf0_("           N/A    "<<setw(12)<<x1<<"  "<<setw(12)<<fx1<<"           N/A\n");
    5353
    5454        //update tolerances
    … …  
    7171                }
    7272                else if (iter>=maxiter){
    73                         if(VerboseControl()) _printf0_("      " << "exiting: Maximum number of iterations has been exceeded  - increase 'maxiter'" << "\n");
     73                        if(VerboseControl()) _printf0_("      " << "exiting: Maximum number of iterations has been exceeded  - increase 'maxiter'\n");
    7474                        loop=false;
    7575                }

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

    r15100 r15104  
    3030        IssmDouble yts;
    3131
    32         if(VerboseModule()) _printf0_("   checking convergence" << "\n");
     32        if(VerboseModule()) _printf0_("   checking convergence\n");
    3333
    3434        /*If uf is NULL in input, f-set is nil, model is fully constrained, therefore converged from
    … …  
    7373        res=nKUoldF/nF;
    7474        if (xIsNan<IssmDouble>(res)){
    75                 _printf0_("norm nf = " << nF << "f and norm kuold = " << nKUoldF << "f" << "\n");
     75                _printf0_("norm nf = " << nF << "f and norm kuold = " << nKUoldF << "f\n");
    7676                _error_("mechanical equilibrium convergence criterion is NaN!");
    7777        }
    … …  
    8383        //print
    8484        if(res<eps_res){
    85                 if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   mechanical equilibrium convergence criterion"<<res*100<< " < "<<eps_res*100<<" %" << "\n");
     85                if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   mechanical equilibrium convergence criterion"<<res*100<< " < "<<eps_res*100<<" %\n");
    8686                converged=true;
    8787        }
    8888        else{
    89                 if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   mechanical equilibrium convergence criterion"<<res*100<<" > "<<eps_res*100<<" %" << "\n");
     89                if(VerboseConvergence()) _printf0_(setw(50)<<left<<"   mechanical equilibrium convergence criterion"<<res*100<<" > "<<eps_res*100<<" %\n");
    9090                converged=false;
    9191        }
    … …  
    106106                if (!xIsNan<IssmDouble>(eps_rel)){
    107107                        if((ndu/nu)<eps_rel){
    108                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " < " << eps_rel*100 << " %" << "\n");
     108                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " < " << eps_rel*100 << " %\n");
    109109                        }
    110110                        else{
    111                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " > " << eps_rel*100 << " %" << "\n");
     111                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " > " << eps_rel*100 << " %\n");
    112112                                converged=false;
    113113                        }
    114114                }
    115                 else _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " %" << "\n");
     115                else _printf0_(setw(50) << left << "   Convergence criterion: norm(du)/norm(u)" << ndu/nu*100 << " %\n");
    116116
    117117        }
    … …  
    131131                if (!xIsNan<IssmDouble>(eps_abs)){
    132132                        if ((nduinf*yts)<eps_abs){
    133                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " < " << eps_abs << " m/yr" << "\n");
     133                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " < " << eps_abs << " m/yr\n");
    134134                        }
    135135                        else{
    136                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " > " << eps_abs << " m/yr" << "\n");
     136                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " > " << eps_abs << " m/yr\n");
    137137                                converged=false;
    138138                        }
    139139                }
    140                 else  _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " m/yr" << "\n");
     140                else  _printf0_(setw(50) << left << "   Convergence criterion: max(du)" << nduinf*yts << " m/yr\n");
    141141
    142142        }

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

    r15100 r15104  
    164164                        if (!xIsNan<IssmDouble>(eps_hyd)){
    165165                                if((ndu_sed/nu_sed)<eps_hyd){
    166                                         if(VerboseConvergence()) _printf0_(setw(50) << left << "   Sediment Convergence criterion:" << ndu_sed/nu_sed*100 << " < " << eps_hyd*100 << " %" << "\n");
     166                                        if(VerboseConvergence()) _printf0_(setw(50) << left << "   Sediment Convergence criterion:" << ndu_sed/nu_sed*100 << " < " << eps_hyd*100 << " %\n");
    167167                                        hydroconverged=true;
    168168                                }
    169169                                else{
    170                                         if(VerboseConvergence()) _printf0_(setw(50) << left << "   Sediment Convergence criterion:" << ndu_sed/nu_sed*100 << " > " << eps_hyd*100 << " %" << "\n");
     170                                        if(VerboseConvergence()) _printf0_(setw(50) << left << "   Sediment Convergence criterion:" << ndu_sed/nu_sed*100 << " > " << eps_hyd*100 << " %\n");
    171171                                        hydroconverged=false;
    172172                                }
    173173                                if(isefficientlayer){
    174174                                        if((ndu_epl/nu_epl)<eps_hyd){
    175                                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   EPL Convergence criterion:" << ndu_epl/nu_epl*100 << " < " << eps_hyd*100 << " %" << "\n");
     175                                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   EPL Convergence criterion:" << ndu_epl/nu_epl*100 << " < " << eps_hyd*100 << " %\n");
    176176                                        }
    177177                                        else{
    178                                                 if(VerboseConvergence()) _printf0_(setw(50) << left << "   EPL Convergence criterion:" << ndu_epl/nu_epl*100 << " > " << eps_hyd*100 << " %" << "\n");
     178                                                if(VerboseConvergence()) _printf0_(setw(50) << left << "   EPL Convergence criterion:" << ndu_epl/nu_epl*100 << " > " << eps_hyd*100 << " %\n");
    179179                                                hydroconverged=false;
    180180                                        }
    181181                                }
    182182                        }
    183                         else _printf0_(setw(50) << left << "   Convergence criterion:" << ndu_sed/nu_sed*100 << " %" << "\n");
     183                        else _printf0_(setw(50) << left << "   Convergence criterion:" << ndu_sed/nu_sed*100 << " %\n");
    184184                        if (hydrocount>=hydro_maxiter){
    185185                                _error_("   maximum number for hydrological global iterations (" << hydro_maxiter << ") exceeded");

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

    r15100 r15104  
    9393                }
    9494                if(count>=max_nonlinear_iterations){
    95                         _printf0_("   maximum number of Newton iterations (" << max_nonlinear_iterations << ") exceeded" << "\n");
     95                        _printf0_("   maximum number of Newton iterations (" << max_nonlinear_iterations << ") exceeded\n");
    9696                        bool max_iteration_state=true;
    9797                        int tempStep=1;

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

    r15100 r15104  
    9090                }
    9191                if(count>=max_nonlinear_iterations){
    92                         _printf0_("   maximum number of nonlinear iterations (" << max_nonlinear_iterations << ") exceeded" << "\n");
     92                        _printf0_("   maximum number of nonlinear iterations (" << max_nonlinear_iterations << ") exceeded\n");
    9393                        converged=true;
    9494                        InputUpdateFromConstantx( femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,converged,ConvergedEnum);

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

    r15100 r15104  
    8989                if(converged==true)break;
    9090                if(count>=max_nonlinear_iterations){
    91                         _printf0_("   maximum number of iterations (" << max_nonlinear_iterations << ") exceeded" << "\n");
     91                        _printf0_("   maximum number of iterations (" << max_nonlinear_iterations << ") exceeded\n");
    9292                        break;
    9393                }

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

    r15100 r15104  
    6363                        if (count>=thermal_maxiter){
    6464                                converged=true;
    65                                 _printf0_("   maximum number of iterations (" << thermal_maxiter << ") exceeded" << "\n");
     65                                _printf0_("   maximum number of iterations (" << thermal_maxiter << ") exceeded\n");
    6666                        }
    6767                }

  • issm/trunk-jpl/src/c/toolkits/issm/Bucket.h

    r15100 r15104  
    8888                /*object virtual functions definitions:*/
    8989                void    Echo(){ /*{{{*/
    90                         _printf_("Bucket echo (cpu #: "<<IssmComm::GetRank()<<")" << "\n");
     90                        _printf_("Bucket echo (cpu #: "<<IssmComm::GetRank()<<")\n");
    9191                        _printf_("bucket type: " << type << "\n");
    9292                        _printf_("num rows: "<<this->m<<" num cols: "<<this->n << "\n");
    … …  
    9595                        int i,j;
    9696
    97                         _printf_("Bucket echo (cpu #: "<<IssmComm::GetRank()<<")" << "\n");
     97                        _printf_("Bucket echo (cpu #: "<<IssmComm::GetRank()<<")\n");
    9898                        _printf_("bucket type: " << type << "\n");
    9999                        _printf_("num rows: "<<this->m<<" num cols: "<<this->n << "\n");
    100100                        if(type==MATRIX_BUCKET){
    101101                                for (i=0;i<this->m;i++){
    102                                         _printf_("row "<<this->idxm[i]<<", column indices: " << "\n");
     102                                        _printf_("row "<<this->idxm[i]<<", column indices: \n");
    103103                                        for (j=0;j<this->n;j++){
    104104                                                _printf_(" "<<this->idxn[j] << "\n");
    105105                                        }
    106                                         _printf_("values: " << "\n");
     106                                        _printf_("values: \n");
    107107                                        for (j=0;j<this->n;j++){
    108108                                                _printf_(" "<<this->values[m*i+j] << "\n");

  • issm/trunk-jpl/src/c/toolkits/issm/IssmDenseMat.h

    r15100 r15104  
    120120                                        _printf_(this->matrix[N*i+j] << " ");
    121121                                }
    122                                 _printf_("" << "\n");
     122                                _printf_("\n");
    123123                        }
    124124                }

  • issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h

    r15100 r15104  
    148148                                                _printf_("row " << j << "\n");
    149149                                                for (k=0;k<this->N;k++){
    150                                                         _printf_(this->matrix[j*this->N+k] << " " << "\n");
     150                                                        _printf_(this->matrix[j*this->N+k] << " \n");
    151151                                                }
    152                                                 _printf_("" << "\n");
     152                                                _printf_("\n");
    153153                                        }
    154154                                }

  • issm/trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h

    r15100 r15104  
    136136                                                _printf_("row " << j << " " << this->vector[j] << "\n");
    137137                                        }
    138                                         _printf_("" << "\n");
     138                                        _printf_("\n");
    139139                                }
    140140                                MPI_Barrier(IssmComm::GetComm());

  • issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp

    r15100 r15104  
    9393        MatInfo info;
    9494        MatGetInfo(this->matrix,MAT_GLOBAL_SUM,&info);
    95         _printf0_("=========================== Stiffness matrix allocation info ===========================" << "\n");
    96         _printf0_("" << "\n");
     95        _printf0_("=========================== Stiffness matrix allocation info ===========================\n");
     96        _printf0_("\n");
    9797   _printf0_(" Block size  : "<<info.block_size << "\n");
    9898        _printf0_(" nz_allocated: "<<info.nz_allocated << "\n");
    9999        _printf0_(" nz_used     : "<<info.nz_used << "\n");
    100         _printf0_(" nz_unneeded : "<<info.nz_unneeded<<" ("<<double(info.nz_unneeded)/double(info.nz_allocated)*100.<<"%)" << "\n");
    101         _printf0_("" << "\n");
    102         _printf0_("========================================================================================" << "\n");
     100        _printf0_(" nz_unneeded : "<<info.nz_unneeded<<" ("<<double(info.nz_unneeded)/double(info.nz_allocated)*100.<<"%)\n");
     101        _printf0_("\n");
     102        _printf0_("========================================================================================\n");
    103103}
    104104/*}}}*/
Note: See TracChangeset for help on using the changeset viewer.