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()
163 template <
class RhsValueType>
164 OPM_HOST_DEVICE
static Evaluation createVariable(
const RhsValueType& value,
int varPos)
171 template <
class RhsValueType>
172 OPM_HOST_DEVICE
static Evaluation createVariable(
int nVars,
const RhsValueType& value,
int varPos)
175 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
176 " with 9 derivatives");
183 template <
class RhsValueType>
184 OPM_HOST_DEVICE
static Evaluation createVariable(
const Evaluation&,
const RhsValueType& value,
int varPos)
194 template <
class RhsValueType>
195 OPM_HOST_DEVICE
static Evaluation createConstant(
int nVars,
const RhsValueType& value)
198 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
199 " with 9 derivatives");
205 template <
class RhsValueType>
206 OPM_HOST_DEVICE
static Evaluation createConstant(
const RhsValueType& value)
213 template <
class RhsValueType>
220 OPM_HOST_DEVICE
void copyDerivatives(
const Evaluation& other)
222 assert(
size() == other.size());
224 data_[1] = other.data_[1];
225 data_[2] = other.data_[2];
226 data_[3] = other.data_[3];
227 data_[4] = other.data_[4];
228 data_[5] = other.data_[5];
229 data_[6] = other.data_[6];
230 data_[7] = other.data_[7];
231 data_[8] = other.data_[8];
232 data_[9] = other.data_[9];
239 assert(
size() == other.size());
241 data_[0] += other.data_[0];
242 data_[1] += other.data_[1];
243 data_[2] += other.data_[2];
244 data_[3] += other.data_[3];
245 data_[4] += other.data_[4];
246 data_[5] += other.data_[5];
247 data_[6] += other.data_[6];
248 data_[7] += other.data_[7];
249 data_[8] += other.data_[8];
250 data_[9] += other.data_[9];
256 template <
class RhsValueType>
257 OPM_HOST_DEVICE
Evaluation& operator+=(
const RhsValueType& other)
268 assert(
size() == other.size());
270 data_[0] -= other.data_[0];
271 data_[1] -= other.data_[1];
272 data_[2] -= other.data_[2];
273 data_[3] -= other.data_[3];
274 data_[4] -= other.data_[4];
275 data_[5] -= other.data_[5];
276 data_[6] -= other.data_[6];
277 data_[7] -= other.data_[7];
278 data_[8] -= other.data_[8];
279 data_[9] -= other.data_[9];
285 template <
class RhsValueType>
286 OPM_HOST_DEVICE
Evaluation& operator-=(
const RhsValueType& other)
297 assert(
size() == other.size());
308 data_[1] = data_[1] * v + other.data_[1] * u;
309 data_[2] = data_[2] * v + other.data_[2] * u;
310 data_[3] = data_[3] * v + other.data_[3] * u;
311 data_[4] = data_[4] * v + other.data_[4] * u;
312 data_[5] = data_[5] * v + other.data_[5] * u;
313 data_[6] = data_[6] * v + other.data_[6] * u;
314 data_[7] = data_[7] * v + other.data_[7] * u;
315 data_[8] = data_[8] * v + other.data_[8] * u;
316 data_[9] = data_[9] * v + other.data_[9] * u;
322 template <
class RhsValueType>
323 OPM_HOST_DEVICE
Evaluation& operator*=(
const RhsValueType& other)
342 assert(
size() == other.size());
348 data_[1] = (v*data_[1] - u*other.data_[1])/(v*v);
349 data_[2] = (v*data_[2] - u*other.data_[2])/(v*v);
350 data_[3] = (v*data_[3] - u*other.data_[3])/(v*v);
351 data_[4] = (v*data_[4] - u*other.data_[4])/(v*v);
352 data_[5] = (v*data_[5] - u*other.data_[5])/(v*v);
353 data_[6] = (v*data_[6] - u*other.data_[6])/(v*v);
354 data_[7] = (v*data_[7] - u*other.data_[7])/(v*v);
355 data_[8] = (v*data_[8] - u*other.data_[8])/(v*v);
356 data_[9] = (v*data_[9] - u*other.data_[9])/(v*v);
363 template <
class RhsValueType>
364 OPM_HOST_DEVICE
Evaluation& operator/=(
const RhsValueType& other)
385 assert(
size() == other.size());
395 template <
class RhsValueType>
396 OPM_HOST_DEVICE
Evaluation operator+(
const RhsValueType& other)
const
408 assert(
size() == other.size());
418 template <
class RhsValueType>
419 OPM_HOST_DEVICE
Evaluation operator-(
const RhsValueType& other)
const
434 result.data_[0] = - data_[0];
435 result.data_[1] = - data_[1];
436 result.data_[2] = - data_[2];
437 result.data_[3] = - data_[3];
438 result.data_[4] = - data_[4];
439 result.data_[5] = - data_[5];
440 result.data_[6] = - data_[6];
441 result.data_[7] = - data_[7];
442 result.data_[8] = - data_[8];
443 result.data_[9] = - data_[9];
450 assert(
size() == other.size());
459 template <
class RhsValueType>
460 OPM_HOST_DEVICE
Evaluation operator*(
const RhsValueType& other)
const
471 assert(
size() == other.size());
480 template <
class RhsValueType>
481 OPM_HOST_DEVICE
Evaluation operator/(
const RhsValueType& other)
const
490 template <
class RhsValueType>
491 OPM_HOST_DEVICE
Evaluation& operator=(
const RhsValueType& other)
502 template <
class RhsValueType>
503 OPM_HOST_DEVICE
bool operator==(
const RhsValueType& other)
const
504 {
return value() == other; }
506 OPM_HOST_DEVICE
bool operator==(
const Evaluation& other)
const
508 assert(
size() == other.size());
510 for (
int idx = 0; idx <
length_(); ++idx) {
511 if (data_[idx] != other.data_[idx]) {
518 OPM_HOST_DEVICE
bool operator!=(
const Evaluation& other)
const
519 {
return !operator==(other); }
521 template <
class RhsValueType>
522 OPM_HOST_DEVICE
bool operator!=(
const RhsValueType& other)
const
523 {
return !operator==(other); }
525 template <
class RhsValueType>
526 OPM_HOST_DEVICE
bool operator>(RhsValueType other)
const
527 {
return value() > other; }
529 OPM_HOST_DEVICE
bool operator>(
const Evaluation& other)
const
531 assert(
size() == other.size());
533 return value() > other.value();
536 template <
class RhsValueType>
537 OPM_HOST_DEVICE
bool operator<(RhsValueType other)
const
538 {
return value() < other; }
540 OPM_HOST_DEVICE
bool operator<(
const Evaluation& other)
const
542 assert(
size() == other.size());
544 return value() < other.value();
547 template <
class RhsValueType>
548 OPM_HOST_DEVICE
bool operator>=(RhsValueType other)
const
549 {
return value() >= other; }
551 OPM_HOST_DEVICE
bool operator>=(
const Evaluation& other)
const
553 assert(
size() == other.size());
555 return value() >= other.value();
558 template <
class RhsValueType>
559 OPM_HOST_DEVICE
bool operator<=(RhsValueType other)
const
560 {
return value() <= other; }
562 OPM_HOST_DEVICE
bool operator<=(
const Evaluation& other)
const
564 assert(
size() == other.size());
566 return value() <= other.value();
570 OPM_HOST_DEVICE
const ValueType& value()
const
574 template <
class RhsValueType>
575 OPM_HOST_DEVICE
void setValue(
const RhsValueType& val)
579 OPM_HOST_DEVICE
const ValueType& derivative(
int varIdx)
const
581 assert(0 <= varIdx && varIdx <
size());
583 return data_[
dstart_() + varIdx];
587 OPM_HOST_DEVICE
void setDerivative(
int varIdx,
const ValueType& derVal)
589 assert(0 <= varIdx && varIdx <
size());
591 data_[
dstart_() + varIdx] = derVal;
594 template<
class Serializer>
595 OPM_HOST_DEVICE
void serializeOp(Serializer& serializer)
601 std::array<ValueT, 10> data_;