From 6738ff87b9562b31d6f7cb8247872de5d7255461 Mon Sep 17 00:00:00 2001
From: Louis Denis <louis.denis@student.uliege.be>
Date: Tue, 22 Feb 2022 15:34:48 +0100
Subject: [PATCH] le meilleur projet c'est le notre
---
examples/gmsh_api/my_code.cpp | 14 ++++++++++----
1 file changed, 10 insertions(+), 4 deletions(-)
diff --git a/examples/gmsh_api/my_code.cpp b/examples/gmsh_api/my_code.cpp
index 1cf0149..2f49a60 100644
--- a/examples/gmsh_api/my_code.cpp
+++ b/examples/gmsh_api/my_code.cpp
@@ -192,16 +192,22 @@ int main(int argc, char **argv)
// convert first jacobian to Eigen format
Eigen::Map<Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> >
jacob(&jacobians[0], 3, 3);
+ // passage en 2D xD
+ Eigen::Matrix<double, 2, 2> jacob2D;
+ jacob2D(0,0) = jacob(0,0);
+ jacob2D(1,0) = jacob(1,0);
+ jacob2D(0,1) = jacob(0,1);
+ jacob2D(1,1) = jacob(1,1);
std::cout << "\tfirst jacobian (Eigen):\n"
- << jacob.format(fmt) << '\n';
-
+ << jacob2D.format(fmt) << '\n';
+
// compute the inverse with Eigen
- Eigen::Matrix<double, 3, 3> jacobinv = jacob.inverse();
+ Eigen::Matrix<double, 2, 2> jacobinv = jacob2D.inverse();
std::cout << "\tinverse of the first jacobian (Eigen):\n"
<< jacobinv.format(fmt) << "\n";
// compute the transpose of the inverse with Eigen
- Eigen::Matrix<double, 3,3> jacobinvtrans = jacobinv.transpose();
+ Eigen::Matrix<double, 2, 2> jacobinvtrans = jacobinv.transpose();
std::cout << "\ttranspose of the inverse of the first jacobian (Eigen):\n"
<< jacobinvtrans.format(fmt) << "\n";
--
GitLab