60 OPM_HOST_DEVICE
constexpr int size()
const
65 OPM_HOST_DEVICE
constexpr int length_()
const
66 {
return size() + 1; }
73 OPM_HOST_DEVICE
constexpr int dstart_()
const
76 OPM_HOST_DEVICE
constexpr int dend_()
const
84 for (
const auto& v: data_)
85 Valgrind::CheckDefined(v);
102 template <
class RhsValueType>
103 OPM_HOST_DEVICE
Evaluation(
const RhsValueType& c)
115 template <
class RhsValueType>
116 OPM_HOST_DEVICE
Evaluation(
const RhsValueType& c,
int varPos)
119 assert(0 <= varPos && varPos <
size());
124 data_[varPos +
dstart_()] = 1.0;
130 OPM_HOST_DEVICE
void clearDerivatives()
162 template <
class RhsValueType>
163 OPM_HOST_DEVICE
static Evaluation createVariable(
const RhsValueType& value,
int varPos)
170 template <
class RhsValueType>
171 OPM_HOST_DEVICE
static Evaluation createVariable(
int nVars,
const RhsValueType& value,
int varPos)
174 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
175 " with 8 derivatives");
182 template <
class RhsValueType>
183 OPM_HOST_DEVICE
static Evaluation createVariable(
const Evaluation&,
const RhsValueType& value,
int varPos)
193 template <
class RhsValueType>
194 OPM_HOST_DEVICE
static Evaluation createConstant(
int nVars,
const RhsValueType& value)
197 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
198 " with 8 derivatives");
204 template <
class RhsValueType>
205 OPM_HOST_DEVICE
static Evaluation createConstant(
const RhsValueType& value)
212 template <
class RhsValueType>
219 OPM_HOST_DEVICE
void copyDerivatives(
const Evaluation& other)
221 assert(
size() == other.size());
223 data_[1] = other.data_[1];
224 data_[2] = other.data_[2];
225 data_[3] = other.data_[3];
226 data_[4] = other.data_[4];
227 data_[5] = other.data_[5];
228 data_[6] = other.data_[6];
229 data_[7] = other.data_[7];
230 data_[8] = other.data_[8];
237 assert(
size() == other.size());
239 data_[0] += other.data_[0];
240 data_[1] += other.data_[1];
241 data_[2] += other.data_[2];
242 data_[3] += other.data_[3];
243 data_[4] += other.data_[4];
244 data_[5] += other.data_[5];
245 data_[6] += other.data_[6];
246 data_[7] += other.data_[7];
247 data_[8] += other.data_[8];
253 template <
class RhsValueType>
254 OPM_HOST_DEVICE
Evaluation& operator+=(
const RhsValueType& other)
265 assert(
size() == other.size());
267 data_[0] -= other.data_[0];
268 data_[1] -= other.data_[1];
269 data_[2] -= other.data_[2];
270 data_[3] -= other.data_[3];
271 data_[4] -= other.data_[4];
272 data_[5] -= other.data_[5];
273 data_[6] -= other.data_[6];
274 data_[7] -= other.data_[7];
275 data_[8] -= other.data_[8];
281 template <
class RhsValueType>
282 OPM_HOST_DEVICE
Evaluation& operator-=(
const RhsValueType& other)
293 assert(
size() == other.size());
304 data_[1] = data_[1] * v + other.data_[1] * u;
305 data_[2] = data_[2] * v + other.data_[2] * u;
306 data_[3] = data_[3] * v + other.data_[3] * u;
307 data_[4] = data_[4] * v + other.data_[4] * u;
308 data_[5] = data_[5] * v + other.data_[5] * u;
309 data_[6] = data_[6] * v + other.data_[6] * u;
310 data_[7] = data_[7] * v + other.data_[7] * u;
311 data_[8] = data_[8] * v + other.data_[8] * u;
317 template <
class RhsValueType>
318 OPM_HOST_DEVICE
Evaluation& operator*=(
const RhsValueType& other)
336 assert(
size() == other.size());
342 data_[1] = (v*data_[1] - u*other.data_[1])/(v*v);
343 data_[2] = (v*data_[2] - u*other.data_[2])/(v*v);
344 data_[3] = (v*data_[3] - u*other.data_[3])/(v*v);
345 data_[4] = (v*data_[4] - u*other.data_[4])/(v*v);
346 data_[5] = (v*data_[5] - u*other.data_[5])/(v*v);
347 data_[6] = (v*data_[6] - u*other.data_[6])/(v*v);
348 data_[7] = (v*data_[7] - u*other.data_[7])/(v*v);
349 data_[8] = (v*data_[8] - u*other.data_[8])/(v*v);
356 template <
class RhsValueType>
357 OPM_HOST_DEVICE
Evaluation& operator/=(
const RhsValueType& other)
377 assert(
size() == other.size());
387 template <
class RhsValueType>
388 OPM_HOST_DEVICE
Evaluation operator+(
const RhsValueType& other)
const
400 assert(
size() == other.size());
410 template <
class RhsValueType>
411 OPM_HOST_DEVICE
Evaluation operator-(
const RhsValueType& other)
const
426 result.data_[0] = - data_[0];
427 result.data_[1] = - data_[1];
428 result.data_[2] = - data_[2];
429 result.data_[3] = - data_[3];
430 result.data_[4] = - data_[4];
431 result.data_[5] = - data_[5];
432 result.data_[6] = - data_[6];
433 result.data_[7] = - data_[7];
434 result.data_[8] = - data_[8];
441 assert(
size() == other.size());
450 template <
class RhsValueType>
451 OPM_HOST_DEVICE
Evaluation operator*(
const RhsValueType& other)
const
462 assert(
size() == other.size());
471 template <
class RhsValueType>
472 OPM_HOST_DEVICE
Evaluation operator/(
const RhsValueType& other)
const
481 template <
class RhsValueType>
482 OPM_HOST_DEVICE
Evaluation& operator=(
const RhsValueType& other)
493 template <
class RhsValueType>
494 OPM_HOST_DEVICE
bool operator==(
const RhsValueType& other)
const
495 {
return value() == other; }
497 OPM_HOST_DEVICE
bool operator==(
const Evaluation& other)
const
499 assert(
size() == other.size());
501 for (
int idx = 0; idx <
length_(); ++idx) {
502 if (data_[idx] != other.data_[idx]) {
509 OPM_HOST_DEVICE
bool operator!=(
const Evaluation& other)
const
510 {
return !operator==(other); }
512 template <
class RhsValueType>
513 OPM_HOST_DEVICE
bool operator!=(
const RhsValueType& other)
const
514 {
return !operator==(other); }
516 template <
class RhsValueType>
517 OPM_HOST_DEVICE
bool operator>(RhsValueType other)
const
518 {
return value() > other; }
520 OPM_HOST_DEVICE
bool operator>(
const Evaluation& other)
const
522 assert(
size() == other.size());
524 return value() > other.value();
527 template <
class RhsValueType>
528 OPM_HOST_DEVICE
bool operator<(RhsValueType other)
const
529 {
return value() < other; }
531 OPM_HOST_DEVICE
bool operator<(
const Evaluation& other)
const
533 assert(
size() == other.size());
535 return value() < other.value();
538 template <
class RhsValueType>
539 OPM_HOST_DEVICE
bool operator>=(RhsValueType other)
const
540 {
return value() >= other; }
542 OPM_HOST_DEVICE
bool operator>=(
const Evaluation& other)
const
544 assert(
size() == other.size());
546 return value() >= other.value();
549 template <
class RhsValueType>
550 OPM_HOST_DEVICE
bool operator<=(RhsValueType other)
const
551 {
return value() <= other; }
553 OPM_HOST_DEVICE
bool operator<=(
const Evaluation& other)
const
555 assert(
size() == other.size());
557 return value() <= other.value();
561 OPM_HOST_DEVICE
const ValueType& value()
const
565 template <
class RhsValueType>
566 OPM_HOST_DEVICE
void setValue(
const RhsValueType& val)
570 OPM_HOST_DEVICE
const ValueType& derivative(
int varIdx)
const
572 assert(0 <= varIdx && varIdx <
size());
574 return data_[
dstart_() + varIdx];
578 OPM_HOST_DEVICE
void setDerivative(
int varIdx,
const ValueType& derVal)
580 assert(0 <= varIdx && varIdx <
size());
582 data_[
dstart_() + varIdx] = derVal;
585 template<
class Serializer>
586 OPM_HOST_DEVICE
void serializeOp(Serializer& serializer)
592 std::array<ValueT, 9> data_;