1 #include <FCMangle.h> 2 #define symtri FortranCInterface_GLOBAL_(symtri,SYMTRI) 3 4 typedef double DARR3[3]; 5 6 int triIntPnt(int, DARR3**,double**); 7 8 void symtri(int *n1, double pt[][4], double wt[], int *err) 9 { 10 double *lwt; 11 DARR3 *lpt; 12 int i,j; 13 *err = triIntPnt(*n1, &lpt, &lwt); 14 for(i=0; i < *n1; i++){ 15 wt[i] = lwt[i]; 16 for(j=0; j < 3; j++) 17 pt[i][j] = lpt[i][j]; 18 } 19 } 20 21 #include <stdio.h> 22 23 #define A 0.3333333333333330000000000000000000 24 #define B 0.666666666666667 25 #define C 0.166666666666667 26 #define D 0.562500000000000 27 #define E 0.520833333333333 28 29 #define F 0.109951743655322 30 #define G 0.223381589678011 31 #define H 0.816847572980459 32 #define I 0.091576213509771 33 #define J 0.108103018168070 34 #define K 0.445948490915965 35 36 #define L 0.225000000000000 37 #define M 0.125939180544827 38 #define N 0.132394152788506 39 40 #define O 0.797426985353087 41 #define P 0.101286507323456 42 #define Q 0.470142064105115 43 #define R 0.059715871789770 44 45 #define S 0.050844906370207 46 #define T 0.116786275726379 47 #define U 0.082851075618374 48 49 #define V 0.873821971016996 50 #define W 0.063089014491502 51 #define X 0.501426509658179 52 #define Y 0.249286745170910 53 #define Z 0.636502499121399 54 #define AA 0.310352451033785 55 #define BB 0.053145049844816 56 57 /* typedef double DARR3[3] ; */ 58 59 static double rst1[][3] = {{A,A,A}}; 60 static double wt1[] = {1.00000000000000000}; 61 62 static double rst3[][3] = {{B,C,C},{C,B,C},{C,C,B}}; 63 static double wt3[] = {A,A,A}; 64 65 static double rst4[][3] = {{A,A,A}, 66 {0.6,0.2,0.2},{0.2,0.6,0.2},{0.2,0.2,0.6}}; 67 static double wt4[] = {-D,E,E,E}; 68 69 static double rst6[][3] = {{H,I,I},{I,H,I},{I,I,H}, 70 {J,K,K},{K,J,K},{K,K,J}}; 71 static double wt6[] = {F,F,F,G,G,G}; 72 73 static double rst7[][3] = {{A,A,A}, 74 {O,P,P},{P,O,P},{P,P,O}, 75 {Q,Q,R},{Q,R,Q},{R,Q,Q}}; 76 static double wt7[] = {L,M,M,M,N,N,N}; 77 78 static double rst12[][3] = {{V,W,1.0-V-W},{W,V,1.0-V-W},{W,W,1.0-W-W}, 79 {X,Y,1.0-X-Y},{Y,X,1.0-X-Y},{Y,Y,1.0-Y-Y}, 80 {Z,AA,1.0-Z-AA},{AA,Z,1.0-Z-AA},{Z,BB,1.0-Z-BB}, 81 {AA,BB,1.0-AA-BB},{BB,AA,1.0-BB-AA},{BB,Z,1.0-BB-Z}}; 82 83 static double wt12[] = {S,S,S,T,T,T,U,U,U,U,U,U}; 84 85 int triIntPnt(int nint, DARR3 **bcord, double **wt) 86 { 87 int retval = 1 ; 88 if( nint == 3 ) {*bcord = rst3 ; *wt = wt3; } 89 else if( nint == 1 ){*bcord = rst1 ; *wt = wt1; } 90 else if( nint == 4 ){*bcord = rst4 ; *wt = wt4; } 91 else if( nint == 6 ){*bcord = rst6 ; *wt = wt6; } 92 else if( nint == 7 ){*bcord = rst7 ; *wt = wt7; } 93 else if( nint == 12){*bcord = rst12; *wt = wt12;} 94 else 95 { 96 fprintf(stderr,"\n%d integration points unsupported in symtri.c; give {1,3,4,6,7,12}\n",nint); 97 retval = 0; 98 } 99 return retval ; 100 } 101 102 103 104