ParMooN
 All Classes Functions Variables Friends Pages
BF_N_T_BDDF3_3D.h
1 // ***********************************************************************
2 // Brezzi-Douglas-Duran-Fortin element of third order on tetrahedra, 3D
3 // ***********************************************************************
4 
5 static double N_T_BDDF3_3D_CM[3600] = {
6  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,0,-0,0,0,0,0,0,-0,-10,10,-5,1,10,-5,1,-5,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
7  0,0,0,0,0,0,0,0,0,0,-10,10,-5,1,10,-5,1,-5,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
8  -10,10,-5,1,10,-5,1,-5,1,1,-0,0,-0,0,-0,0,-0,0,-0,-0,0,-0,0,0,0,-0,0,-0,-0,0,-0,0,-0,0,0,-0,0,-0,-0,0,0,-0,-0,-0,0,-0,-0,-0,0,-0,-0,-0,-0,0,0,-0,-0,-0,0,-0,
9  -3.26388898,9.750000085,-9.375000077,3.777777814,3.250000159,-4.666666723,0.1249999565,-1.500000264,2.500000142,-0.1805555173,-3.263888735,3.249999962,-1.499999823,-0.180555641,9.749999832,-4.666666707,2.499999976,-9.374999852,0.1250000018,3.777777725,-0.3055554459,6.916666651,-4.333333282,8.611111124,0.4166665162,-4.000000021,-4.333333393,0.416666829,6.916666734,-0.3055556694,62.52777771,-48.99999991,28.87499995,-7.597222151,-48.99999983,12.33333318,-2.625000184,28.8749999,-2.624999729,-7.597222284,1296.666668,-1843.333336,-1621.666673,-1621.666666,251.666674,-571.666677,-338.3333461,-326.6666724,251.6666642,-571.666664,-326.6666689,-338.3333265,6999.999997,-5179.999983,16520.00006,6999.999981,-5180.000015,30239.99979,13720.00018,8119.999869,
10  0,0,0,0,0,0,0,0,0,0,47,-27,6,0,-87,36,-6,51,-9,-11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
11  47,-87,51,-11,-27,36,-9,6,-6,-0,0,-0,0,-0,0,-0,0,-0,0,0,-0,0,-0,-0,-0,0,-0,0,-0,-0,0,-0,0,-0,-0,0,-0,0,-0,-0,-0,0,0,0,-0,0,0,0,-0,0,0,0,0,-0,-0,0,0,1e-10,-0,0,
12  0,0,-0,0,0,0,0,0,0,0,-0,0,0,0,0,0,0,0,0,0,-0,0,0,0,0,0,0,0,0,0,47,-87,51,-11,-27,36,-9,6,-6,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
13  -3.263888803,3.250000001,-1.500000079,-0.1805555274,9.749999426,-4.666666364,2.499999955,-9.374999572,0.1249998604,3.777777652,62.52777729,-48.99999988,28.87499982,-7.597222183,-48.99999964,12.33333347,-2.624999973,28.87499968,-2.625000043,-7.597222116,8.61111089,-4.333333252,6.916666565,-0.3055555447,-4.333333161,-3.999999973,0.4166666329,6.91666655,0.4166667029,-0.3055556081,-3.263888748,9.749999846,-9.374999856,3.777777641,3.249999632,-4.666666553,0.1250002085,-1.499999647,2.499999765,-0.1805556193,251.6666618,-338.3333274,-571.6666553,-326.6666595,1296.666645,-1621.666641,-1843.3333,-1621.666642,251.6666575,-326.666656,-571.6666532,-338.333319,7000.000061,-5179.999986,-8260.00031,-13999.99973,10359.9999,30239.99981,-52219.99955,21139.99976,
14  47,-27,6,0,-87,36,-6,51,-9,-11,0,-0,0,-0,-0,0,0,0,-0,-0,0,-0,0,-0,-0,-0,-0,0,0,-0,0,0,-0,0,-0,-0,-0,0,0,-0,-0,0,0,0,-0,0,0,0,-0,0,0,0,1e-10,-0,-1e-10,-0,0,1e-10,-1e-10,1e-10,
15  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,47,-27,6,-0,-87,36,-6,51,-9,-11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
16  0,0,0,0,0,0,0,0,0,0,47,-87,51,-11,-27,36,-9,6,-6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
17  62.52777777,-49.00000008,28.87500015,-7.597222285,-48.99999958,12.33333309,-2.62499991,28.87499983,-2.625000004,-7.597222133,-3.263888557,9.749999933,-9.375000011,3.777777825,3.249999804,-4.666666761,0.1250000004,-1.499999825,2.50000004,-0.1805556098,-0.3055554453,0.4166666008,0.4166667174,-0.3055555783,6.916666644,-4.000000005,6.916666754,-4.333333375,-4.33333343,8.611111268,-3.263888965,3.250000057,-1.50000009,-0.1805554904,9.750000213,-4.66666663,2.499999972,-9.375000271,0.1249999684,3.777777905,251.6666699,-338.3333369,-326.666672,-571.6666744,251.6666806,-326.666682,-338.3333541,-571.6666859,1296.666678,-1621.666679,-1621.666677,-1843.333354,-14000.00005,10359.99996,-8259.999748,6999.999739,-5179.999886,-60479.99958,38499.99935,-29259.99962,
18  9.758333594,-37.65000021,47.62500018,-22.66666675,-3.150000525,3.400000188,4.125000104,1.50000075,-4.500000371,0.7083332121,9.758332968,-3.149999997,1.499999613,0.7083335229,-37.64999957,3.400000134,-4.499999935,47.62499962,4.124999978,-22.66666653,4.383333047,-25.74999996,25.99999987,-42.3666667,-0.2499996065,6.00000005,26.00000014,-0.2500004085,-25.75000016,4.383333606,-110.1166665,75.59999974,-56.62499984,14.45833314,75.59999957,-2.599999637,1.875000448,-56.62499975,1.874999353,14.45833347,-3370.000004,5642.000006,3745.000015,3744.999998,-655.0000193,1855.000027,707.0000332,700.0000149,-654.9999949,1854.999995,700.0000066,706.9999848,-13271.99997,12179.99995,-28056.00019,-13271.99993,12180.00003,-72575.99949,-44520.00042,-13943.99968,
19  0,0,0,0,0,0,0,0,0,0,-72,18,-0,0,180,-36,0,-144,18,36,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
20  -72,180,-144,36,18,-36,18,-0,0,0,-0,0,-0,0,-0,0,-0,0,-0,-0,-0,0,-0,0,0,-0,0,-0,0,0,-0,0,-0,0,0,-0,0,-0,0,0,0,-0,-0,-0,0,-0,-0,-0,0,-0,-0,-0,-1e-10,0,1e-10,0,-0,-2e-10,1e-10,-1e-10,
21  11.53333362,-21.50909119,9.204545726,-0.6742425536,-27.69090958,37.60000015,-5.613636199,17.79545542,-21.38636412,0.3409089918,8.116665772,-5.009090377,-2.795455413,0.4090912566,-20.19090833,11.60000008,0.3863637228,8.045453891,-4.136363611,0.6742426225,-11.63333375,-17.77272716,-8.227272948,-9.300000027,13.77272774,18.00000007,8.227272902,-18.27272778,-13.72727292,6.533333665,-163.0499998,228.9909089,-95.0454545,32.99242405,81.80909022,-80.39999938,-4.863635829,-30.20454507,32.61363549,-0.7424240733,-2250.000006,2898.000008,4747.27275,2182.727272,-950.0000341,1692.727316,1918.000058,897.2727582,-759.9999935,1457.27272,1622.727286,727.9999799,-12768.00001,13439.99997,-71064.00025,-44687.99983,24360,-108863.9995,-37800.00035,-43481.45429,
22  11.53333293,-27.69090906,17.79545481,0.3409090085,-21.5090886,37.59999889,-21.38636352,9.204543638,-5.613635785,-0.6742419268,-163.0499982,81.80909044,-30.20454473,-0.7424244183,228.9909077,-80.40000053,32.61363627,-95.04545333,-4.863636198,32.99242383,-9.299999139,-8.22727305,-17.77272686,-11.63333338,8.227272066,17.99999986,13.77272737,-13.72727222,-18.27272735,6.533333453,8.116666059,-20.19090838,8.045453959,0.6742429261,-5.0090894,11.59999935,-4.136364349,-2.795455884,0.3863645919,0.4090911147,-949.9999813,1917.999978,1692.727228,897.2727015,-2249.999919,4747.27263,2897.999871,2182.727179,-759.9999673,1622.727235,1457.272675,727.9999508,-12768.00019,13439.99994,28056.00109,39311.99899,-37799.99965,-108863.9994,225959.9984,-37830.54465,
23  -144,198,-54,0,198,-252,54,-54,54,-0,-0,0,-0,0,0,0,0,-0,-0,0,-0,-0,-0,0,0,-0,-0,-0,0,0,-0,0,-0,0,0,-0,-0,-0,-0,0,-0,0,0,-0,0,-0,-0,-0,0,-0,-0,-0,-1e-10,0,1e-10,-0,-0,-1e-10,0,-0,
24  8.116666957,-20.1909094,8.045454852,0.6742422908,-5.009091277,11.60000015,-4.136363486,-2.795453751,0.3863631597,0.4090907914,11.53333311,-27.69090917,17.79545407,0.34090939,-21.50909052,37.60000014,-21.38636357,9.2045451,-5.613636367,-0.6742422751,6.533333007,-13.72727274,8.227272593,-9.300000065,-18.27272672,18.00000009,-8.227272475,13.77272665,-17.77272757,-11.63333284,-163.0499998,81.80909057,-30.20454525,-0.742424533,228.9909086,-80.39999952,32.61363714,-95.04545432,-4.863637448,32.99242456,-2250.000004,2898.000006,2182.72729,4747.272721,-760.0000159,1457.272753,728.0000289,1622.727281,-949.9999869,1692.727259,897.272728,1917.999969,-44688,24359.99991,-71064.00006,-12768.00007,13440.00011,-108863.999,-37800.001,-26070.54478,
25  0,0,0,0,0,0,0,0,0,0,-144,198,-54,0,198,-252,54,-54,54,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
26  -163.0500001,228.9909095,-95.04545515,32.99242449,81.80908965,-80.39999927,-4.863636691,-30.20454514,32.61363654,-0.7424245,11.53333222,-21.50909048,9.204545524,-0.6742426297,-27.69090858,37.60000026,-5.613636407,17.79545409,-21.38636376,0.3409092187,6.533333047,-18.27272702,13.77272711,-11.63333323,-13.72727285,18.00000004,-17.77272764,8.227273031,-8.227272343,-9.300000646,8.116666923,-5.009091063,-2.795454312,0.4090907906,-20.19090984,11.59999994,0.3863634585,8.045455473,-4.136363309,0.6742419197,-950.0000104,1918.00001,897.2727415,1692.727301,-760.0000426,1622.727316,728.0000626,1457.272791,-2250.000042,4747.272775,2182.727311,2898.000078,39312.00012,-37799.99984,28055.99929,-12767.99911,13439.99958,254015.9983,-136919.9975,65718.54398,
27  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-72,180,-144,36,18,-36,18,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
28  9.758333125,-3.150000009,1.500000208,0.7083332581,-37.64999861,3.399999216,-4.499999878,47.62499905,4.125000355,-22.6666664,-110.1166654,75.59999955,-56.62499943,14.45833316,75.59999915,-2.600000306,1.874999969,-56.62499924,1.875000066,14.4583331,-42.36666617,25.99999982,-25.74999978,4.38333331,25.99999959,5.999999949,-0.2499999157,-25.74999972,-0.2500000898,4.383333456,9.758333064,-37.64999978,47.62499979,-22.6666664,-3.149999199,3.399999801,4.124999534,1.499999206,-4.49999951,0.7083334877,-654.9999886,706.9999862,1854.999974,699.9999827,-3369.999949,3744.999939,5641.99992,3744.999938,-654.9999777,699.9999738,1854.999969,706.9999641,-13272.00016,12179.99996,-4115.999192,44687.99933,-24359.99977,-72575.99954,112979.9989,-58043.99938,
29  -72,18,-0,-0,180,-36,0,-144,18,36,-0,0,-0,0,0,-0,-0,-0,0,0,-0,0,-0,0,0,0,-0,-0,-0,0,-0,-0,0,-0,0,0,0,-0,-0,0,0,-0,0,-0,-0,-0,0,-0,0,-0,-0,-0,-1e-10,0,2e-10,0,-0,-1e-10,1e-10,-1e-10,
30  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-144,198,-54,-0,198,-252,54,-54,54,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
31  8.116666457,-5.009090934,-2.795454277,0.4090908041,-20.19090729,11.59999905,0.3863638144,8.045453118,-4.13636325,0.6742429166,-163.0499985,228.9909091,-95.04545433,32.99242432,81.80908963,-80.40000045,-4.863636533,-30.20454437,32.61363658,-0.742424633,-9.299999208,8.227272439,-13.72727237,6.533333291,-8.227273316,17.99999992,-18.27272713,-17.77272694,13.7727271,-11.63333307,11.53333278,-21.50909021,9.204544733,-0.6742418169,-27.69090779,37.59999969,-5.613637221,17.79545323,-21.38636276,0.3409093416,-759.9999826,727.999979,1457.272687,1622.727247,-2249.999928,2182.727185,2897.999887,4747.272645,-949.9999674,897.272689,1692.727225,1917.999948,-44688.00024,24359.99996,79296.00102,39311.99905,-37799.99965,-108863.9993,174719.9984,-99761.45377,
32  -163.0499998,81.80909119,-30.20454613,-0.742423965,228.9909068,-80.39999858,32.61363595,-95.04545352,-4.86363658,32.99242381,8.116665109,-20.1909088,8.045454336,0.6742424085,-5.009089961,11.60000054,-4.136363665,-2.795455443,0.3863634767,0.4090911932,-11.6333339,13.77272752,-18.27272749,6.533333404,-17.77272693,18.00000003,-13.72727306,-8.22727279,8.227273087,-9.300000559,11.53333364,-27.69090928,17.7954549,0.3409087473,-21.50909186,37.59999979,-21.38636323,9.204546681,-5.613636484,-0.6742429042,-760.0000162,728.0000187,1622.727301,1457.272762,-950.0000697,897.2728067,1918.000103,1692.727367,-2250.00005,2182.727331,4747.272774,2898.00009,39312.00031,-37799.99989,79295.99862,-44687.99881,24359.99954,254015.9985,-188159.9974,159569.453,
33  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,0,-0,0,0,0,0,0,0,-72,18,0,0,180,-36,-0,-144,18,36,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
34  0,0,0,0,0,0,0,0,0,0,-72,180,-144,36,18,-36,18,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
35  -110.1166667,75.60000014,-56.62500026,14.45833345,75.59999927,-2.5999996,1.874999835,-56.62499967,1.875000046,14.45833314,9.758332679,-37.64999997,47.62500015,-22.66666685,-3.149999554,3.40000016,4.125000037,1.499999631,-4.5000001,0.7083334531,4.383333101,-0.2499998521,-0.2500001131,4.383333384,-25.74999998,5.999999995,-25.75000018,26.00000014,26.00000022,-42.36666702,9.75833351,-3.150000173,1.500000251,0.7083331739,-37.65000044,3.399999912,-4.499999952,47.62500056,4.125000091,-22.66666695,-655.0000065,707.000007,700.0000111,1855.000015,-655.0000277,700.0000307,707.0000421,1855.000037,-3370.000021,3745.000024,3745.000022,5642.000039,44688.00009,-24359.99992,-4116.000469,-13271.99948,12179.99976,108863.9991,-68459.99865,62915.99924,
36  -7.494444613,32.90000013,-48.2500001,28.88888894,-1.099999633,6.266666534,-14.25000006,-1.000000486,7.00000023,-1.527777695,-7.494444234,-1.099999965,-0.9999997897,-1.527777882,32.89999974,6.266666573,6.999999959,-48.24999976,-14.24999998,28.8888888,-5.077777601,23.83333331,-31.66666659,43.75555558,-1.16666691,2.999999971,-31.66666675,-1.16666642,23.83333343,-5.077777936,57.58888877,-36.59999983,32.74999989,-7.861110994,-36.59999975,-4.733333544,-0.2500002641,32.74999985,-0.2499996239,-7.861111188,2073.333336,-3798.66667,-2123.333343,-2123.333332,403.3333452,-1283.33335,-368.6666871,-373.3333425,403.3333307,-1283.333331,-373.3333377,-368.6666583,6271.999978,-6999.999972,11536.00014,6271.999953,-7000.000019,42335.99971,30800.00024,5823.99981,
37  0,0,0,0,0,0,0,0,0,0,36,-0,0,0,-108,0,-0,108,-0,-36,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
38  36,-108,108,-36,-0,0,0,0,-0,-0,0,-0,0,-0,0,-0,0,-0,0,-0,0,-0,0,-0,-0,0,-0,0,0,-0,0,-0,0,-0,-0,0,-0,0,0,-0,-0,0,0,0,-0,0,-0,0,-0,0,0,0,1e-10,-0,-0,-0,0,1e-10,-1e-10,1e-10,
39  -15.90000037,41.92727307,-30.98863667,-0.5189392464,26.47272797,-52.20000025,39.96590891,-9.511364748,0.534091505,2.018939545,-7.441665699,2.177272293,5.011364577,2.18939357,24.22272638,-9.200000155,-11.03409104,-12.88636286,16.65909092,-14.14772753,11.18333385,7.568181698,31.43181843,17.01666671,-9.568182432,-30.00000007,-23.4318184,18.31818245,29.68181841,-11.4000004,134.6416664,-158.3227269,66.8863635,-35.35227246,-64.27727191,39.79999931,19.34090841,39.98863588,-33.96590805,-0.6893941293,3055.000007,-4571.00001,-5631.818209,-2803.181817,1265.00004,-3188.181871,-1981.000068,-976.818216,899.9999926,-2096.818174,-1683.181832,-755.9999768,12095.99998,-15119.99995,63588.00041,43175.99979,-28140.00001,127007.9993,63420.00044,39808.36329,
40  -10.13333301,31.17272722,-26.38636381,-0.8560605571,12.627271,-33.39999917,37.15909084,-0.6136349625,-10.15909137,2.189393561,104.6166653,-27.32727241,3.613635774,2.393939534,-201.8727262,50.60000041,-10.840909,99.1363627,-12.40909105,-13.81060574,5.366665997,-1.318181579,27.3181815,12.36666669,-2.681817634,-20.99999989,-19.31818186,9.181817756,22.81818185,-8.133333379,-8.383332843,19.17272669,-7.636363173,6.643939,3.627271538,-9.399999445,-1.090908518,2.386364661,4.840908304,-1.060606197,669.9999861,-1861.999983,-738.1817838,-451.8181638,789.9999389,-2341.818107,-1021.999902,-528.1817498,519.9999772,-1368.181792,-871.8181433,-391.9999668,5712.000145,-8399.999965,-17304.00087,-24527.99925,29399.99975,72575.99958,-182279.9989,18067.63582,
41  108,-216,108,-0,-108,216,-108,0,-0,-0,0,-0,0,-0,0,-0,-0,-0,0,0,-0,0,-0,-0,-0,0,0,0,-0,0,0,-0,0,-0,-0,0,0,0,-0,-0,0,-0,-0,-0,-0,0,0,0,0,-0,-0,-0,-1e-10,0,-1e-10,1e-10,-0,-2e-10,2e-10,-2e-10,
42  -7.441667105,24.22272765,-12.88636397,-14.14772712,2.177273514,-9.200000279,16.65909072,5.011362412,-11.03409028,2.189394139,-15.89999968,26.47272748,-9.511363128,2.018939061,41.92727219,-52.20000022,0.5340908323,-30.98863587,39.96590911,-0.5189395958,-11.39999957,29.68181817,-23.431818,17.01666674,18.31818113,-30.00000011,31.4318179,-9.568181061,7.568182157,11.18333278,134.6416664,-64.27727225,39.98863604,-0.6893935829,-158.3227267,39.79999942,-33.96590997,66.88636334,19.34091032,-35.35227307,3055.000006,-4571.000009,-2803.181842,-5631.818176,900.0000237,-2096.818219,-756.000042,-1683.181833,1264.999987,-3188.181805,-976.8181875,-1980.999967,43175.99998,-28139.99989,63588.00013,12096.00003,-15120.00011,127007.9989,63420.00108,22855.6356,
43  0,0,0,0,0,0,0,0,0,0,108,-108,-0,0,-216,216,0,108,-108,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
44  104.6166668,-201.8727276,99.1363641,-13.81060625,-27.3272721,50.59999957,-12.40909065,3.613636421,-10.84090932,2.393939525,-10.13333259,12.62727232,-0.6136364158,2.189394113,31.17272701,-33.40000013,-10.15909088,-26.38636341,37.15909099,-0.8560606632,-8.133333195,22.81818164,-19.31818173,12.36666658,9.181818362,-21.00000002,27.31818211,-2.681818493,-1.318182118,5.366667179,-8.38333349,3.627272764,2.386363543,-1.060606027,19.17272778,-9.399999937,4.840909313,-7.636364279,-1.090909429,6.643939781,670.0000063,-1862.000006,-451.8181881,-738.1818386,520.0000254,-1368.181842,-392.0000364,-871.8182226,790.0000302,-2341.818216,-528.181843,-1022.000057,-24528.00007,29399.99988,-17303.99954,5711.999384,-8399.999692,-199583.9987,89879.99817,-35539.63524,
45  -10.13333353,12.62727291,-0.6136365457,2.189394033,31.17272767,-33.40000009,-10.15909105,-26.38636434,37.15909127,-0.8560605342,-8.383332443,3.627272131,2.386364408,-1.060606339,19.17272661,-9.400000072,4.840909047,-7.636363054,-1.090909122,6.643939238,12.36666701,27.31818167,-1.318181605,5.366666669,-19.31818209,-21.00000007,-2.681818287,22.81818215,9.181818307,-8.133333535,104.6166665,-201.8727272,99.13636365,-13.81060598,-27.32727218,50.5999995,-12.40909114,3.61363605,-10.84090859,2.393939323,790.0000058,-1022.000007,-2341.818201,-528.1818206,670.0000324,-738.1818567,-1862.000055,-451.8182134,519.9999979,-871.8181789,-1368.181834,-391.9999889,5712.000026,-8399.999994,53256.00017,46031.99982,-20999.99997,72575.99966,19320.00011,35020.36354,
46  -15.89999956,26.47272727,-9.511363966,2.018939495,41.92727011,-52.19999871,0.5340907765,-30.98863436,39.96590842,-0.5189399125,134.6416644,-64.27727195,39.98863534,-0.68939365,-158.3227257,39.8000006,-33.96590904,66.88636222,19.34090895,-35.35227228,17.01666571,31.43181857,7.568181338,11.18333339,-23.4318175,-29.99999984,-9.568181939,29.68181764,18.31818191,-11.40000015,-7.441666073,24.22272667,-12.88636315,-14.14772773,2.177271115,-9.199999342,16.65909155,5.011365085,-11.03409179,2.189393671,1264.999979,-1980.999975,-3188.181768,-976.8181511,3054.999904,-5631.818069,-4570.999848,-2803.181706,899.9999614,-1683.181774,-2096.81812,-755.9999415,12096.00021,-15119.99993,-6132.001256,-64343.99884,43259.9996,127007.9993,-229739.9982,53935.63537,
47  108,-108,0,-0,-216,216,-0,108,-108,0,0,-0,0,-0,-0,0,-0,0,-0,-0,0,-0,0,-0,-0,-0,-0,0,0,-0,0,0,-0,0,-0,-0,-0,0,0,-0,0,-0,-0,0,0,-0,-0,0,-0,0,0,0,1e-10,-1e-10,-0,-0,0,2e-10,-1e-10,1e-10,
48  -8.400000501,16.20000057,-5.537e-07,-2.999999761,16.20000046,-25.20000021,-2.477e-07,-1.1053e-06,7.885e-07,-2.999999896,-8.399999232,16.19999953,1.1027e-06,-3.000000581,16.19999909,-25.19999997,-1.599e-07,6.975e-07,-8.4e-08,-3.000000252,9.600000519,3.37e-08,1.931e-07,9.600000101,-9.355e-07,17.99999991,-3.711e-07,9.976e-07,4.006e-07,9.59999928,207.5999996,-226.7999995,26.99999984,-2.999999527,-226.7999991,244.799999,-27.00000145,26.99999955,-26.99999799,-3.000000539,1440.000006,-2016.00001,-2520.000028,-2519.99999,1440.000029,-2520.000046,-2016.000049,-2520.000021,1439.999978,-2519.999976,-2519.999999,-2015.99995,42336.00005,-30239.99992,157248.0003,42335.99996,-30240.0001,217727.9989,60480.0011,84671.99926,
49  -8.399999342,16.2,-5.323e-07,-2.999999821,16.199996,-25.1999981,-2.528e-07,3.1108e-06,-8.666e-07,-3.000000933,207.5999971,-226.7999996,26.99999936,-2.999999926,-226.7999977,244.8000008,-26.99999983,26.99999804,-27.0000003,-2.999999317,9.599998559,5.245e-07,-6.901e-07,9.600000087,1.0745e-06,18.00000026,-2.446e-07,-7.721e-07,2.323e-07,9.599999652,-8.399998893,16.19999855,1.2743e-06,-3.000001021,16.19999742,-25.19999896,1.4166e-06,2.3728e-06,-1.7856e-06,-3.000000358,1439.999967,-2015.999961,-2519.999924,-2519.999953,1439.999868,-2519.999839,-2015.999791,-2519.999843,1439.99994,-2519.99993,-2519.999908,-2015.999906,42336.0003,-30239.99988,-114912.0017,-48383.99824,60479.99934,217727.9986,-393119.997,84671.99849,
50  207.5999998,-226.8000007,27.00000126,-3.000000497,-226.7999962,244.7999979,-26.99999937,26.99999823,-26.99999988,-2.99999926,-8.399997434,16.19999939,1.085e-07,-2.999999821,16.19999867,-25.2000009,1.606e-07,1.3111e-06,2.588e-07,-3.00000039,9.600000873,-5.506e-07,4.6e-07,9.599999811,-1.26e-07,17.99999983,6.649e-07,-1.782e-07,-6.419e-07,9.600001086,-8.40000059,16.20000049,-6.687e-07,-2.999999615,16.20000169,-25.19999995,7.08e-08,-2.0153e-06,-3.182e-07,-2.999999063,1440.000028,-2016.00003,-2520.000046,-2520.000063,1440.000112,-2520.000121,-2016.000169,-2520.000157,1440.000091,-2520.000104,-2520.000093,-2016.000162,-48384.00035,60479.99976,-114911.9983,42335.99808,-30239.99916,-508031.9972,332639.9954,-187487.9973,
51  -8.383333498,19.17272748,-7.636363869,6.643939491,3.627272863,-9.400000065,-1.090909199,2.386363135,4.840909413,-1.060605981,-10.1333333,31.17272743,-26.38636333,-0.8560608176,12.62727253,-33.40000013,37.15909086,-0.6136361649,-10.15909089,2.189393842,-8.133333107,9.181818193,-2.681818086,5.36666672,22.81818142,-21.00000008,-1.318182025,-19.31818135,27.31818208,12.36666627,104.6166665,-27.32727249,3.613636195,2.393939606,-201.872727,50.59999969,-10.84090965,99.13636352,-12.40909013,-13.81060632,790.0000023,-1022.000004,-528.1818305,-2341.818176,520.0000087,-871.8181977,-392.0000173,-1368.181819,669.999989,-738.1818072,-451.8181807,-1861.999974,46031.99999,-20999.99991,53255.99997,5712.000122,-8400.000115,72575.99912,19320.00093,13027.63575,
52  0,0,0,0,0,0,0,0,0,0,108,-216,108,-0,-108,216,-108,0,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
53  134.6416667,-158.3227276,66.88636421,-35.35227296,-64.27727151,39.79999933,19.34090941,39.98863602,-33.96590929,-0.6893936697,-15.89999886,41.92727236,-30.98863657,-0.5189391032,26.47272667,-52.2000002,39.9659091,-9.511363151,0.5340910328,2.018939255,-11.39999971,18.31818157,-9.568181656,11.18333322,29.68181832,-30.00000002,7.568182206,-23.43181854,31.43181775,17.01666737,-7.441666959,2.17727295,5.011363339,2.189394094,24.22272805,-9.199999939,-11.03409076,-12.8863646,16.65909059,-14.14772674,1265.00001,-1981.00001,-976.8181959,-3188.181847,900.0000422,-1683.181862,-756.0000621,-2096.818245,3055.000041,-5631.818228,-2803.181857,-4571.000077,-64344.0001,43259.99982,-6131.999299,12095.99908,-15119.99955,-235871.9981,133139.9974,-65191.63482,
54  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,0,-0,0,0,0,0,0,0,36,-108,108,-36,-0,0,0,0,-0,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
55  -7.494444322,-1.099999992,-1.000000129,-1.527777731,32.89999918,6.266667148,6.999999922,-48.24999948,-14.25000022,28.88888875,57.58888815,-36.59999967,32.74999961,-7.861110975,-36.59999952,-4.733333159,-0.2499999961,32.74999955,-0.2500000237,-7.861110979,43.75555528,-31.66666657,23.83333321,-5.077777766,-31.66666643,3.000000024,-1.166666717,23.83333317,-1.166666613,-5.077777848,-7.494444317,32.89999994,-48.24999993,28.88888876,-1.100000433,6.266666753,-14.24999974,-0.9999995593,6.999999744,-1.527777868,403.3333269,-368.6666588,-1283.333319,-373.3333232,2073.333303,-2123.333298,-3798.66662,-2123.333296,403.3333202,-373.3333178,-1283.333316,-368.6666451,6272.000094,-6999.999978,12375.9995,-30687.9996,13999.99987,42335.99973,-60759.99935,36903.99961,
56  36,-0,0,0,-108,0,-0,108,-0,-36,0,-0,0,-0,-0,0,0,0,-0,-0,-0,-0,0,-0,-0,-0,0,0,0,-0,0,0,-0,0,-0,-0,-0,0,0,-0,-0,0,-0,0,0,-0,-0,-0,-0,0,-0,0,1e-10,-0,-1e-10,-0,0,-0,-0,0,
57  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,108,-216,108,0,-108,216,-108,-0,0,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
58  -7.441666381,2.177272747,5.011363281,2.189394082,24.22272505,-9.199998683,-11.03409117,-12.88636212,16.65909038,-14.14772777,134.6416648,-158.3227269,66.88636307,-35.35227263,-64.27727126,39.80000051,19.34090923,39.98863509,-33.96590928,-0.6893935061,17.0166658,-23.4318179,29.68181783,-11.39999997,31.43181894,-29.99999995,18.31818165,7.568181373,-9.568181608,11.18333305,-15.89999948,41.92727219,-30.98863577,-0.5189400246,26.47272588,-52.19999978,39.96591016,-9.511362173,0.5340898899,2.018939125,899.9999806,-755.9999761,-2096.818138,-1683.181789,3054.999916,-2803.181716,-4570.999871,-5631.818081,1264.999961,-976.8181359,-3188.181768,-1980.999936,43176.00031,-28139.99995,-39312.00142,-64343.99884,43259.9996,127007.9991,-196559.998,130948.3626,
59  104.6166665,-27.32727292,3.613636925,2.39393916,-201.8727252,50.59999865,-10.84090872,99.13636278,-12.40909062,-13.81060573,-8.383332003,19.17272695,-7.636363282,6.643939285,3.627271913,-9.400000463,-1.09090906,2.386364414,4.840909199,-1.060606307,12.36666714,-19.31818197,22.81818196,-8.13333337,27.31818139,-21,9.181818434,-1.318181648,-2.681818461,5.366667072,-10.13333355,31.17272734,-26.38636384,-0.8560602985,12.62727349,-33.39999977,37.15909036,-0.6136373836,-10.15909064,2.189394283,520.0000134,-392.0000161,-1368.181842,-871.8182095,670.0000592,-451.8182507,-1862.000085,-738.1818991,790.000041,-528.181867,-2341.818217,-1022.000075,-24528.0003,29399.99992,-90383.99866,46031.99898,-20999.99963,-199583.9988,162959.9979,-145732.3624,
60  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,0,-0,0,0,0,0,0,0,108,-108,-0,0,-216,216,0,108,-108,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
61  -8.383333274,3.627272759,2.386363478,-1.060605993,19.17272623,-9.399999463,4.840908972,-7.636362676,-1.090909288,6.643939015,104.6166657,-201.8727275,99.13636372,-13.81060627,-27.32727178,50.60000033,-12.40909073,3.613635582,-10.8409093,2.39393969,5.366666075,-2.681817946,9.181817902,-8.133333297,-1.318181435,-20.99999996,22.81818171,27.31818166,-19.31818167,12.36666643,-10.13333291,12.62727214,-0.6136357118,2.189393447,31.17272633,-33.39999983,-10.1590903,-26.38636267,37.15909033,-0.8560608211,519.9999874,-391.9999848,-871.8181519,-1368.1818,789.9999497,-528.1817557,-1021.999918,-2341.818127,669.9999774,-451.8181552,-738.1817826,-1861.999965,46032.00018,-20999.99999,-90384.00065,-24527.99936,29399.99976,72575.99953,-109199.9989,87940.36317,
62  134.6416666,-64.277273,39.98863697,-0.6893941961,-158.3227253,39.79999879,-33.9659087,66.88636266,19.34090922,-35.35227227,-7.441665162,24.22272723,-12.88636368,-14.14772708,2.177271703,-9.200000525,16.65909086,5.011364567,-11.03409069,2.189393632,11.18333394,-9.56818211,18.31818207,-11.40000008,7.568181532,-30.00000001,29.68181852,31.43181814,-23.43181857,17.01666728,-15.90000035,26.47272757,-9.511364135,2.018939788,41.92727369,-52.19999975,0.5340905525,-30.98863762,39.96590912,-0.518938861,900.0000162,-756.0000187,-1683.181848,-2096.818215,1265.000069,-976.8182614,-1981.000105,-3188.181908,3055.000048,-2803.181875,-5631.818231,-4571.000086,-64344.00031,43259.9999,-39311.99871,43175.99884,-28139.99954,-235871.9985,166319.9974,-173284.3622,
63  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,36,0,-0,-0,-108,-0,0,108,0,-36,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,
64  0,0,0,0,0,0,0,0,0,0,36,-108,108,-36,-0,0,0,0,-0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-0,0,0,
65  57.5888889,-36.60000006,32.75000011,-7.861111162,-36.59999969,-4.733333488,-0.2499999243,32.74999984,-0.2500000418,-7.86111101,-7.494444123,32.90000003,-48.25000014,28.88888902,-1.10000025,6.2666666,-14.25000004,-0.9999998051,7.00000006,-1.527777843,-5.077777656,-1.166666749,-1.166666604,-5.077777805,23.83333334,3.00000001,23.83333343,-31.66666676,-31.66666679,43.75555575,-7.494444545,-1.099999885,-1.00000016,-1.527777684,32.90000022,6.266666718,6.99999998,-48.25000029,-14.25000006,28.88888904,403.3333366,-368.6666701,-373.333339,-1283.333341,403.333347,-373.3333487,-368.666688,-1283.333351,2073.333343,-2123.333345,-2123.333345,-3798.666685,-30688.00004,13999.99996,12376.00022,6271.999741,-6999.999877,-48383.99952,29959.99929,-33655.99962
66 };
67 
68 static void N_T_BDDF3_3D_Funct(double xi, double eta, double zeta,
69  double *values)
70 {
71  int nBF = 60; // number of basis functions
72  // monomials x-component, y-component and z-component
73  double mon_x[]={1,0,0,xi,0,0,eta,0,0,zeta,0,0,
74  xi*xi,0,0,xi*eta,0,0,xi*zeta,0,0,
75  eta*eta,0,0,eta*zeta,0,0,zeta*zeta,0,0,
76  xi*xi*xi,0,0,xi*xi*eta,0,0,xi*xi*zeta,0,0,
77  xi*eta*eta,0,0,xi*eta*zeta,0,0,xi*zeta*zeta,0,0,
78  eta*eta*eta,0,0,eta*eta*zeta,0,0,eta*zeta*zeta,0,0,
79  zeta*zeta*zeta,0,0};
80  double mon_y[]={0,1,0,0,xi,0,0,eta,0,0,zeta,0,
81  0,xi*xi,0,0,xi*eta,0,0,xi*zeta,0,
82  0,eta*eta,0,0,eta*zeta,0,0,zeta*zeta,0,
83  0,xi*xi*xi,0,0,xi*xi*eta,0,0,xi*xi*zeta,0,
84  0,xi*eta*eta,0,0,xi*eta*zeta,0,0,xi*zeta*zeta,0,
85  0,eta*eta*eta,0,0,eta*eta*zeta,0,0,eta*zeta*zeta,0,
86  0,zeta*zeta*zeta,0};
87  double mon_z[]={0,0,1,0,0,xi,0,0,eta,0,0,zeta,
88  0,0,xi*xi,0,0,xi*eta,0,0,xi*zeta,
89  0,0,eta*eta,0,0,eta*zeta,0,0,zeta*zeta,
90  0,0,xi*xi*xi,0,0,xi*xi*eta,0,0,xi*xi*zeta,
91  0,0,xi*eta*eta,0,0,xi*eta*zeta,0,0,xi*zeta*zeta,
92  0,0,eta*eta*eta,0,0,eta*eta*zeta,0,0,eta*zeta*zeta,
93  0,0,zeta*zeta*zeta};
94 
95  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
96  for(int i=0; i<nBF; i++)
97  {
98  for(int j=0; j<nBF; j++)
99  {
100  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
101  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
102  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
103  }
104  }
105 }
106 
107 static void N_T_BDDF3_3D_DeriveXi(double xi, double eta, double zeta,
108  double *values)
109 {
110  int nBF = 60; // number of basis functions
111  // monomials x-component, y-component and z-component
112  double mon_x[]={0,0,0,1,0,0,0,0,0,0,0,0,
113  2*xi,0,0,eta,0,0,zeta,0,0,
114  0,0,0,0,0,0,0,0,0,
115  3*xi*xi,0,0,2*xi*eta,0,0,2*xi*zeta,0,0,
116  eta*eta,0,0,eta*zeta,0,0,zeta*zeta,0,0,
117  0,0,0,0,0,0,0,0,0,
118  0,0,0};
119  double mon_y[]={0,0,0,0,1,0,0,0,0,0,0,0,
120  0,2*xi,0,0,eta,0,0,zeta,0,
121  0,0,0,0,0,0,0,0,0,
122  0,3*xi*xi,0,0,2*xi*eta,0,0,2*xi*zeta,0,
123  0,eta*eta,0,0,eta*zeta,0,0,zeta*zeta,0,
124  0,0,0,0,0,0,0,0,0,
125  0,0,0};
126  double mon_z[]={0,0,0,0,0,1,0,0,0,0,0,0,
127  0,0,2*xi,0,0,eta,0,0,zeta,
128  0,0,0,0,0,0,0,0,0,
129  0,0,3*xi*xi,0,0,2*xi*eta,0,0,2*xi*zeta,
130  0,0,eta*eta,0,0,eta*zeta,0,0,zeta*zeta,
131  0,0,0,0,0,0,0,0,0,
132  0,0,0};
133 
134  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
135  for(int i=0; i<nBF; i++)
136  {
137  for(int j=0; j<nBF; j++)
138  {
139  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
140  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
141  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
142  }
143  }
144 }
145 
146 static void N_T_BDDF3_3D_DeriveEta(double xi, double eta, double zeta,
147  double *values)
148 {
149  int nBF = 60; // number of basis functions
150  // monomials x-component, y-component and z-component
151  double mon_x[]={0,0,0,0,0,0,1,0,0,0,0,0,
152  0,0,0,xi,0,0,0,0,0,
153  2*eta,0,0,zeta,0,0,0,0,0,
154  0,0,0,xi*xi,0,0,0,0,0,
155  xi*2*eta,0,0,xi*zeta,0,0,0,0,0,
156  3*eta*eta,0,0,2*eta*zeta,0,0,zeta*zeta,0,0,
157  0,0,0};
158  double mon_y[]={0,0,0,0,0,0,0,1,0,0,0,0,
159  0,0,0,0,xi,0,0,0,0,
160  0,2*eta,0,0,zeta,0,0,0,0,
161  0,0,0,0,xi*xi,0,0,0,0,
162  0,xi*2*eta,0,0,xi*zeta,0,0,0,0,
163  0,3*eta*eta,0,0,2*eta*zeta,0,0,zeta*zeta,0,
164  0,0,0};
165  double mon_z[]={0,0,0,0,0,0,0,0,1,0,0,0,
166  0,0,0,0,0,xi,0,0,0,
167  0,0,2*eta,0,0,zeta,0,0,0,
168  0,0,0,0,0,xi*xi,0,0,0,
169  0,0,xi*2*eta,0,0,xi*zeta,0,0,0,
170  0,0,3*eta*eta,0,0,2*eta*zeta,0,0,zeta*zeta,
171  0,0,0};
172 
173  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
174  for(int i=0; i<nBF; i++)
175  {
176  for(int j=0; j<nBF; j++)
177  {
178  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
179  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
180  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
181  }
182  }
183 }
184 
185 static void N_T_BDDF3_3D_DeriveZeta(double xi, double eta, double zeta,
186  double *values)
187 {
188  int nBF = 60; // number of basis functions
189  // monomials x-component, y-component and z-component
190  double mon_x[]={0,0,0,0,0,0,0,0,0,1,0,0,
191  0,0,0,0,0,0,xi,0,0,
192  0,0,0,eta,0,0,2*zeta,0,0,
193  0,0,0,0,0,0,xi*xi,0,0,
194  0,0,0,xi*eta,0,0,xi*2*zeta,0,0,
195  0,0,0,eta*eta,0,0,eta*2*zeta,0,0,
196  3*zeta*zeta,0,0};
197  double mon_y[]={0,0,0,0,0,0,0,0,0,0,1,0,
198  0,0,0,0,0,0,0,xi,0,
199  0,0,0,0,eta,0,0,2*zeta,0,
200  0,0,0,0,0,0,0,xi*xi,0,
201  0,0,0,0,xi*eta,0,0,xi*2*zeta,0,
202  0,0,0,0,eta*eta,0,0,eta*2*zeta,0,
203  0,3*zeta*zeta,0};
204  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,1,
205  0,0,0,0,0,0,0,0,xi,
206  0,0,0,0,0,eta,0,0,2*zeta,
207  0,0,0,0,0,0,0,0,xi*xi,
208  0,0,0,0,0,xi*eta,0,0,xi*2*zeta,
209  0,0,0,0,0,eta*eta,0,0,eta*2*zeta,
210  0,0,3*zeta*zeta};
211 
212  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
213  for(int i=0; i<nBF; i++)
214  {
215  for(int j=0; j<nBF; j++)
216  {
217  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
218  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
219  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
220  }
221  }
222 }
223 
224 static void N_T_BDDF3_3D_DeriveXiXi(double xi, double eta, double zeta,
225  double *values)
226 {
227  int nBF = 60; // number of basis functions
228  // monomials x-component, y-component and z-component
229  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
230  2,0,0,0,0,0,0,0,0,
231  0,0,0,0,0,0,0,0,0,
232  3*2*xi,0,0,2*eta,0,0,2*zeta,0,0,
233  0,0,0,0,0,0,0,0,0,
234  0,0,0,0,0,0,0,0,0,
235  0,0,0};
236  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
237  0,2,0,0,0,0,0,0,0,
238  0,0,0,0,0,0,0,0,0,
239  0,3*2*xi,0,0,2*eta,0,0,2*zeta,0,
240  0,0,0,0,0,0,0,0,0,
241  0,0,0,0,0,0,0,0,0,
242  0,0,0};
243  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
244  0,0,2,0,0,0,0,0,0,
245  0,0,0,0,0,0,0,0,0,
246  0,0,3*2*xi,0,0,2*eta,0,0,2*zeta,
247  0,0,0,0,0,0,0,0,0,
248  0,0,0,0,0,0,0,0,0,
249  0,0,0};
250 
251  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
252  for(int i=0; i<nBF; i++)
253  {
254  for(int j=0; j<nBF; j++)
255  {
256  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
257  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
258  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
259  }
260  }
261 }
262 
263 static void N_T_BDDF3_3D_DeriveXiEta(double xi, double eta, double zeta,
264  double *values)
265 {
266  int nBF = 60; // number of basis functions
267  // monomials x-component, y-component and z-component
268  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
269  0,0,0,1,0,0,0,0,0,
270  0,0,0,0,0,0,0,0,0,
271  0,0,0,2*xi,0,0,0,0,0,
272  2*eta,0,0,zeta,0,0,0,0,0,
273  0,0,0,0,0,0,0,0,0,
274  0,0,0};
275  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
276  0,0,0,0,1,0,0,0,0,
277  0,0,0,0,0,0,0,0,0,
278  0,0,0,0,2*xi,0,0,0,0,
279  0,2*eta,0,0,zeta,0,0,0,0,
280  0,0,0,0,0,0,0,0,0,
281  0,0,0};
282  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
283  0,0,0,0,0,1,0,0,0,
284  0,0,0,0,0,0,0,0,0,
285  0,0,0,0,0,2*xi,0,0,0,
286  0,0,2*eta,0,0,zeta,0,0,0,
287  0,0,0,0,0,0,0,0,0,
288  0,0,0};
289 
290  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
291  for(int i=0; i<nBF; i++)
292  {
293  for(int j=0; j<nBF; j++)
294  {
295  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
296  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
297  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
298  }
299  }
300 }
301 
302 static void N_T_BDDF3_3D_DeriveXiZeta(double xi, double eta, double zeta,
303  double *values)
304 {
305  int nBF = 60; // number of basis functions
306  // monomials x-component, y-component and z-component
307  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
308  0,0,0,0,0,0,1,0,0,
309  0,0,0,0,0,0,0,0,0,
310  0,0,0,0,0,0,2*xi,0,0,
311  0,0,0,eta,0,0,2*zeta,0,0,
312  0,0,0,0,0,0,0,0,0,
313  0,0,0};
314  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
315  0,0,0,0,0,0,0,1,0,
316  0,0,0,0,0,0,0,0,0,
317  0,0,0,0,0,0,0,2*xi,0,
318  0,0,0,0,eta,0,0,2*zeta,0,
319  0,0,0,0,0,0,0,0,0,
320  0,0,0};
321  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
322  0,0,0,0,0,0,0,0,1,
323  0,0,0,0,0,0,0,0,0,
324  0,0,0,0,0,0,0,0,2*xi,
325  0,0,0,0,0,eta,0,0,2*zeta,
326  0,0,0,0,0,0,0,0,0,
327  0,0,0};
328 
329  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
330  for(int i=0; i<nBF; i++)
331  {
332  for(int j=0; j<nBF; j++)
333  {
334  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
335  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
336  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
337  }
338  }
339 }
340 
341 static void N_T_BDDF3_3D_DeriveEtaEta(double xi, double eta, double zeta,
342  double *values)
343 {
344  int nBF = 60; // number of basis functions
345  // monomials x-component, y-component and z-component
346  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
347  0,0,0,0,0,0,0,0,0,
348  2,0,0,0,0,0,0,0,0,
349  0,0,0,0,0,0,0,0,0,
350  xi*2,0,0,0,0,0,0,0,0,
351  3*2*eta,0,0,2*zeta,0,0,0,0,0,
352  0,0,0};
353  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
354  0,0,0,0,0,0,0,0,0,
355  0,2,0,0,0,0,0,0,0,
356  0,0,0,0,0,0,0,0,0,
357  0,xi*2,0,0,0,0,0,0,0,
358  0,3*2*eta,0,0,2*zeta,0,0,0,0,
359  0,0,0};
360  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
361  0,0,0,0,0,0,0,0,0,
362  0,0,2,0,0,0,0,0,0,
363  0,0,0,0,0,0,0,0,0,
364  0,0,xi*2,0,0,0,0,0,0,
365  0,0,3*2*eta,0,0,2*zeta,0,0,0,
366  0,0,0};
367 
368  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
369  for(int i=0; i<nBF; i++)
370  {
371  for(int j=0; j<nBF; j++)
372  {
373  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
374  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
375  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
376  }
377  }
378 }
379 
380 static void N_T_BDDF3_3D_DeriveEtaZeta(double xi, double eta, double zeta,
381  double *values)
382 {
383  int nBF = 60; // number of basis functions
384  // monomials x-component, y-component and z-component
385  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
386  0,0,0,0,0,0,0,0,0,
387  0,0,0,1,0,0,0,0,0,
388  0,0,0,0,0,0,0,0,0,
389  0,0,0,xi,0,0,0,0,0,
390  0,0,0,2*eta,0,0,2*zeta,0,0,
391  0,0,0};
392  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
393  0,0,0,0,0,0,0,0,0,
394  0,0,0,0,1,0,0,0,0,
395  0,0,0,0,0,0,0,0,0,
396  0,0,0,0,xi,0,0,0,0,
397  0,0,0,0,2*eta,0,0,2*zeta,0,
398  0,0,0};
399  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
400  0,0,0,0,0,0,0,0,0,
401  0,0,0,0,0,1,0,0,0,
402  0,0,0,0,0,0,0,0,0,
403  0,0,0,0,0,xi,0,0,0,
404  0,0,0,0,0,2*eta,0,0,2*zeta,
405  0,0,0};
406 
407  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
408  for(int i=0; i<nBF; i++)
409  {
410  for(int j=0; j<nBF; j++)
411  {
412  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
413  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
414  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
415  }
416  }
417 }
418 
419 static void N_T_BDDF3_3D_DeriveZetaZeta(double xi, double eta, double zeta,
420  double *values)
421 {
422  int nBF = 60; // number of basis functions
423  // monomials x-component, y-component and z-component
424  double mon_x[]={0,0,0,0,0,0,0,0,0,0,0,0,
425  0,0,0,0,0,0,0,0,0,
426  0,0,0,0,0,0,2,0,0,
427  0,0,0,0,0,0,0,0,0,
428  0,0,0,0,0,0,xi*2,0,0,
429  0,0,0,0,0,0,eta*2,0,0,
430  3*2*zeta,0,0};
431  double mon_y[]={0,0,0,0,0,0,0,0,0,0,0,0,
432  0,0,0,0,0,0,0,0,0,
433  0,0,0,0,0,0,0,2,0,
434  0,0,0,0,0,0,0,0,0,
435  0,0,0,0,0,0,0,xi*2,0,
436  0,0,0,0,0,0,0,eta*2,0,
437  0,3*2*zeta,0};
438  double mon_z[]={0,0,0,0,0,0,0,0,0,0,0,0,
439  0,0,0,0,0,0,0,0,0,
440  0,0,0,0,0,0,0,0,2,
441  0,0,0,0,0,0,0,0,0,
442  0,0,0,0,0,0,0,0,xi*2,
443  0,0,0,0,0,0,0,0,eta*2,
444  0,0,3*2*zeta};
445 
446  memset(values, 0.0, 3*nBF*SizeOfDouble); // 3 is the space dimension
447  for(int i=0; i<nBF; i++)
448  {
449  for(int j=0; j<nBF; j++)
450  {
451  values[i ] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_x[j];
452  values[i+ nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_y[j];
453  values[i+2*nBF] += N_T_BDDF3_3D_CM[i+j*nBF]*mon_z[j];
454  }
455  }
456 }
457 
458 //polynomial degree is 3 but somehow 4 is needed in order to work properly
459 TBaseFunct3D *BF_N_T_BDDF3_3D_Obj =
460 new TBaseFunct3D(60, BF_N_T_BDDF3_3D, BFUnitTetrahedron,
461  N_T_BDDF3_3D_Funct, N_T_BDDF3_3D_DeriveXi,
462  N_T_BDDF3_3D_DeriveEta, N_T_BDDF3_3D_DeriveZeta,
463  N_T_BDDF3_3D_DeriveXiXi, N_T_BDDF3_3D_DeriveXiEta,
464  N_T_BDDF3_3D_DeriveXiZeta, N_T_BDDF3_3D_DeriveEtaEta,
465  N_T_BDDF3_3D_DeriveEtaZeta, N_T_BDDF3_3D_DeriveZetaZeta,
466  4, 1,
467  0, NULL, 3);
TBaseFunct3D
Definition: BaseFunct3D.h:27
Generated on Tue May 10 2016 09:08:00 for ParMooN by   doxygen 1.8.5