auto_q_bases_3d_nodes.cpp

Go to the documentation of this file.

#include "auto_q_bases_3d_nodes.hpp"
 
 
namespace polyfem {
namespace autogen {
namespace {
void q_0_nodes_3d(Eigen::MatrixXd &res) {
 res.resize(1, 3); res << 
0.5, 0.5, 0.5;
}
 
void q_1_nodes_3d(Eigen::MatrixXd &res) {
 res.resize(8, 3); res << 
0, 0, 0,
1, 0, 0,
1, 1, 0,
0, 1, 0,
0, 0, 1,
1, 0, 1,
1, 1, 1,
0, 1, 1;
}
 
void q_2_nodes_3d(Eigen::MatrixXd &res) {
 res.resize(27, 3); res << 
0, 0, 0,
1, 0, 0,
1, 1, 0,
0, 1, 0,
0, 0, 1,
1, 0, 1,
1, 1, 1,
0, 1, 1,
1.0/2.0, 0, 0,
1, 1.0/2.0, 0,
1.0/2.0, 1, 0,
0, 1.0/2.0, 0,
0, 0, 1.0/2.0,
1, 0, 1.0/2.0,
1, 1, 1.0/2.0,
0, 1, 1.0/2.0,
1.0/2.0, 0, 1,
1, 1.0/2.0, 1,
1.0/2.0, 1, 1,
0, 1.0/2.0, 1,
0, 1.0/2.0, 1.0/2.0,
1, 1.0/2.0, 1.0/2.0,
1.0/2.0, 0, 1.0/2.0,
1.0/2.0, 1, 1.0/2.0,
1.0/2.0, 1.0/2.0, 0,
1.0/2.0, 1.0/2.0, 1,
1.0/2.0, 1.0/2.0, 1.0/2.0;
}
 
void q_3_nodes_3d(Eigen::MatrixXd &res) {
 res.resize(64, 3); res << 
0, 0, 0,
1, 0, 0,
1, 1, 0,
0, 1, 0,
0, 0, 1,
1, 0, 1,
1, 1, 1,
0, 1, 1,
1.0/3.0, 0, 0,
2.0/3.0, 0, 0,
1, 1.0/3.0, 0,
1, 2.0/3.0, 0,
2.0/3.0, 1, 0,
1.0/3.0, 1, 0,
0, 2.0/3.0, 0,
0, 1.0/3.0, 0,
0, 0, 1.0/3.0,
0, 0, 2.0/3.0,
1, 0, 2.0/3.0,
1, 0, 1.0/3.0,
1, 1, 2.0/3.0,
1, 1, 1.0/3.0,
0, 1, 2.0/3.0,
0, 1, 1.0/3.0,
1.0/3.0, 0, 1,
2.0/3.0, 0, 1,
1, 1.0/3.0, 1,
1, 2.0/3.0, 1,
2.0/3.0, 1, 1,
1.0/3.0, 1, 1,
0, 2.0/3.0, 1,
0, 1.0/3.0, 1,
0, 2.0/3.0, 2.0/3.0,
0, 2.0/3.0, 1.0/3.0,
0, 1.0/3.0, 2.0/3.0,
0, 1.0/3.0, 1.0/3.0,
1, 1.0/3.0, 1.0/3.0,
1, 1.0/3.0, 2.0/3.0,
1, 2.0/3.0, 1.0/3.0,
1, 2.0/3.0, 2.0/3.0,
1.0/3.0, 0, 1.0/3.0,
1.0/3.0, 0, 2.0/3.0,
2.0/3.0, 0, 1.0/3.0,
2.0/3.0, 0, 2.0/3.0,
1.0/3.0, 1, 1.0/3.0,
1.0/3.0, 1, 2.0/3.0,
2.0/3.0, 1, 1.0/3.0,
2.0/3.0, 1, 2.0/3.0,
1.0/3.0, 1.0/3.0, 0,
1.0/3.0, 2.0/3.0, 0,
2.0/3.0, 1.0/3.0, 0,
2.0/3.0, 2.0/3.0, 0,
1.0/3.0, 1.0/3.0, 1,
1.0/3.0, 2.0/3.0, 1,
2.0/3.0, 1.0/3.0, 1,
2.0/3.0, 2.0/3.0, 1,
1.0/3.0, 1.0/3.0, 1.0/3.0,
1.0/3.0, 1.0/3.0, 2.0/3.0,
1.0/3.0, 2.0/3.0, 1.0/3.0,
1.0/3.0, 2.0/3.0, 2.0/3.0,
2.0/3.0, 1.0/3.0, 1.0/3.0,
2.0/3.0, 1.0/3.0, 2.0/3.0,
2.0/3.0, 2.0/3.0, 1.0/3.0,
2.0/3.0, 2.0/3.0, 2.0/3.0;
}
 
void q_m2_nodes_3d(Eigen::MatrixXd &res) {
 res.resize(20, 3); res << 
0.0, 0.0, 0.0,
1.0, 0.0, 0.0,
1.0, 1.0, 0.0,
0.0, 1.0, 0.0,
0.0, 0.0, 1.0,
1.0, 0.0, 1.0,
1.0, 1.0, 1.0,
0.0, 1.0, 1.0,
0.5, 0.0, 0.0,
1.0, 0.5, 0.0,
0.5, 1.0, 0.0,
0.0, 0.5, 0.0,
0.0, 0.0, 0.5,
1.0, 0.0, 0.5,
1.0, 1.0, 0.5,
0.0, 1.0, 0.5,
0.5, 0.0, 1.0,
1.0, 0.5, 1.0,
0.5, 1.0, 1.0,
0.0, 0.5, 1.0;
}
 
}
 
void q_nodes_3d(const int q, Eigen::MatrixXd &val){
switch(q){
    case 0: q_0_nodes_3d(val); break;
    case 1: q_1_nodes_3d(val); break;
    case 2: q_2_nodes_3d(val); break;
    case 3: q_3_nodes_3d(val); break;
    case -2: q_m2_nodes_3d(val); break;
    default: assert(false);
}}
}}