Ignore:
Timestamp:
03/12/12 09:35:33 (13 years ago)
Author:
Eric.Larour
Message:

Large commit, to strip out Petsc from the code. Created new Matrix and Vector
wrapper objects to the Petsc Mat and Vec objects. These wrappers make it possible
to hook up a different package to handle matrix and vector operations.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/objects/Node.cpp

    r11001 r11679  
    594594/*}}}*/
    595595/*FUNCTION Node::CreateNodalConstraints{{{1*/
    596 void  Node::CreateNodalConstraints(Vec ys){
     596void  Node::CreateNodalConstraints(Vector* ys){
    597597
    598598        int i;
     
    613613               
    614614                /*Add values into constraint vector: */
    615                 VecSetValues(ys,this->indexing.ssize,this->indexing.sdoflist,values,INSERT_VALUES);
     615                ys->SetValues(this->indexing.ssize,this->indexing.sdoflist,values,INSERT_VALUES);
    616616        }
    617617
     
    875875/*}}}*/
    876876/*FUNCTION Node::VecMerge {{{1*/
    877 void   Node::VecMerge(Vec ug, double* vector_serial,int setenum){
     877void   Node::VecMerge(Vector* ug, double* vector_serial,int setenum){
    878878
    879879        double* values=NULL;
     
    897897
    898898                        /*Add values into ug: */
    899                         VecSetValues(ug,this->indexing.fsize,indices,(const double*)values,INSERT_VALUES);
     899                        ug->SetValues(this->indexing.fsize,indices,values,INSERT_VALUES);
    900900                }
    901901        }
     
    915915
    916916                        /*Add values into ug: */
    917                         VecSetValues(ug,this->indexing.ssize,indices,(const double*)values,INSERT_VALUES);
     917                        ug->SetValues(this->indexing.ssize,indices,values,INSERT_VALUES);
    918918                }
    919919        }
     
    926926/*}}}*/
    927927/*FUNCTION Node::VecReduce {{{1*/
    928 void   Node::VecReduce(Vec vector, double* ug_serial,int setenum){
     928void   Node::VecReduce(Vector* vector, double* ug_serial,int setenum){
    929929
    930930        double* values=NULL;
     
    945945
    946946                        /*Add values into ug: */
    947                         VecSetValues(vector,this->indexing.fsize,this->indexing.fdoflist,(const double*)values,INSERT_VALUES);
     947                        vector->SetValues(this->indexing.fsize,this->indexing.fdoflist,values,INSERT_VALUES);
    948948                }
    949949        }
     
    961961
    962962                        /*Add values into ug: */
    963                         VecSetValues(vector,this->indexing.ssize,this->indexing.sdoflist,(const double*)values,INSERT_VALUES);
     963                        vector->SetValues(this->indexing.ssize,this->indexing.sdoflist,values,INSERT_VALUES);
    964964                }
    965965        }
Note: See TracChangeset for help on using the changeset viewer.