# Diff of /trunk/finley/src/Quadrature.c

revision 1169 by gross, Thu Mar 29 06:44:30 2007 UTC revision 1170 by btully, Fri May 25 05:24:57 2007 UTC
36    double Q1,Q2;    double Q1,Q2;
37    #define DIM 2    #define DIM 2
38
39    /*  the easy case: */    /*  the easy cases: */
40
41    if (numQuadNodes==1) {    if (numQuadNodes==1) {
45      } else if (numQuadNodes==3){
55      } else if (numQuadNodes==4){
68    } else {    } else {
69
70      /*  get scheme on [0.1]^2 */      /*  get scheme on [0.1]^2 */

72      if (! Finley_noError()) return;      if (! Finley_noError()) return;
73
94    double Q1,Q2,Q3,JA11,JA12,JA13,JA21,JA22,JA23,JA31,JA32,JA33,DET;    double Q1,Q2,Q3,JA11,JA12,JA13,JA21,JA22,JA23,JA31,JA32,JA33,DET;
95    #define DIM 3    #define DIM 3
96
97    /*  the easy case: */    /*  the easy cases: */

98    if (numQuadNodes==1) {    if (numQuadNodes==1) {
103      } else if (numQuadNodes==4){
104        double alpha=0.58541020;
105        double beta =0.13819660;
122      } else if (numQuadNodes==5){
143    } else {    } else {
144
145      /*  get scheme on [0.1]^3 */      /*  get scheme on [0.1]^3 */
# Line 459  int Finley_Quad_getNumNodesTri(int order Line 520  int Finley_Quad_getNumNodesTri(int order
521    if (order<=1) {    if (order<=1) {
522        return 1;        return 1;
523      } else if (order==2){
524          return 3;
525      } else if (order==3){
526          return 4;
527    } else {    } else {
529        if (Finley_noError()) {        if (Finley_noError()) {
# Line 483  int Finley_Quad_getNumNodesTet(int order Line 548  int Finley_Quad_getNumNodesTet(int order
549    if (order<=1) {    if (order<=1) {
550        return 1;        return 1;
551      } else if (order==2){
552          return 4;
553      } else if (order==3){
554          return 5;
555    } else {    } else {