30 template<
int PolySize>
38 rhoPolynomial_(rhoPoly)
44 template<
int PolySize>
51 rhoPolynomial_(ip.rhoPolynomial_)
55 template<
int PolySize>
59 const icoPolynomial<PolySize>& ip
63 rhoPolynomial_(ip.rhoPolynomial_)
67 template<
int PolySize>
78 template<
int PolySize>
88 template<
int PolySize>
91 return rhoPolynomial_.evaluate(T)/this->W();
95 template<
int PolySize>
102 template<
int PolySize>
111 template<
int PolySize>
119 rhoPolynomial_ = ip.rhoPolynomial_;
125 template<
int PolySize>
126 inline void Foam::icoPolynomial<PolySize>::operator+=
131 scalar molr1 = this->nMoles();
135 molr1 /= this->nMoles();
136 scalar molr2 = ip.nMoles()/this->nMoles();
138 rhoPolynomial_ = molr1*rhoPolynomial_ + molr2*ip.rhoPolynomial_;
142 template<
int PolySize>
143 inline void Foam::icoPolynomial<PolySize>::operator-=
148 scalar molr1 = this->nMoles();
152 molr1 /= this->nMoles();
153 scalar molr2 = ip.nMoles()/this->nMoles();
155 rhoPolynomial_ = molr1*rhoPolynomial_ - molr2*ip.rhoPolynomial_;
159 template<
int PolySize>
168 template<
int PolySize>
175 scalar nMoles = ip1.
nMoles() + ip2.nMoles();
176 scalar molr1 = ip1.
nMoles()/nMoles;
177 scalar molr2 = ip2.nMoles()/nMoles;
181 static_cast<const specie&
>(ip1)
182 + static_cast<const specie&>(ip2),
183 molr1*ip1.rhoPolynomial_ + molr2*ip2.rhoPolynomial_
188 template<
int PolySize>
195 scalar nMoles = ip1.
nMoles() + ip2.nMoles();
196 scalar molr1 = ip1.
nMoles()/nMoles;
197 scalar molr2 = ip2.nMoles()/nMoles;
201 static_cast<const specie&
>(ip1)
202 - static_cast<const specie&>(ip2),
203 molr1*ip1.rhoPolynomial_ - molr2*ip2.rhoPolynomial_
208 template<
int PolySize>
217 s*
static_cast<const specie&
>(ip),
223 template<
int PolySize>