element_p4dc.hpp

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

const int TypeOfFE_P4dcLagrange::nn[15][4] =  { 
                 { 0 ,0 ,0 ,0 } ,
                 { 1 ,1 ,1 ,1 } ,
                 { 2 ,2 ,2 ,2 } ,
                 { 1 ,1 ,1 ,2 } ,
                 { 1 ,1 ,2 ,2 } ,
                 { 1 ,2 ,2 ,2 } ,
                 { 0 ,2 ,2 ,2 } ,
                 { 0 ,0 ,2 ,2 } ,
                 { 0 ,0 ,0 ,2 } ,
                 { 0 ,0 ,0 ,1 } ,
                 { 0 ,0 ,1 ,1 } ,
                 { 0 ,1 ,1 ,1 } ,
                 { 0 ,1 ,2 ,2 } ,
                 { 0 ,1 ,1 ,2 } ,
                 { 0 ,0 ,1 ,2 } }
;
const int TypeOfFE_P4dcLagrange::aa[15][4] =  { 
                 { 0 ,1 ,2 ,3 } ,
                 { 0 ,1 ,2 ,3 } ,
                 { 0 ,1 ,2 ,3 } ,
                 { 0 ,1 ,2 ,0 } ,
                 { 0 ,1 ,0 ,1 } ,
                 { 0 ,0 ,1 ,2 } ,
                 { 0 ,0 ,1 ,2 } ,
                 { 0 ,1 ,0 ,1 } ,
                 { 0 ,1 ,2 ,0 } ,
                 { 0 ,1 ,2 ,0 } ,
                 { 0 ,1 ,0 ,1 } ,
                 { 0 ,0 ,1 ,2 } ,
                 { 0 ,0 ,0 ,1 } ,
                 { 0 ,0 ,1 ,0 } ,
                 { 0 ,1 ,0 ,0 } }
;
const int TypeOfFE_P4dcLagrange::ff[15] =  { 24 ,24 ,24 ,6 ,4 ,6 ,6 ,4 ,6 ,6 ,4 ,6 ,2 ,2 ,2 };
const int TypeOfFE_P4dcLagrange::il[15] =  { 4 ,0 ,0 ,0 ,0 ,0 ,1 ,2 ,3 ,3 ,2 ,1 ,1 ,1 ,2 };
const int TypeOfFE_P4dcLagrange::jl[15] =  { 0 ,4 ,0 ,3 ,2 ,1 ,0 ,0 ,0 ,1 ,2 ,3 ,1 ,2 ,1 };
const int TypeOfFE_P4dcLagrange::kl[15] =  { 0 ,0 ,4 ,1 ,2 ,3 ,3 ,2 ,1 ,0 ,0 ,0 ,2 ,1 ,1 };