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
constexpr Evaluation(
const RhsValueType& c): data_{}
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
constexpr void clearDerivatives()
166 template <
class RhsValueType>
167 OPM_HOST_DEVICE
static Evaluation createVariable(
const RhsValueType& value,
int varPos)
174 template <
class RhsValueType>
175 OPM_HOST_DEVICE
static Evaluation createVariable(
int nVars,
const RhsValueType& value,
int varPos)
178 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
179 " with 12 derivatives");
186 template <
class RhsValueType>
187 OPM_HOST_DEVICE
static Evaluation createVariable(
const Evaluation&,
const RhsValueType& value,
int varPos)
197 template <
class RhsValueType>
198 OPM_HOST_DEVICE
static Evaluation createConstant(
int nVars,
const RhsValueType& value)
201 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
202 " with 12 derivatives");
208 template <
class RhsValueType>
209 OPM_HOST_DEVICE
static Evaluation createConstant(
const RhsValueType& value)
216 template <
class RhsValueType>
223 OPM_HOST_DEVICE
void copyDerivatives(
const Evaluation& other)
225 assert(
size() == other.size());
227 data_[1] = other.data_[1];
228 data_[2] = other.data_[2];
229 data_[3] = other.data_[3];
230 data_[4] = other.data_[4];
231 data_[5] = other.data_[5];
232 data_[6] = other.data_[6];
233 data_[7] = other.data_[7];
234 data_[8] = other.data_[8];
235 data_[9] = other.data_[9];
236 data_[10] = other.data_[10];
237 data_[11] = other.data_[11];
238 data_[12] = other.data_[12];
245 assert(
size() == other.size());
247 data_[0] += other.data_[0];
248 data_[1] += other.data_[1];
249 data_[2] += other.data_[2];
250 data_[3] += other.data_[3];
251 data_[4] += other.data_[4];
252 data_[5] += other.data_[5];
253 data_[6] += other.data_[6];
254 data_[7] += other.data_[7];
255 data_[8] += other.data_[8];
256 data_[9] += other.data_[9];
257 data_[10] += other.data_[10];
258 data_[11] += other.data_[11];
259 data_[12] += other.data_[12];
265 template <
class RhsValueType>
266 OPM_HOST_DEVICE
Evaluation& operator+=(
const RhsValueType& other)
277 assert(
size() == other.size());
279 data_[0] -= other.data_[0];
280 data_[1] -= other.data_[1];
281 data_[2] -= other.data_[2];
282 data_[3] -= other.data_[3];
283 data_[4] -= other.data_[4];
284 data_[5] -= other.data_[5];
285 data_[6] -= other.data_[6];
286 data_[7] -= other.data_[7];
287 data_[8] -= other.data_[8];
288 data_[9] -= other.data_[9];
289 data_[10] -= other.data_[10];
290 data_[11] -= other.data_[11];
291 data_[12] -= other.data_[12];
297 template <
class RhsValueType>
298 OPM_HOST_DEVICE
Evaluation& operator-=(
const RhsValueType& other)
309 assert(
size() == other.size());
320 data_[1] = data_[1] * v + other.data_[1] * u;
321 data_[2] = data_[2] * v + other.data_[2] * u;
322 data_[3] = data_[3] * v + other.data_[3] * u;
323 data_[4] = data_[4] * v + other.data_[4] * u;
324 data_[5] = data_[5] * v + other.data_[5] * u;
325 data_[6] = data_[6] * v + other.data_[6] * u;
326 data_[7] = data_[7] * v + other.data_[7] * u;
327 data_[8] = data_[8] * v + other.data_[8] * u;
328 data_[9] = data_[9] * v + other.data_[9] * u;
329 data_[10] = data_[10] * v + other.data_[10] * u;
330 data_[11] = data_[11] * v + other.data_[11] * u;
331 data_[12] = data_[12] * v + other.data_[12] * u;
337 template <
class RhsValueType>
338 OPM_HOST_DEVICE
Evaluation& operator*=(
const RhsValueType& other)
360 assert(
size() == other.size());
366 data_[1] = (v*data_[1] - u*other.data_[1])/(v*v);
367 data_[2] = (v*data_[2] - u*other.data_[2])/(v*v);
368 data_[3] = (v*data_[3] - u*other.data_[3])/(v*v);
369 data_[4] = (v*data_[4] - u*other.data_[4])/(v*v);
370 data_[5] = (v*data_[5] - u*other.data_[5])/(v*v);
371 data_[6] = (v*data_[6] - u*other.data_[6])/(v*v);
372 data_[7] = (v*data_[7] - u*other.data_[7])/(v*v);
373 data_[8] = (v*data_[8] - u*other.data_[8])/(v*v);
374 data_[9] = (v*data_[9] - u*other.data_[9])/(v*v);
375 data_[10] = (v*data_[10] - u*other.data_[10])/(v*v);
376 data_[11] = (v*data_[11] - u*other.data_[11])/(v*v);
377 data_[12] = (v*data_[12] - u*other.data_[12])/(v*v);
384 template <
class RhsValueType>
385 OPM_HOST_DEVICE
Evaluation& operator/=(
const RhsValueType& other)
409 assert(
size() == other.size());
419 template <
class RhsValueType>
420 OPM_HOST_DEVICE
Evaluation operator+(
const RhsValueType& other)
const
432 assert(
size() == other.size());
442 template <
class RhsValueType>
443 OPM_HOST_DEVICE
Evaluation operator-(
const RhsValueType& other)
const
458 result.data_[0] = - data_[0];
459 result.data_[1] = - data_[1];
460 result.data_[2] = - data_[2];
461 result.data_[3] = - data_[3];
462 result.data_[4] = - data_[4];
463 result.data_[5] = - data_[5];
464 result.data_[6] = - data_[6];
465 result.data_[7] = - data_[7];
466 result.data_[8] = - data_[8];
467 result.data_[9] = - data_[9];
468 result.data_[10] = - data_[10];
469 result.data_[11] = - data_[11];
470 result.data_[12] = - data_[12];
477 assert(
size() == other.size());
486 template <
class RhsValueType>
487 OPM_HOST_DEVICE
Evaluation operator*(
const RhsValueType& other)
const
498 assert(
size() == other.size());
507 template <
class RhsValueType>
508 OPM_HOST_DEVICE
Evaluation operator/(
const RhsValueType& other)
const
517 template <
class RhsValueType>
518 OPM_HOST_DEVICE
Evaluation& operator=(
const RhsValueType& other)
529 template <
class RhsValueType>
530 OPM_HOST_DEVICE
bool operator==(
const RhsValueType& other)
const
531 {
return value() == other; }
533 OPM_HOST_DEVICE
bool operator==(
const Evaluation& other)
const
535 assert(
size() == other.size());
537 for (
int idx = 0; idx <
length_(); ++idx) {
538 if (data_[idx] != other.data_[idx]) {
545 OPM_HOST_DEVICE
bool operator!=(
const Evaluation& other)
const
546 {
return !operator==(other); }
548 template <
class RhsValueType>
549 OPM_HOST_DEVICE
bool operator!=(
const RhsValueType& other)
const
550 {
return !operator==(other); }
552 template <
class RhsValueType>
553 OPM_HOST_DEVICE
bool operator>(RhsValueType other)
const
554 {
return value() > other; }
556 OPM_HOST_DEVICE
bool operator>(
const Evaluation& other)
const
558 assert(
size() == other.size());
560 return value() > other.value();
563 template <
class RhsValueType>
564 OPM_HOST_DEVICE
bool operator<(RhsValueType other)
const
565 {
return value() < other; }
567 OPM_HOST_DEVICE
bool operator<(
const Evaluation& other)
const
569 assert(
size() == other.size());
571 return value() < other.value();
574 template <
class RhsValueType>
575 OPM_HOST_DEVICE
bool operator>=(RhsValueType other)
const
576 {
return value() >= other; }
578 OPM_HOST_DEVICE
bool operator>=(
const Evaluation& other)
const
580 assert(
size() == other.size());
582 return value() >= other.value();
585 template <
class RhsValueType>
586 OPM_HOST_DEVICE
bool operator<=(RhsValueType other)
const
587 {
return value() <= other; }
589 OPM_HOST_DEVICE
bool operator<=(
const Evaluation& other)
const
591 assert(
size() == other.size());
593 return value() <= other.value();
597 OPM_HOST_DEVICE
const ValueType& value()
const
601 template <
class RhsValueType>
602 OPM_HOST_DEVICE
constexpr void setValue(
const RhsValueType& val)
606 OPM_HOST_DEVICE
const ValueType& derivative(
int varIdx)
const
608 assert(0 <= varIdx && varIdx <
size());
610 return data_[
dstart_() + varIdx];
614 OPM_HOST_DEVICE
void setDerivative(
int varIdx,
const ValueType& derVal)
616 assert(0 <= varIdx && varIdx <
size());
618 data_[
dstart_() + varIdx] = derVal;
621 template<
class Serializer>
622 OPM_HOST_DEVICE
void serializeOp(Serializer& serializer)
628 std::array<ValueT, 13> data_;