diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2021-11-08 23:19:54 +0100 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2021-11-08 23:19:54 +0100 |
| commit | 54416465850a9280a121127e0b5301d97ad13a61 (patch) | |
| tree | dd6cdf2fdb63e6dd268db45dfb7300a079041ebc /stokes.hpp | |
| parent | 6b42ce01cd289567a6109d831fe86a0667d18f37 (diff) | |
| download | code-54416465850a9280a121127e0b5301d97ad13a61.tar.gz code-54416465850a9280a121127e0b5301d97ad13a61.tar.bz2 code-54416465850a9280a121127e0b5301d97ad13a61.zip | |
Removed some unused files.
Diffstat (limited to 'stokes.hpp')
| -rw-r--r-- | stokes.hpp | 2 |
1 files changed, 1 insertions, 1 deletions
@@ -336,7 +336,7 @@ public: std::pair<Vector<Scalar>, Matrix<Scalar>> gsGradHess(const Tensor<Scalar, p>& J, Real t) const { Vector<Scalar> z = f(t); auto [H, dH, ddH] = hamGradHess(J, z); - Tensor<Scalar, 3> dddH = ((p - 2) * (p - 1) * p / (Real)factorial(p)) * J; + Tensor<Scalar, 3> dddH = ((p - 2) * (p - 1) * p / (Real)factorial(p)) * J; // BUG IF P != 3 Vector<Scalar> ż = zDot(z, dH); Tensor<Scalar, 1> żT = Eigen::TensorMap<Tensor<const Scalar, 1>>(ż.data(), {z.size()}); Vector<Scalar> dz = df(t); |