inf_fe_lagrange_eval.c

一个用来实现偏微分方程中网格的计算库

        case 6:
          return 134521.003125*(v+1.)*v*(v+.8888888888888889)*(v+.7777777777777778)*(v+.6666666666666667)*(v+.5555555555555556)*(v+.3333333333333333)*(v+.2222222222222222)*(v+.1111111111111111);

        case 7:
          return -89680.66875000004*(v+1.)*v*(v+.8888888888888889)*(v+.7777777777777778)*(v+.6666666666666667)*(v+.5555555555555556)*(v+.4444444444444444)*(v+.2222222222222222)*(v+.1111111111111111);

        case 8:
          return 38434.57232142857*(v+1.)*v*(v+.8888888888888889)*(v+.7777777777777778)*(v+.6666666666666667)*(v+.5555555555555556)*(v+.4444444444444444)*(v+.3333333333333333)*(v+.1111111111111111);

        case 9:
          return -9608.643080357148*(v+1.)*v*(v+.8888888888888889)*(v+.7777777777777778)*(v+.6666666666666667)*(v+.5555555555555556)*(v+.4444444444444444)*(v+.3333333333333333)*(v+.2222222222222222);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    case TENTH:
      switch (i)
        {
        case 0:
          return 2755.73192239859*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 1:
          return 2755.73192239859*(v+1.)*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 2:
          return -27557.31922398588*(v+1.)*v*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 3:
          return 124007.9365079364*(v+1.)*v*(v+.9)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 4:
          return -330687.8306878306*(v+1.)*v*(v+.9)*(v+.8)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 5:
          return 578703.7037037039*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.5)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 6:
          return -694444.4444444449*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.4)*(v+.3)*(v+.2)*(v+.1);

        case 7:
          return 578703.7037037039*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.3)*(v+.2)*(v+.1);

        case 8:
          return -330687.8306878306*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.2)*(v+.1);

        case 9:
          return 124007.9365079364*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.1);

        case 10:
          return -27557.31922398588*(v+1.)*v*(v+.9)*(v+.8)*(v+.7)*(v+.6)*(v+.5)*(v+.4)*(v+.3)*(v+.2);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    case ELEVENTH:
      switch (i)
        {
        case 0:
          return -7147.658895778214*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 1:
          return 7147.65889577822*(v+1.)*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 2:
          return 78624.2478535604*(v+1.)*v*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 3:
          return -393121.2392678026*(v+1.)*v*(v+.9090909090909091)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 4:
          return 1179363.717803407*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 5:
          return -2358727.435606812*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 6:
          return 3302218.409849539*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 7:
          return -3302218.40984954*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 8:
          return 2358727.435606811*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.2727272727272727)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 9:
          return -1179363.717803406*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.1818181818181818)*(v+.9090909090909091e-1);

        case 10:
          return 393121.2392678019*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.9090909090909091e-1);

        case 11:
          return -78624.24785356052*(v+1.)*v*(v+.9090909090909091)*(v+.8181818181818182)*(v+.7272727272727273)*(v+.6363636363636364)*(v+.5454545454545455)*(v+.4545454545454545)*(v+.3636363636363636)*(v+.2727272727272727)*(v+.1818181818181818);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    case TWELFTH:
      switch (i)
        {
        case 0:
          return 18613.92623376623*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 1:
          return 18613.92623376623*(v+1.)*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 2:
          return -223367.1148051951*(v+1.)*v*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 3:
          return 1228519.131428571*(v+1.)*v*(v+.9166666666666667)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 4:
          return -4095063.771428572*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 5:
          return 9213893.485714271*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 6:
          return -14742229.57714284*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 7:
          return 17199267.84*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 8:
          return -14742229.57714286*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 9:
          return 9213893.485714291*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.25)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 10:
          return -4095063.77142857*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.1666666666666667)*(v+.8333333333333333e-1);

        case 11:
          return 1228519.131428571*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.8333333333333333e-1);

        case 12:
          return -223367.1148051951*(v+1.)*v*(v+.9166666666666667)*(v+.8333333333333333)*(v+.75)*(v+.6666666666666667)*(v+.5833333333333333)*(v+.5)*(v+.4166666666666667)*(v+.3333333333333333)*(v+.25)*(v+.1666666666666667);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    case THIRTEENTH:
      switch (i)
        {
        case 0:
          return -48638.84613847011*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 1:
          return 48638.84613847013*(v+1.)*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 2:
          return 632304.9998001099*(v+1.)*v*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 3:
          return -3793829.998800662*(v+1.)*v*(v+.9230769230769231)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 4:
          return 13910709.99560246*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 5:
          return -34776774.98900616*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 6:
          return 62598194.98021102*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 7:
          return -83464259.97361468*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 8:
          return 83464259.97361466*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 9:
          return -62598194.98021099*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 10:
          return 34776774.98900612*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.2307692307692308)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 11:
          return -13910709.99560244*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.1538461538461538)*(v+.7692307692307692e-1);

        case 12:
          return 3793829.998800668*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.7692307692307692e-1);

        case 13:
          return -632304.9998001105*(v+1.)*v*(v+.9230769230769231)*(v+.8461538461538462)*(v+.7692307692307692)*(v+.6923076923076923)*(v+.6153846153846154)*(v+.5384615384615385)*(v+.4615384615384615)*(v+.3846153846153846)*(v+.3076923076923077)*(v+.2307692307692308)*(v+.1538461538461538);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    case FOURTEENTH:
      switch (i)
        {
        case 0:
          return 127463.0033762123*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 1:
          return 127463.0033762123*(v+1.)*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 2:
          return -1784482.047266971*(v+1.)*v*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 3:
          return 11599133.30723529*(v+1.)*v*(v+.9285714285714286)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 4:
          return -46396533.2289412*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 5:
          return 127590466.3795883*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 6:
          return -255180932.7591769*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 7:
          return 382771399.1387658*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 8:
          return -437453027.5871608*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 9:
          return 382771399.1387654*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 10:
          return -255180932.7591768*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 11:
          return 127590466.3795886*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2142857142857143)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 12:
          return -46396533.22894132*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.1428571428571429)*(v+.7142857142857143e-1);

        case 13:
          return 11599133.30723532*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.7142857142857143e-1);

        case 14:
          return -1784482.047266969*(v+1.)*v*(v+.9285714285714286)*(v+.8571428571428571)*(v+.7857142857142857)*(v+.7142857142857143)*(v+.6428571428571429)*(v+.5714285714285714)*(v+.5)*(v+.4285714285714286)*(v+.3571428571428571)*(v+.2857142857142857)*(v+.2142857142857143)*(v+.1428571428571429);

        default:
          std::cerr << "bad index i = " << i << std::endl;
          libmesh_error();
        }


    // provided the constructor worked fine, we should never end up here.
    default:
      std::cerr << "Lagrange polynomials only defined up to 15. "
		<< std::endl;
      libmesh_error();
    }

  // we never end up here.
  return 0.;
}