Changeset 61


Ignore:
Timestamp:
04/27/09 10:48:29 (16 years ago)
Author:
Eric.Larour
Message:

Wrong solver_string instead of solverstring

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/parallel/GradJCompute.cpp

    r46 r61  
    2020        /*intermediary: */
    2121        int analysis_type;
    22         char* solver_string=NULL;
     22        char* solverstring=NULL;
    2323        char* control_type=NULL;
    2424       
     
    4242        /*some parameters:*/
    4343        femmodel->parameters->FindParam((void*)&analysis_type,"analysis_type");
    44         femmodel->parameters->FindParam((void*)&solver_string,"solver_string");
     44        femmodel->parameters->FindParam((void*)&solverstring,"solverstring");
    4545        femmodel->parameters->FindParam((void*)&control_type,"control_type");
    4646
    47         _printf_("%s\n","   Recover solution for this stiffness and right hand side: \n");
     47        _printf_("%s\n","      recover solution for this stiffness and right hand side:");
    4848        diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,inputs,femmodel);
    4949        VecToMPISerial(&u_g_double,u_g); VecFree(&u_g);
    5050
    51         _printf_("%s\n","   Buid Du, difference between observed velocity and model velocity:\n");
     51        _printf_("%s\n","      buid Du, difference between observed velocity and model velocity:");
    5252        Dux( &du_g, femmodel->elements,femmodel->nodes,femmodel->loads,femmodel->materials,u_g_double,u_g_obs, inputs,analysis_type);
    5353
    54         _printf_("%s\n","   Reduce adjoint load from g-set to f-set:\n");
     54        _printf_("%s\n","      reduce adjoint load from g-set to f-set:");
    5555        Reduceloadfromgtofx(&du_f, du_g, femmodel->Gmn, K_fs0, femmodel->ys0, femmodel->nodesets);
    5656        VecFree(&du_g);MatFree(&K_fs0);
    5757
    58         _printf_("%s\n","   Solve for adjoint vector: \n");
    59         Solverx(&lambda_f, K_ff0, du_f, NULL, solver_string);
     58        _printf_("%s\n","      solve for adjoint vector:");
     59        Solverx(&lambda_f, K_ff0, du_f, NULL, solverstring);
    6060        VecFree(&du_f); MatFree(&K_ff0);
    6161       
    62         _printf_("%s\n","   Merge back to g set:\n");
     62        _printf_("%s\n","      merge back to g set:");
    6363        Mergesolutionfromftogx(&lambda_g, lambda_f,femmodel->Gmn,femmodel->ys0,femmodel->nodesets);
    6464        VecFree(&lambda_f);
    6565
    66         _printf_("%s\n","   Compute gradJ:\n");
     66        _printf_("%s\n","      compute gradJ:");
    6767        VecToMPISerial(&lambda_g_double,lambda_g);VecFree(&lambda_g);
    6868
     
    7171       
    7272        /*Free ressources:*/
    73         xfree((void**)&solver_string);
     73        xfree((void**)&solverstring);
    7474        xfree((void**)&control_type);
    7575        xfree((void**)&u_g_double);
Note: See TracChangeset for help on using the changeset viewer.