Commit 641803cb authored by p7k's avatar p7k

Remove fixed size Eigen::Vector3d and Eigen::Vector4d from tests due to...

Remove fixed size Eigen::Vector3d and Eigen::Vector4d from tests due to possible SSE alignment issues on Titan
parent 0969b17c
#!/usr/bin/env bash
module unload PrgEnv-intel PrgEnv-gnu PrgEnv-pgi PrgEnv-cray
module load PrgEnv-gnu
module swap gcc gcc/7.2.0
module load cmake3
\ No newline at end of file
......@@ -489,7 +489,8 @@ void calculate_torque_speed_curve(Eigen::MatrixXd &torque, std::array<Eigen::Mat
}
Eigen::MatrixXd XY = Eigen::MatrixXd::Zero(9,9);
Eigen::Vector3d xy; xy << -1.0, 0.0, 1.0;
Eigen::VectorXd xy(3);
xy << -1.0, 0.0, 1.0;
Eigen::VectorXd ztorque = Eigen::VectorXd::Zero(9);
Eigen::VectorXd zspeed = Eigen::VectorXd::Zero(9);
for (size_t i = 0; i != 3; ++i) {
......@@ -666,10 +667,10 @@ double_t get_maximum(Eigen::VectorXd value) {
+0.0, +0.0, +1.0,
+1.0, +1.0, +1.0;
Eigen::Vector3d y;
Eigen::VectorXd y(3);
y << value(i - 1), value(i), value(i + 1);
Eigen::Vector3d p = X.colPivHouseholderQr().solve(y);
Eigen::VectorXd p = X.colPivHouseholderQr().solve(y);
double_t maximum = p(2) - p(1) * p(1) / (4 * p(0));
......
......@@ -279,10 +279,10 @@ public:
std::vector<std::shared_ptr<const Contour>> creg;
Eigen::Vector4d vx_m;
Eigen::Vector4d vx_p;
Eigen::Vector4d vy_m;
Eigen::Vector4d vy_p;
Eigen::VectorXd vx_m = Eigen::VectorXd::Zero(4);
Eigen::VectorXd vx_p = Eigen::VectorXd::Zero(4);
Eigen::VectorXd vy_m = Eigen::VectorXd::Zero(4);
Eigen::VectorXd vy_p = Eigen::VectorXd::Zero(4);
template<size_t QuadratureOrder>
std::shared_ptr<FiniteElementMesh<1, QuadratureOrder>> get_mesh() {
......@@ -293,7 +293,7 @@ public:
TEST_F(SimpleSquare, derivative_1_temp) {
auto df = get_mesh<1>()->derivative();
Eigen::Vector2d tmp;
Eigen::VectorXd tmp;
tmp = df.dx(0).transpose() * vx_p;
EXPECT_NEAR(tmp(0), 1.0, FLT_EPSILON);
......@@ -331,7 +331,7 @@ TEST_F(SimpleSquare, derivative_1_temp) {
TEST_F(SimpleSquare, derivative_2_temp) {
auto df = get_mesh<2>()->derivative();
Eigen::Vector2d tmp;
Eigen::VectorXd tmp;
for (size_t i = 0; i != 3; ++i) {
tmp = df.dx(i).transpose() * vx_p;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment