eschnett/cactusamrex - Change 3YL5CNXXPPYWRQRXMW36VFFJP5S4RD45ISQ3TWAUBX44HIM7ECTQC

Weyl: Use tensor norm when checking symmetry/antisymmetry

Created by Erik Schnetter on April 3, 2020

3YL5CNXXPPYWRQRXMW36VFFJP5S4RD45ISQ3TWAUBX44HIM7ECTQC

Dependencies

In channels

master

main

Change contents

Replacement in Weyl/src/tensor.hxx at line 334 [3.15983]

B:BD[3.25449] → [3.25449:25624]

    const auto is_approx{[](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <=
             1.0e-12 * (1 + norm1<T>()(fgood) + norm1<T>()(fother));

[3.25449]

[3.25624]

    const auto scale = norm1<T>()(elts[0]) + norm1<T>()(elts[1]) +
                       norm1<T>()(elts[2]) + norm1<T>()(elts[3]) +
                       norm1<T>()(elts[4]) + norm1<T>()(elts[5]);
    const auto is_approx{[scale](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <= 1.0e-12 * (1 + scale);

Replacement in Weyl/src/tensor.hxx at line 804 [3.15983]

B:BD[3.41222] → [3.41222:41397]

    const auto is_approx{[](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <=
             1.0e-12 * (1 + norm1<T>()(fgood) + norm1<T>()(fother));

[3.41222]

[3.41397]

    const auto scale = norm1<T>()(elts[0]) + norm1<T>()(elts[1]) +
                       norm1<T>()(elts[2]) + norm1<T>()(elts[3]) +
                       norm1<T>()(elts[4]) + norm1<T>()(elts[5]) +
                       norm1<T>()(elts[6]) + norm1<T>()(elts[7]) +
                       norm1<T>()(elts[8]) + norm1<T>()(elts[9]);
    const auto is_approx{[scale](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <= 1.0e-12 * (1 + scale);

Replacement in Weyl/src/tensor.hxx at line 1183 [3.15983]
B:BD[3.53701] → [3.53701:53723]
```
    // Check symmetry
```
[3.53701]
[2.7999]
```
    // Check antisymmetry
```

Replacement in Weyl/src/tensor.hxx at line 1194 [3.15983]

∅:D[2.8110] → [3.53885:54060]

B:BD[3.53885] → [3.53885:54060]

    const auto is_approx{[](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <=
             1.0e-12 * (1 + norm1<T>()(fgood) + norm1<T>()(fother));

[2.8110]

[3.54060]

    const auto scale = norm1<T>()(elts[0]) + norm1<T>()(elts[1]) +
                       norm1<T>()(elts[2]) + norm1<T>()(elts[3]) +
                       norm1<T>()(elts[4]) + norm1<T>()(elts[5]);
    const auto is_approx{[scale](const T &fgood, const T &fother) {
      return norm1<T>()(fother - fgood) <= 1.0e-12 * (1 + scale);