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()
164 template <
class RhsValueType>
165 OPM_HOST_DEVICE
static Evaluation createVariable(
const RhsValueType& value,
int varPos)
172 template <
class RhsValueType>
173 OPM_HOST_DEVICE
static Evaluation createVariable(
int nVars,
const RhsValueType& value,
int varPos)
176 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
177 " with 10 derivatives");
184 template <
class RhsValueType>
185 OPM_HOST_DEVICE
static Evaluation createVariable(
const Evaluation&,
const RhsValueType& value,
int varPos)
195 template <
class RhsValueType>
196 OPM_HOST_DEVICE
static Evaluation createConstant(
int nVars,
const RhsValueType& value)
199 throw std::logic_error(
"This statically-sized evaluation can only represent objects"
200 " with 10 derivatives");
206 template <
class RhsValueType>
207 OPM_HOST_DEVICE
static Evaluation createConstant(
const RhsValueType& value)
214 template <
class RhsValueType>
221 OPM_HOST_DEVICE
void copyDerivatives(
const Evaluation& other)
223 assert(
size() == other.size());
225 data_[1] = other.data_[1];
226 data_[2] = other.data_[2];
227 data_[3] = other.data_[3];
228 data_[4] = other.data_[4];
229 data_[5] = other.data_[5];
230 data_[6] = other.data_[6];
231 data_[7] = other.data_[7];
232 data_[8] = other.data_[8];
233 data_[9] = other.data_[9];
234 data_[10] = other.data_[10];
241 assert(
size() == other.size());
243 data_[0] += other.data_[0];
244 data_[1] += other.data_[1];
245 data_[2] += other.data_[2];
246 data_[3] += other.data_[3];
247 data_[4] += other.data_[4];
248 data_[5] += other.data_[5];
249 data_[6] += other.data_[6];
250 data_[7] += other.data_[7];
251 data_[8] += other.data_[8];
252 data_[9] += other.data_[9];
253 data_[10] += other.data_[10];
259 template <
class RhsValueType>
260 OPM_HOST_DEVICE
Evaluation& operator+=(
const RhsValueType& other)
271 assert(
size() == other.size());
273 data_[0] -= other.data_[0];
274 data_[1] -= other.data_[1];
275 data_[2] -= other.data_[2];
276 data_[3] -= other.data_[3];
277 data_[4] -= other.data_[4];
278 data_[5] -= other.data_[5];
279 data_[6] -= other.data_[6];
280 data_[7] -= other.data_[7];
281 data_[8] -= other.data_[8];
282 data_[9] -= other.data_[9];
283 data_[10] -= other.data_[10];
289 template <
class RhsValueType>
290 OPM_HOST_DEVICE
Evaluation& operator-=(
const RhsValueType& other)
301 assert(
size() == other.size());
312 data_[1] = data_[1] * v + other.data_[1] * u;
313 data_[2] = data_[2] * v + other.data_[2] * u;
314 data_[3] = data_[3] * v + other.data_[3] * u;
315 data_[4] = data_[4] * v + other.data_[4] * u;
316 data_[5] = data_[5] * v + other.data_[5] * u;
317 data_[6] = data_[6] * v + other.data_[6] * u;
318 data_[7] = data_[7] * v + other.data_[7] * u;
319 data_[8] = data_[8] * v + other.data_[8] * u;
320 data_[9] = data_[9] * v + other.data_[9] * u;
321 data_[10] = data_[10] * v + other.data_[10] * u;
327 template <
class RhsValueType>
328 OPM_HOST_DEVICE
Evaluation& operator*=(
const RhsValueType& other)
348 assert(
size() == other.size());
354 data_[1] = (v*data_[1] - u*other.data_[1])/(v*v);
355 data_[2] = (v*data_[2] - u*other.data_[2])/(v*v);
356 data_[3] = (v*data_[3] - u*other.data_[3])/(v*v);
357 data_[4] = (v*data_[4] - u*other.data_[4])/(v*v);
358 data_[5] = (v*data_[5] - u*other.data_[5])/(v*v);
359 data_[6] = (v*data_[6] - u*other.data_[6])/(v*v);
360 data_[7] = (v*data_[7] - u*other.data_[7])/(v*v);
361 data_[8] = (v*data_[8] - u*other.data_[8])/(v*v);
362 data_[9] = (v*data_[9] - u*other.data_[9])/(v*v);
363 data_[10] = (v*data_[10] - u*other.data_[10])/(v*v);
370 template <
class RhsValueType>
371 OPM_HOST_DEVICE
Evaluation& operator/=(
const RhsValueType& other)
393 assert(
size() == other.size());
403 template <
class RhsValueType>
404 OPM_HOST_DEVICE
Evaluation operator+(
const RhsValueType& other)
const
416 assert(
size() == other.size());
426 template <
class RhsValueType>
427 OPM_HOST_DEVICE
Evaluation operator-(
const RhsValueType& other)
const
442 result.data_[0] = - data_[0];
443 result.data_[1] = - data_[1];
444 result.data_[2] = - data_[2];
445 result.data_[3] = - data_[3];
446 result.data_[4] = - data_[4];
447 result.data_[5] = - data_[5];
448 result.data_[6] = - data_[6];
449 result.data_[7] = - data_[7];
450 result.data_[8] = - data_[8];
451 result.data_[9] = - data_[9];
452 result.data_[10] = - data_[10];
459 assert(
size() == other.size());
468 template <
class RhsValueType>
469 OPM_HOST_DEVICE
Evaluation operator*(
const RhsValueType& other)
const
480 assert(
size() == other.size());
489 template <
class RhsValueType>
490 OPM_HOST_DEVICE
Evaluation operator/(
const RhsValueType& other)
const
499 template <
class RhsValueType>
500 OPM_HOST_DEVICE
Evaluation& operator=(
const RhsValueType& other)
511 template <
class RhsValueType>
512 OPM_HOST_DEVICE
bool operator==(
const RhsValueType& other)
const
513 {
return value() == other; }
515 OPM_HOST_DEVICE
bool operator==(
const Evaluation& other)
const
517 assert(
size() == other.size());
519 for (
int idx = 0; idx <
length_(); ++idx) {
520 if (data_[idx] != other.data_[idx]) {
527 OPM_HOST_DEVICE
bool operator!=(
const Evaluation& other)
const
528 {
return !operator==(other); }
530 template <
class RhsValueType>
531 OPM_HOST_DEVICE
bool operator!=(
const RhsValueType& other)
const
532 {
return !operator==(other); }
534 template <
class RhsValueType>
535 OPM_HOST_DEVICE
bool operator>(RhsValueType other)
const
536 {
return value() > other; }
538 OPM_HOST_DEVICE
bool operator>(
const Evaluation& other)
const
540 assert(
size() == other.size());
542 return value() > other.value();
545 template <
class RhsValueType>
546 OPM_HOST_DEVICE
bool operator<(RhsValueType other)
const
547 {
return value() < other; }
549 OPM_HOST_DEVICE
bool operator<(
const Evaluation& other)
const
551 assert(
size() == other.size());
553 return value() < other.value();
556 template <
class RhsValueType>
557 OPM_HOST_DEVICE
bool operator>=(RhsValueType other)
const
558 {
return value() >= other; }
560 OPM_HOST_DEVICE
bool operator>=(
const Evaluation& other)
const
562 assert(
size() == other.size());
564 return value() >= other.value();
567 template <
class RhsValueType>
568 OPM_HOST_DEVICE
bool operator<=(RhsValueType other)
const
569 {
return value() <= other; }
571 OPM_HOST_DEVICE
bool operator<=(
const Evaluation& other)
const
573 assert(
size() == other.size());
575 return value() <= other.value();
579 OPM_HOST_DEVICE
const ValueType& value()
const
583 template <
class RhsValueType>
584 OPM_HOST_DEVICE
constexpr void setValue(
const RhsValueType& val)
588 OPM_HOST_DEVICE
const ValueType& derivative(
int varIdx)
const
590 assert(0 <= varIdx && varIdx <
size());
592 return data_[
dstart_() + varIdx];
596 OPM_HOST_DEVICE
void setDerivative(
int varIdx,
const ValueType& derVal)
598 assert(0 <= varIdx && varIdx <
size());
600 data_[
dstart_() + varIdx] = derVal;
603 template<
class Serializer>
604 OPM_HOST_DEVICE
void serializeOp(Serializer& serializer)
610 std::array<ValueT, 11> data_;