p5tet.h

finite element library for mathematic majored research

    {-0.02147880466, 0.01509376536, 0.01307158424, 0.009243005859, 0, 0.01093647202, 0.009471262597, 0.008929592012, 0.1237320564, -0.06015059611, 0, 0, 0.006508728684, 0.006764069264, 0, -0.05749458874, 0.01269202093, 0.0266937239, 0.00637722566, 0.0005636724387, 0, 0, 0, 0, 0, 0, 0, -0.04317402181, -0.04486775961, 0, -0.01994122649, 0.02158701091, 0.01424670034, -0.06169002209, 0.04611408627, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.04512991348, -0.04690038185, 0, 0.03752030548, 0.01624676885, -0.01191368933, 0.02063511523, -0.01823903739},
    {0.07875561708, -0.03018753073, 0.01742877899, 0.05853903711, 0, -0.02187294404, -0.01894252519, -0.01785918402, 0.0103110047, 0.08384628549, 0, 0, 0, 0, 0, -0.005260942761, -0.004556110079, -0.0533874478, 0.03082325736, -0.03194143819, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.03489714636, -0.03022181527, -0.02849340069, 0.01645067256, -0.06439354388, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.04377368973, -0.03790912732, 0.02382737866, -0.01375674349, 0.03039839564},
    {0.07875561708, -0.03773441341, -0.05664353172, -0.003081001953, 0.042356774, 0.0546823601, 0.03788505039, 0.04464796006, 0.005155502352, -0.007290981347, 0.04175743862, -0.00504163956, -0.024733169, -0.02029220779, 0.02959305156, 0.0101461039, 0.002169576228, 0.0266937239, 0.01328588679, 0.01427970178, -0.05573742249, -0.02457789311, -0.005193027733, 0.004022501985, 0.004644785208, -0.002198313589, -0.01300539858, -0.01583047467, -0.01080149768, -0.01689450834, -0.0267960231, -0.01942830982, -0.02951102214, -0.02173838874, -0.0178640154, 0, 0, 0, 0, 0, 0, -0.05127097164, -0.02260837341, -0.004776890744, 0.00370016366, 0.00427258097, -0.03676644349, -0.007768329036, 0.006017321797, 0.006948204719, -0.009082219892, 0.007035057278, 0.008123384426, 0.006807822476, 0.007860996278, 0.006079679129},
    {-0.02147880466, 0, -0.004357194748, -0.02156701367, -0.127070322, 0, 0.02841378779, 0, 0.03608851646, 0.03098667072, 0, 0.04537475604, 0, -0.02029220779, -0.04220648338, 0, -0.00618329225, 0, 0.001594306415, 0.003945707071, 0.08360613374, 0, -0.01038605547, 0, 0.004644785208, 0, 0.02415288308, 0, -0.01080149768, 0.04344302144, 0, -0.01942830982, 0, -0.02173838874, -0.0178640154, 0, 0, 0, 0, 0, 0, 0.07690645746, 0, -0.009553781488, 0, 0.00427258097, 0, -0.01553665807, 0, 0.006948204719, -0.01816443978, 0, 0.008123384426, 0, 0.007860996278, 0.006079679129},
    {0.07875561708, 0.03773441341, -0.05664353172, -0.003081001953, 0.042356774, -0.0546823601, 0.03788505039, -0.04464796006, 0.005155502352, -0.007290981347, -0.04175743862, -0.00504163956, 0.024733169, -0.02029220779, 0.02959305156, -0.0101461039, 0.002169576228, -0.0266937239, 0.01328588679, 0.01427970178, -0.05573742249, 0.02457789311, -0.005193027733, -0.004022501985, 0.004644785208, 0.002198313589, -0.01300539858, 0.01583047467, -0.01080149768, -0.01689450834, 0.0267960231, -0.01942830982, 0.02951102214, -0.02173838874, -0.0178640154, 0, 0, 0, 0, 0, 0, -0.05127097164, 0.02260837341, -0.004776890744, -0.00370016366, 0.00427258097, 0.03676644349, -0.007768329036, -0.006017321797, 0.006948204719, -0.009082219892, -0.007035057278, 0.008123384426, -0.006807822476, 0.007860996278, 0.006079679129},
    {-0.02147880466, -0.0188672067, 0.006535792122, 0.009243005859, 0.01411892467, 0, -0.003157087532, -0.1116199001, -0.05413277469, -0.06015059611, -0.008948022562, -0.00252081978, 0.0006508728684, 0.001127344877, -0.02765252359, 0.03607503608, 0.03742518993, -0.01886970138, 0.01992883019, 0.0005636724387, 0, 0, 0, 0, 0, 0.05935446691, 0.01672122675, -0.004317402181, -0.007477959935, 0.003218001588, 0.0224338798, 0.01870874279, 0.04630177612, 0.04318301547, 0.04611408627, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.06204337339, 0.01747874033, -0.004512991348, -0.007816730308, 0.03632887957, -0.00938007637, -0.01624676885, -0.01191368933, -0.02063511523, -0.01823903739},
    {-0.02147880466, 0.0188672067, 0.006535792122, 0.009243005859, 0.01411892467, 0, -0.003157087532, 0.1116199001, -0.05413277469, -0.06015059611, 0.008948022562, -0.00252081978, -0.0006508728684, 0.001127344877, -0.02765252359, -0.03607503608, 0.03742518993, 0.01886970138, 0.01992883019, 0.0005636724387, 0, 0, 0, 0, 0, -0.05935446691, 0.01672122675, 0.004317402181, -0.007477959935, 0.003218001588, -0.0224338798, 0.01870874279, -0.04630177612, 0.04318301547, 0.04611408627, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.06204337339, 0.01747874033, 0.004512991348, -0.007816730308, 0.03632887957, 0.00938007637, -0.01624676885, 0.01191368933, -0.02063511523, -0.01823903739},
    {0.07875561708, 0, -0.03485755798, 0.05853903711, -0.02823784933, 0, 0.006314175065, 0, -0.02062200941, 0.08384628549, 0, 0, 0, 0, -0.006791847899, 0, 0.00151870336, 0, -0.06164651471, -0.03194143819, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.04505202223, 0, 0.01007393842, 0, -0.03290134512, -0.06439354388, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.05651159044, 0, 0.01263637577, 0, 0.02751348697, 0.03039839564},
    {0.07875561708, -0.03773441341, -0.02178597374, -0.0523770332, 0.042356774, 0.0546823601, 0, 0.04464796006, 0.02577751176, 0.02916392539, 0.04175743862, 0.02856929084, 0.02733666047, 0.01878908129, -0.009217507864, -0.007891414141, -0.003471321965, -0.02485277742, -0.01434875774, -0.01503126503, -0.05573742249, -0.004915578621, -0.01973350538, -0.03942051945, -0.03251349645, -0.02418144948, -0.009289570415, -0.01583047467, -0.01412503543, -0.002413501191, -0.0006231633279, -0.002158701091, 0.007123350172, 0.00411266814, 0.005400748842, 0, -0.05105866857, -0.04052656179, -0.01141706538, 0.01031755761, 0.01276466714, 0.004660997422, -0.01438714671, -0.005645416334, 0.005718434747, 0.007379912584, 0.003676644349, -0.001553665807, 0.002406928719, 0.003474102359, 0.001009135544, 0.0002605576769, 0.0009025982696, -0.0008509778095, -0.0004913122674, -0.0008685255898},
    {-0.02147880466, 0, -0.02178597374, 0.003081001953, -0.127070322, 0, 0.04735631299, 0, 0.02577751176, 0.01275921736, 0, -0.03529147692, 0, 0.002254689755, 0.05093885925, 0, -0.009437656592, 0, -0.01434875774, -0.01296446609, 0.08360613374, 0, 0.04777585514, 0, -0.03251349645, 0, 0.009289570415, 0, -0.01412503543, -0.01448100715, 0, -0.002158701091, 0, 0.00411266814, 0.005400748842, 0, 0.07658800285, 0, -0.02283413075, 0, 0.01276466714, -0.006991496133, 0, -0.01129083267, 0, 0.007379912584, 0, -0.003107331614, 0, 0.003474102359, 0.002018271087, 0, 0.0009025982696, 0, -0.0004913122674, -0.0008685255898},
    {0.07875561708, 0.03773441341, -0.02178597374, -0.0523770332, 0.042356774, -0.0546823601, 0, -0.04464796006, 0.02577751176, 0.02916392539, -0.04175743862, 0.02856929084, -0.02733666047, 0.01878908129, -0.009217507864, 0.007891414141, -0.003471321965, 0.02485277742, -0.01434875774, -0.01503126503, -0.05573742249, 0.004915578621, -0.01973350538, 0.03942051945, -0.03251349645, 0.02418144948, -0.009289570415, 0.01583047467, -0.01412503543, -0.002413501191, 0.0006231633279, -0.002158701091, -0.007123350172, 0.00411266814, 0.005400748842, 0, -0.05105866857, 0.04052656179, -0.01141706538, -0.01031755761, 0.01276466714, 0.004660997422, 0.01438714671, -0.005645416334, -0.005718434747, 0.007379912584, -0.003676644349, -0.001553665807, -0.002406928719, 0.003474102359, 0.001009135544, -0.0002605576769, 0.0009025982696, 0.0008509778095, -0.0004913122674, -0.0008685255898},
    {-0.02147880466, -0.0188672067, 0.01089298687, 0.003081001953, 0.01411892467, -0.1093647202, -0.07892718831, 0.02232398003, -0.01288875588, 0.01275921736, -0.008948022562, -0.02772901758, -0.005857855816, 0.01916486291, 0.001455395978, 0.03832972583, 0.03482169846, -0.01242638871, 0.007174378868, -0.01296446609, 0, 0.05898694345, 0.007270238826, 0.05068352501, 0.06502699291, -0.006594940767, 0.01486331266, -0.004317402181, -0.001661768874, -0.004022501985, -0.008101123263, -0.01151307248, 0.003561675086, -0.00205633407, 0.005400748842, 0, 0, 0.06078984269, 0.03995972881, -0.01547633642, -0.03191166785, 0, -0.01233184004, 0.003039839564, -0.003027406631, -0.007768329036, 0.001378741631, -0.00427258097, 0.0009025982696, 0.0008685255898, 0.001009135544, 0.0007816730308, 0.001805196539, -0.0004254889047, 0.0002456561337, -0.0008685255898},
    {-0.02147880466, 0.0188672067, 0.01089298687, 0.003081001953, 0.01411892467, 0.1093647202, -0.07892718831, -0.02232398003, -0.01288875588, 0.01275921736, 0.008948022562, -0.02772901758, 0.005857855816, 0.01916486291, 0.001455395978, -0.03832972583, 0.03482169846, 0.01242638871, 0.007174378868, -0.01296446609, 0, -0.05898694345, 0.007270238826, -0.05068352501, 0.06502699291, 0.006594940767, 0.01486331266, 0.004317402181, -0.001661768874, -0.004022501985, 0.008101123263, -0.01151307248, -0.003561675086, -0.00205633407, 0.005400748842, 0, 0, -0.06078984269, 0.03995972881, 0.01547633642, -0.03191166785, 0, 0.01233184004, 0.003039839564, 0.003027406631, -0.007768329036, -0.001378741631, -0.00427258097, -0.0009025982696, 0.0008685255898, 0.001009135544, -0.0007816730308, 0.001805196539, 0.0004254889047, 0.0002456561337, -0.0008685255898},
    {0.07875561708, 0, 0.04357194748, -0.0523770332, -0.02823784933, 0, 0.06314175065, 0, -0.05155502352, 0.02916392539, 0, -0.00672218608, 0, -0.06012506013, 0.0009702639856, 0, -0.01258354212, 0, 0.02869751547, -0.01503126503, 0, 0, -0.04362143295, 0, -0.06502699291, 0, 0.0111474845, 0, 0.03157360861, -0.001609000794, 0, -0.002878268121, 0, -0.008225336279, 0.005400748842, 0, 0, 0, -0.05327963842, 0, 0.04254889047, 0, 0, 0.01823903739, 0, -0.003884164518, 0, -0.004660997422, 0, -0.003474102359, 0.000672757029, 0, 0.001203464359, 0, 0.0009826245347, -0.0008685255898},
    {0.1073940233, -0.188672067, -0.1089298687, -0.07702504883, -0.07059462333, -0.0546823601, -0.03157087532, -0.04464796006, -0.02577751176, -0.01822745337, 0.08948022562, -0.01512491868, 0.05857855816, 0.0744047619, -0.01309856381, -0.02818362193, 0.009763093026, 0.06903549284, 0.03985766038, 0.0597492785, 0, 0.07373367932, 0.03635119413, -0.01206750595, -0.02322392604, 0.06594940767, 0.007431656332, 0.04317402181, 0.04652952848, 0.03459351707, 0.0155790832, 0.02518484606, -0.02544053633, -0.0146881005, -0.03281993527, 0, 0, 0, 0, 0, 0, 0, 0.06782512022, 0.03343823521, -0.01110049098, -0.02136290485, -0.01378741631, 0.01942082259, -0.009025982696, -0.01910756298, -0.009082219892, -0.002345019093, -0.008123384426, 0.003403911238, 0.001965249069, 0.006079679129},
    {0.1073940233, 0.188672067, -0.1089298687, -0.07702504883, -0.07059462333, 0.0546823601, -0.03157087532, 0.04464796006, -0.02577751176, -0.01822745337, -0.08948022562, -0.01512491868, -0.05857855816, 0.0744047619, -0.01309856381, 0.02818362193, 0.009763093026, -0.06903549284, 0.03985766038, 0.0597492785, 0, -0.07373367932, 0.03635119413, 0.01206750595, -0.02322392604, -0.06594940767, 0.007431656332, -0.04317402181, 0.04652952848, 0.03459351707, -0.0155790832, 0.02518484606, 0.02544053633, -0.0146881005, -0.03281993527, 0, 0, 0, 0, 0, 0, 0, -0.06782512022, 0.03343823521, 0.01110049098, -0.02136290485, 0.01378741631, 0.01942082259, 0.009025982696, -0.01910756298, -0.009082219892, 0.002345019093, -0.008123384426, -0.003403911238, 0.001965249069, 0.006079679129},
    {0.1073940233, 0, 0.2178597374, -0.07702504883, 0, 0, -0.09471262597, 0, 0.05155502352, -0.01822745337, 0, -0.09074951207, 0, -0.0946969697, 0.02328633566, 0, -0.02278055039, 0, -0.07971532075, 0.0597492785, 0, 0, -0.07270238826, 0, 0.04644785208, 0, -0.04830576616, 0, -0.0781031371, 0.01448100715, 0, 0.04317402181, 0, 0.029376201, -0.03281993527, 0, 0, 0, 0, 0, 0, 0, 0, -0.06687647042, 0, 0.0427258097, 0, 0.03107331614, 0, 0.006948204719, -0.006054813261, 0, -0.01083117923, 0, -0.003930498139, 0.006079679129},
    {0.1073940233, 0, 0, 0.2310751465, 0, 0, 0, 0, 0, -0.1093647202, 0, 0.01008327912, 0, -0.004509379509, -0.09314534262, 0, -0.1041396589, 0, 0, -0.09244227994, 0, 0, 0, 0, 0, 0, -0.06688490699, 0, 0.02991183974, -0.05792402858, 0, -0.06476103272, 0, 0, 0.08350388594, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.06991496133, 0, 0.03126692123, 0.02421925305, 0, 0.02707794809, 0, 0, -0.01823903739}};
    
    // Extract relevant coefficients
    const double coeff0_0 = coefficients0[dof][0];
    const double coeff0_1 = coefficients0[dof][1];
    const double coeff0_2 = coefficients0[dof][2];
    const double coeff0_3 = coefficients0[dof][3];
    const double coeff0_4 = coefficients0[dof][4];
    const double coeff0_5 = coefficients0[dof][5];
    const double coeff0_6 = coefficients0[dof][6];
    const double coeff0_7 = coefficients0[dof][7];
    const double coeff0_8 = coefficients0[dof][8];
    const double coeff0_9 = coefficients0[dof][9];
    const double coeff0_10 = coefficients0[dof][10];
    const double coeff0_11 = coefficients0[dof][11];
    const double coeff0_12 = coefficients0[dof][12];
    const double coeff0_13 = coefficients0[dof][13];
    const double coeff0_14 = coefficients0[dof][14];
    const double coeff0_15 = coefficients0[dof][15];
    const double coeff0_16 = coefficients0[dof][16];
    const double coeff0_17 = coefficients0[dof][17];
    const double coeff0_18 = coefficients0[dof][18];
    const double coeff0_19 = coefficients0[dof][19];
    const double coeff0_20 = coefficients0[dof][20];
    const double coeff0_21 = coefficients0[dof][21];
    const double coeff0_22 = coefficients0[dof][22];
    const double coeff0_23 = coefficients0[dof][23];
    const double coeff0_24 = coefficients0[dof][24];
    const double coeff0_25 = coefficients0[dof][25];
    const double coeff0_26 = coefficients0[dof][26];
    const double coeff0_27 = coefficients0[dof][27];
    const double coeff0_28 = coefficients0[dof][28];
    const double coeff0_29 = coefficients0[dof][29];
    const double coeff0_30 = coefficients0[dof][30];
    const double coeff0_31 = coefficients0[dof][31];
    const double coeff0_32 = coefficients0[dof][32];
    const double coeff0_33 = coefficients0[dof][33];
    const double coeff0_34 = coefficients0[dof][34];
    const double coeff0_35 = coefficients0[dof][35];
    const double coeff0_36 = coefficients0[dof][36];
    const double coeff0_37 = coefficients0[dof][37];
    const double coeff0_38 = coefficients0[dof][38];
    const double coeff0_39 = coefficients0[dof][39];
    const double coeff0_40 = coefficients0[dof][40];
    const double coeff0_41 = coefficients0[dof][41];
    const double coeff0_42 = coefficients0[dof][42];
    const double coeff0_43 = coefficients0[dof][43];
    const double coeff0_44 = coefficients0[dof][44];
    const double coeff0_45 = coefficients0[dof][45];
    const double coeff0_46 = coefficients0[dof][46];
    const double coeff0_47 = coefficients0[dof][47];
    const double coeff0_48 = coefficients0[dof][48];
    const double coeff0_49 = coefficients0[dof][49];
    const double coeff0_50 = coefficients0[dof][50];
    const double coeff0_51 = coefficients0[dof][51];
    const double coeff0_52 = coefficients0[dof][52];
    const double coeff0_53 = coefficients0[dof][53];
    const double coeff0_54 = coefficients0[dof][54];
    const double coeff0_55 = coefficients0[dof][55];
    
    // Compute value(s)
    *values = coeff0_0*basisvalue0 + coeff0_1*basisvalue1 + coeff0_2*basisvalue2 + coeff0_3*basisvalue3 + coeff0_4*basisvalue4 + coeff0_5*basisvalue5 + coeff0_6*basisvalue6 + coeff0_7*basisvalue7 + coeff0_8*basisvalue8 + coeff0_9*basisvalue9 + coeff0_10*basisvalue10 + coeff0_11*basisvalue11 + coeff0_12*basisvalue12 + coeff0_13*basisvalue13 + coeff0_14*basisvalue14 + coeff0_15*basisvalue15 + coeff0_16*basisvalue16 + coeff0_17*basisvalue17 + coeff0_18*basisvalue18 + coeff0_19*basisvalue19 + coeff0_20*basisvalue20 + coeff0_21*basisvalue21 + coeff0_22*basisvalue22 + coeff0_23*basisvalue23 + coeff0_24*basisvalue24 + coeff0_25*basisvalue25 + coeff0_26*basisvalue26 + coeff0_27*basisvalue27 + coeff0_28*basisvalue28 + coeff0_29*basisvalue29 + coeff0_30*basisvalue30 + coeff0_31*basisvalue31 + coeff0_32*basisvalue32 + coeff0_33*basisvalue33 + coeff0_34*basisvalue34 + coeff0_35*basisvalue35 + coeff0_36*basisvalue36 + coeff0_37*basisvalue37 + coeff0_38*basisvalue38 + coeff0_39*basisvalue39 + coeff0_40*basisvalue40 + coeff0_41*basisvalue41 + coeff0_42*basisvalue42 + coeff0_43*basisvalue43 + coeff0_44*basisvalue44 + coeff0_45*basisvalue45 + coeff0_46*basisvalue46 + coeff0_47*basisvalue47 + coeff0_48*basisvalue48 + coeff0_49*basisvalue49 + coeff0_50*basisvalue50 + coeff0_51*basisvalue51 + coeff0_52*basisvalue52 + coeff0_53*basisvalue53 + coeff0_54*basisvalue54 + coeff0_55*basisvalue55;
  }

  /// Evaluate order n derivatives of basis function i at given point in cell
  virtual void evaluate_basis_derivatives(unsigned int i,
                                          unsigned int n,
                                          double* values,
                                          const double* coordinates,
                                          const ufc::cell& c) const
  {
    // Extract vertex coordinates
    const double * const * element_coordinates = c.coordinates;
    
    // Compute Jacobian of affine map from reference cell
    const double J_00 = element_coordinates[1][0] - element_coordinates[0][0];
    const double J_01 = element_coordinates[2][0] - element_coordinates[0][0];
    const double J_02 = element_coordinates[3][0] - element_coordinates[0][0];
    const double J_10 = element_coordinates[1][1] - element_coordinates[0][1];
    const double J_11 = element_coordinates[2][1] - element_coordinates[0][1];
    const double J_12 = element_coordinates[3][1] - element_coordinates[0][1];
    const double J_20 = element_coordinates[1][2] - element_coordinates[0][2];
    const double J_21 = element_coordinates[2][2] - element_coordinates[0][2];
    const double J_22 = element_coordinates[3][2] - element_coordinates[0][2];
      
    // Compute sub determinants
    const double d00 = J_11*J_22 - J_12*J_21;
    const double d01 = J_12*J_20 - J_10*J_22;
    const double d02 = J_10*J_21 - J_11*J_20;
    
    const double d10 = J_02*J_21 - J_01*J_22;
    const double d11 = J_00*J_22 - J_02*J_20;
    const double d12 = J_01*J_20 - J_00*J_21;