26 : id_( 0 ), n_( 1, 0.0 ), nInit_( 1, 0.0 ),
27 diffConst_( 1.0e-12 ), motorConst_( 0.0 )
34 assert( voxel <
nInit_.size() );
40 assert( voxel <
nInit_.size() );
46 assert( voxel <
n_.size() );
52 assert( voxel <
n_.size() );
58 assert( voxel <
n_.size() );
59 return prev_[ voxel ];
69 assert( vec.size() ==
n_.size() );
74 vector< double >::const_iterator q )
76 assert( start + num <=
n_.size() );
77 vector< double >::iterator p =
n_.begin() + start;
78 for (
unsigned int i = 0; i < num; ++i )
109 nInit_.resize( num, 0.0 );
110 n_.resize( num, 0.0 );
129 const vector< double >& diagVal )
131 if ( ops.size() > 0 )
133 assert( diagVal.size() ==
n_.size() );
146 if (
ops_.size() == 0 )
return;
148 i =
ops_.begin(); i !=
ops_.end(); ++i )
149 n_[i->c_] -=
n_[i->b_] * i->a_;
152 vector< double >::iterator iy =
n_.begin();
153 for ( vector< double >::const_iterator
160 assert(
n_.size() ==
nInit_.size() );
vector< double > prev_
Number of molecules of pool in each voxel.
void setOps(const vector< Triplet< double > > &ops_, const vector< double > &diagVal_)
Assigns prev_ = n_.
const vector< double > & getNvec() const
Used by parent solver to manipulate 'n'.
void setN(unsigned int vox, double value)
unsigned int getNumVoxels() const
unsigned int getId() const
unsigned int id_
Assign operations.
void setNvec(const vector< double > &n)
Used by parent solver to manipulate 'n'.
void setId(unsigned int id)
vector< double > n_
Integer conversion of Id of pool handled.
double getDiffConst() const
double getNinit(unsigned int vox) const
vector< Triplet< double > > ops_
Motor const, ie, transport rate.
void setNinit(unsigned int vox, double value)
vector< double > diagVal_
double diffConst_
Boundary condition: Initial 'n'.
void setMotorConst(double value)
double motorConst_
Diffusion const, assumed uniform.
void setNumVoxels(unsigned int num)
double getMotorConst() const
void setDiffConst(double value)
vector< double > nInit_
molecules of pool on previous timestep
double getN(unsigned int vox) const
double getPrev(unsigned int vox) const