Changeset 7639


Ignore:
Timestamp:
03/15/11 09:00:52 (14 years ago)
Author:
Eric.Larour
Message:

added handy routines, like SetValue and SetDiag

Location:
issm/trunk/src/c/objects/Numerics
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Numerics/ElementMatrix.cpp

    r6412 r7639  
    460460}
    461461/*}}}*/
     462/*FUNCTION ElementMatrix::SetDiag{{{1*/
     463void ElementMatrix::SetDiag(double scalar){
     464
     465        int i;
     466
     467        if(this->nrows!=this->ncols)_error_("need square matrix in input!");
     468
     469        for(i=0;i<this->nrows;i++){
     470                this->values[this->ncols*i+i]=scalar;
     471        }
     472}
     473/*}}}*/
  • issm/trunk/src/c/objects/Numerics/ElementMatrix.h

    r6029 r7639  
    6464                void Transpose(void);
    6565                void Init(ElementMatrix* Ke);
     66                void SetDiag(double scalar);
    6667                /*}}}*/
    6768};
  • issm/trunk/src/c/objects/Numerics/ElementVector.cpp

    r7391 r7639  
    273273}
    274274/*}}}*/
     275/*FUNCTION ElementVector::SetValue{{{1*/
     276void ElementVector::SetValue(double scalar){
     277
     278        int i;
     279
     280        for(i=0;i<this->nrows;i++)this->values[i]=scalar;
     281
     282}
     283/*}}}*/
  • issm/trunk/src/c/objects/Numerics/ElementVector.h

    r7391 r7639  
    4545                void Echo(void);
    4646                void Init(ElementVector* pe);
     47                void SetValue(double scalar);
    4748                /*}}}*/
    4849};
Note: See TracChangeset for help on using the changeset viewer.