Changeset 7191 for issm/trunk


Ignore:
Timestamp:
01/26/11 16:20:51 (14 years ago)
Author:
Mathieu Morlighem
Message:

Fixed reduce load, look at Kfs size instead of what it is pointing to

Location:
issm/trunk/src/c
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r6580 r7191  
    1717        Vec         Kfsy_s = NULL;
    1818        int         Kfsm,Kfsn;
     19        int         global_m,global_n;
    1920        PetscScalar a;
    2021        bool        fromlocalsize = true;
    … …  
    2223
    2324        _printf_(VerboseModule(),"   Dirichlet lifting applied to load vector\n");
    24        
    25         if(pf && Kfs){
     25
     26        MatGetSize(Kfs,&global_m,&global_n);
     27        if(pf && global_m*global_n){
    2628
    2729                /*pf = pf - Kfs * y_s;*/

  • issm/trunk/src/c/toolkits/petsc/patches/MatMultPatch.cpp

    r464 r7191  
    2424        int m,n;
    2525        Vec X_rel=NULL;
     26
     27        _assert_(A); _assert_(X);
    2628
    2729        if (MatMultCompatible(A,X)){
Note: See TracChangeset for help on using the changeset viewer.