| /phasta/shapeFunction/src/ |
| H A D | Fn.c | 18 double Fn(int i, int j, double r, double s) { in Fn() argument 29 f = s-1.0/3.0; in Fn() 35 t1 = s*s; in Fn() 36 f = t1-3.0/4.0*s+3.0/28.0; in Fn() 39 f = r*s-r/4-s/4+1.0/14.0; in Fn() 46 t1 = s*s; in Fn() 47 f = t1*s-6.0/5.0*t1+2.0/5.0*s-1.0/30.0; in Fn() 50 t1 = s*s; in Fn() 51 f = r*t1-3.0/5.0*r*s+r/15-t1/5+2.0/15.0*s-1.0/60.0; in Fn() 55 f = t1*s-t1/5-3.0/5.0*r*s+2.0/15.0*r+s/15-1.0/60.0; in Fn() [all …]
|
| H A D | FnDrv.c | 18 int FnDrv(int i,int j, double r, double s, double drv[2]) { in FnDrv() argument 39 drv[1] = 2.0*s-3.0/4.0; in FnDrv() 42 drv[0] = s-1.0/4.0; in FnDrv() 50 t1 = s*s; in FnDrv() 52 drv[1] = 3.0*t1-12.0/5.0*s+2.0/5.0; in FnDrv() 55 t1 = s*s; in FnDrv() 56 drv[0] = t1-3.0/5.0*s+1.0/15.0; in FnDrv() 57 drv[1] = 2.0*r*s-3.0/5.0*r-2.0/5.0*s+2.0/15.0; in FnDrv() 61 drv[0] = 2.0*r*s-2.0/5.0*r-3.0/5.0*s+2.0/15.0; in FnDrv() 70 t1 = s*s; in FnDrv() [all …]
|
| H A D | En.c | 17 double En(int ip, double r, double s) { in En() argument 27 f = s-r; in En() 30 t1 = s*s; in En() 32 f = t1-3.0*r*s+t3; in En() 35 t1 = s*s; in En() 37 f = t1*s-6.0*r*t1+6.0*t4*s-t4*r; in En() 40 t1 = s*s; in En() 44 f = t2-10.0*r*t1*s+20.0*t5*t1-10.0*t5*r*s+t9; in En() 47 t1 = s*s; in En() 51 f = t2*s-15.0*t2*r+50.0*t5*t1*s-50.0*t5*r*t1+15.0*t10*s-t10*r; in En() [all …]
|
| H A D | Bn.c | 19 double Bn(int i, int j, int k, double r, double s, double t) { in Bn() argument 30 f = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0; in Bn() 33 f = s+2.0/5.0*r-2.0/5.0; in Bn() 40 t2 = s+r-1.0; in Bn() 45 f = (s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn() 48 t1 = s*s; in Bn() 51 f = t1+6.0/7.0*t2*s+t4/7; in Bn() 54 f = (r-2.0/9.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn() 57 f = (r-2.0/9.0)*(s+2.0/5.0*r-2.0/5.0); in Bn() 65 t3 = s+r-1.0; in Bn() [all …]
|
| H A D | EnDrv.c | 17 int EnDrv(int ip, double r, double s, double drv[2]) { in EnDrv() argument 31 drv[0] = -3.0*s+2.0*r; in EnDrv() 32 drv[1] = 2.0*s-3.0*r; in EnDrv() 35 t1 = s*s; in EnDrv() 36 t2 = r*s; in EnDrv() 42 t1 = s*s; in EnDrv() 43 t2 = t1*s; in EnDrv() 46 t5 = t4*s; in EnDrv() 52 t1 = s*s; in EnDrv() 54 t4 = r*t1*s; in EnDrv() [all …]
|
| H A D | Lagrange.c | 24 double r=par[0],s=par[1],t=par[2],t8,t4; in TetQuadLag() local 25 t8 = 1.0-r-s-t; in TetQuadLag() 27 N[1] = (2.0*s-1.0)*s; in TetQuadLag() 29 N[3] = (1.0-2.0*r-2.0*s-2.0*t)*t8; in TetQuadLag() 30 N[4] = 4.0*r*s; in TetQuadLag() 31 N[5] = 4.0*s*t; in TetQuadLag() 34 N[8] = 4.0*s*t8; in TetQuadLag() 37 t4 = -3.0+4.0*r+4.0*s+4.0*t; in TetQuadLag() 42 dN[1][1] = 4.0*s-1.0; in TetQuadLag() 50 dN[4][0] = 4.0*s; in TetQuadLag() [all …]
|
| H A D | BnDrv.c | 19 int BnDrv(int i,int j, int k, double r, double s, double t, double drv[3]) { in BnDrv() argument 46 t1 = 6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0; in BnDrv() 49 drv[2] = 2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0; in BnDrv() 52 drv[0] = 2.0/7.0*t+6.0/7.0*s+8.0/21.0*r-8.0/21.0; in BnDrv() 53 drv[1] = t+4.0/3.0*s+6.0/7.0*r-6.0/7.0; in BnDrv() 54 drv[2] = s+2.0/7.0*r-2.0/7.0; in BnDrv() 57 drv[0] = 6.0/7.0*s+2.0/7.0*r-2.0/7.0; in BnDrv() 58 drv[1] = 2.0*s+6.0/7.0*r-6.0/7.0; in BnDrv() 62 drv[0] = t+2.0/3.0*s+4.0/3.0*r-22.0/27.0; in BnDrv() 67 drv[0] = s+4.0/5.0*r-22.0/45.0; in BnDrv() [all …]
|
| H A D | blendDrv.c | 55 double r,s; in F_edgeBlendTriDrv() local 58 r = L[0] ; s = L[1] ; in F_edgeBlendTriDrv() 63 drv[0] = -2.0*s; in F_edgeBlendTriDrv() 67 drv[0] = 2.0*s; in F_edgeBlendTriDrv() 68 drv[1] = -2.0+2.0*r+4.0*s; in F_edgeBlendTriDrv() 71 drv[0] = 4.0*r-2.0+2.0*s; in F_edgeBlendTriDrv() 96 double r=L[0],s=L[1],t=L[2]; in R_edgeBlendTetDrv() local 101 drv[0] = -2.0*s; in R_edgeBlendTetDrv() 106 drv[0] = -2.0+4.0*r+2.0*s+2.0*t; in R_edgeBlendTetDrv() 113 drv[2] = -2.0*s; in R_edgeBlendTetDrv() [all …]
|
| H A D | shapeFuncInternals.h | 55 double En(int ip, double r, double s); 56 int EnDrv(int ip, double r, double s, double *drv); 57 double Fn(int i, int j, double r, double s); 58 int FnDrv(int i,int j, double r, double s, double drv[2]); 59 double Bn(int i, int j, int k, double r, double s, double t); 60 int BnDrv(int i,int j, int k, double r, double s, double t, double drv[3]);
|
| /phasta/svLS/ |
| H A D | INCOMMU.f | 50 INTEGER i, j, k, s, e, ierr, nTasks, tF, stat(MPI_STATUS_SIZE) local 68 s = cS(i)%ptr 69 e = s + cS(i)%n - 1 71 CALL MPI_IRECV(rTmp(:,s:e), cS(i)%n*dof, mpreal, i-1, 74 CALL MPI_ISEND(R(:,s:e), cS(i)%n*dof, mpreal, i-1, 85 s = cS(i)%blPtr(j) 86 e = s + cS(i)%blN(j) - 1 87 R(:,s:e) = R(:,s:e) + rTmp(:,k:k+e-s) 96 s = cS(i)%blPtr(j) 97 e = s + cS(i)%blN(j) - 1 [all …]
|
| H A D | LHS.f | 52 INTEGER i, j, k, a, Ac, ai, s, e, nTasks, tF, maxnNo, ierr, local 83 s = rowPtr(Ac) 85 DO i=s, e 93 lhs%rowPtr(1,Ac) = s 218 s = lhs%cS(i)%ptr 219 e = s + lhs%cS(i)%n - 1 221 CALL MPI_RECV(part(s:e), lhs%cS(i)%n, mpint, i-1, 225 DO j=s, e 227 IF (part(j).NE.ltg(k) .OR. j.EQ.s) THEN 239 DO j=s, e [all …]
|
| H A D | GMRES.f | 59 2 s(:), err(:) local 62 2 y(ls%sD), c(ls%sD), s(ls%sD), err(ls%sD+1)) 122 tmp = c(j)*h(j,i) + s(j)*h(j+1,i) 123 h(j+1,i) = -s(j)*h(j,i) + c(j)*h(j+1,i) 128 s(i) = h(i+1,i)/tmp 131 err(i+1) = -s(i)*err(i) 180 REAL*8, ALLOCATABLE :: u(:,:,:), h(:,:), X(:,:), y(:), c(:), s(:), local 184 2 y(ls%sD), c(ls%sD), s(ls%sD), err(ls%sD+1)) 225 tmp = c(j)*h(j,i) + s(j)*h(j+1,i) 226 h(j+1,i) = -s(j)*h(j,i) + c(j)*h(j+1,i) [all …]
|
| /phasta/phSolver/common/ |
| H A D | ValType.cc | 82 char* s = (char*) malloc(str.size()+1); in get_vector() local 83 strcpy(s,str.c_str()); in get_vector() 84 char *strTok = strtok(s," "); in get_vector() 90 free(s); in get_vector() 96 char* s = (char*) malloc(str.size()+1); in get_ivector() local 97 strcpy(s,str.c_str()); in get_ivector() 98 char *strTok = strtok(s," "); in get_ivector() 104 free(s); in get_ivector()
|
| H A D | phString.cc | 14 std::string s = ss.str(); in phstr_appendInt() local 15 strcpy(dest, s.c_str()); in phstr_appendInt() 21 std::string s = ss.str(); in phstr_appendDbl() local 22 strcpy(dest, s.c_str()); in phstr_appendDbl() 28 std::string s = ss.str(); in phstr_appendStr() local 29 strcpy(dest, s.c_str()); in phstr_appendStr()
|
| H A D | GaussLegendreSimplex.c | 9 double *r, double *s, double *t, double *J); 12 double *r, double *s,double *J); 57 double *r, double *s, double *t, double *J) { in brickToTet() argument 61 *s = 0.5e0*(1.0e0+eta)*r1; in brickToTet() 62 rs1 = 1.0e0-(*r)-(*s); in brickToTet() 67 void quadToTri(double xi,double eta,double *r, double *s, double *J) { in quadToTri() argument 71 *s = 0.5e0*(1.0e0+eta)*r1; in quadToTri()
|
| H A D | phComm.c | 11 int s; in phcomm_size() local 12 MPI_Comm_size(MPI_COMM_WORLD, &s); in phcomm_size() 13 return s; in phcomm_size()
|
| H A D | renum.f | 33 integer temp, etmp, s local 145 s = (thmax-thmin)*radcyl/ds 146 allocate (xsinfin(nfint,s+1,nsd)) 147 allocate (elcnfin(nfint,s+1)) 217 integer temp, etmp, s local 281 s = (zmax - zmin) / ds 282 allocate (xsinfin(nfint,s+1,nsd)) 283 allocate (elcnfin(nfint,s+1)) 295 imax(jj) = s
|
| H A D | ValType.h | 11 ValType(const string & s) in ValType() argument 12 : str(s) in ValType()
|
| H A D | fillsparse.f | 87 integer f, g, r, s, t local 111 s = (b-1)*nflow !starting index of the bth node's 121 lhsK(t+f,k) = lhsK(t+f,k) + EGmass(e,r+f,s+g) 149 integer f, g, r, s, t local
|
| /phasta/phSolver/compressible/ |
| H A D | get_h.c | 89 void u_line(double s, double xc[3], double udir[3], double pnt[3]); 108 double xr_xc[3], s, pt[3]; in intersect() local 150 s = dotProd(normal,xr_xc)/n_dot_u; in intersect() 152 u_line(s,xc,u,pt); /* find intersection point */ in intersect() 173 void u_line(double s, double xc[3], double udir[3], double pnt[3]){ in u_line() argument 176 pnt[i] = xc[i]+s*udir[i]; in u_line()
|
| /phasta/phSolver/testIC-DSIM/ |
| H A D | geom.sms | 162 1 24-q+s-j-g 163 1 24-r+i-k-s 175 0 34+o+w-s-x
|
| H A D | geom.sdm | 161 1 24-q+s-j-g 162 1 24-r+i-k-s
|
| /phasta/phSolver/common/test/ |
| H A D | phIOreadIlwork.cc | 99 std::string s = ss.str(); in main() local 100 int size = s.size(); in main() 104 int ret = MPI_File_write_at(outfile,offset,(void*)s.c_str(),s.size(),MPI_BYTE,&status); in main()
|
| /phasta/AcuStat/src/ |
| H A D | phasta.cc | 111 std::string s = exec_string.str(); in phasta() local 112 system( s.c_str() ); in phasta()
|
| /phasta/M2NFixBnd/src/ |
| H A D | phasta.cc | 110 std::string s = exec_string.str(); in phasta() local 111 system( s.c_str() ); in phasta()
|