47 #include <Xyce_config.h>
52 #include <N_UTL_fwd.h>
56 #include "N_LAS_Vector.h"
57 #include "N_LAS_Matrix.h"
58 #include "Epetra_MapColoring.h"
59 #include "N_ERH_ErrorMgr.h"
75 (IO::InitialConditionsData::NodeNamePairMap & op_in,
76 const Teuchos::RCP <Epetra_MapColoring>& color_map,
77 Linear::Vector* cloneVector
82 color_map_ = color_map;
83 tmp_vector_ptr_ =
new Linear::Vector(*cloneVector);
108 (
const Linear::Vector * solution, Linear::Vector * residual_vector)
111 dout() <<
"Inside AugmentLinSysIC::augmentResidual:" << std::endl;
114 IO::InitialConditionsData::NodeNamePairMap::iterator op_i = op_.begin();
115 IO::InitialConditionsData::NodeNamePairMap::iterator op_end = op_.end();
116 for ( ; op_i != op_end ; ++op_i)
118 int row = (*op_i).second.first;
121 if ( (*color_map_)[row] == 0)
123 (*residual_vector)[row] = 0.0;
140 dout() <<
"Inside AugmentLinSysIC::augmentJacobian:" << std::endl;
143 std::vector<int> col;
144 std::vector<double> val;
145 IO::InitialConditionsData::NodeNamePairMap::iterator op_i =
op_.begin();
146 IO::InitialConditionsData::NodeNamePairMap::iterator op_end =
op_.end();
150 for ( ; op_i != op_end ; ++op_i)
152 int row = (*op_i).second.first;
158 rowLen = jacobian->getLocalRowLength(row);
160 col.resize(rowLen,0);
161 val.resize(rowLen,0.0);
162 jacobian->getLocalRowCopy(row, rowLen, numEntries, &val[0], &col[0]);
165 for (
int i=0;i<val.size();++i) val[i] = 0.0;
167 jacobian->putLocalRow(row, rowLen, &val[0], &col[0]);
170 (*tmp_vector_ptr_)[row] = 1.0;
void augmentResidual(const Xyce::Linear::Vector *solution, Xyce::Linear::Vector *residual_vector)
Augments the Residual.
Xyce::IO::InitialConditionsData::NodeNamePairMap & op_
map of specified variables
Pure virtual class to augment a linear system.
AugmentLinSysIC(Xyce::IO::InitialConditionsData::NodeNamePairMap &op_in, const Teuchos::RCP< Epetra_MapColoring > &color_map, Xyce::Linear::Vector *cloneVector)
Ctor.
Xyce::Linear::Vector * tmp_vector_ptr_
Temporary vector used to store diagonal.
Teuchos::RCP< Epetra_MapColoring > color_map_
Color 0 are the voltage unknowns.
void augmentJacobian(Xyce::Linear::Matrix *jacobian)
Augments the Jacobian.