| /phasta/shapeFunction/src/ |
| 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() 39 t1 = t*t; in Bn() 42 f = t1+6.0/5.0*t*t2+3.0/10.0*t4; 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() 54 f = (r-2.0/9.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn() 64 t1 = t*t; in Bn() 67 f = t1*t+12.0/7.0*t3*t1+6.0/7.0*t*t5+4.0/35.0*t3*t5; in Bn() 70 t2 = t*t; in Bn() 73 f = (s+2.0/9.0*r-2.0/9.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5); in Bn() [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() 28 N[2] = (2.0*t-1.0)*t; in TetQuadLag() 29 N[3] = (1.0-2.0*r-2.0*s-2.0*t)*t8; in TetQuadLag() 31 N[5] = 4.0*s*t; in TetQuadLag() 32 N[6] = 4.0*t*r; in TetQuadLag() 35 N[9] = 4.0*t*t8; in TetQuadLag() 37 t4 = -3.0+4.0*r+4.0*s+4.0*t; in TetQuadLag() 46 dN[2][2] = 4.0*t-1.0; in TetQuadLag() 54 dN[5][1] = 4.0*t; 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() 62 drv[0] = t+2.0/3.0*s+4.0/3.0*r-22.0/27.0; in BnDrv() 77 t1 = t*t; in BnDrv() 79 t3 = t*t2; in BnDrv() 87 t1 = t*t; in BnDrv() 89 t3 = t*t2; in BnDrv() [all …]
|
| H A D | blendDrv.c | 96 double r=L[0],s=L[1],t=L[2]; in R_edgeBlendTetDrv() local 106 drv[0] = -2.0+4.0*r+2.0*s+2.0*t; in R_edgeBlendTetDrv() 112 drv[1] = -2.0*t; in R_edgeBlendTetDrv() 117 drv[1] = -2.0+2.0*r+4.0*s+2.0*t; in R_edgeBlendTetDrv() 121 drv[0] = -2.0*t; in R_edgeBlendTetDrv() 126 drv[0] = 2.0*t; in R_edgeBlendTetDrv() 127 drv[1] = 2.0*t; in R_edgeBlendTetDrv() 128 drv[2] = -2.0+2.0*r+2.0*s+4.0*t; in R_edgeBlendTetDrv()
|
| H A D | shapeFuncInternals.h | 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]);
|
| H A D | modeShapeDrv.c | 30 double rs,rst2,P1P2,t,P1,P2,dP1,dP2,t1,t2; in F_modeShapeTriDrv() local
|
| /phasta/phSolver/compressible/ |
| H A D | e3eig1.f | 268 dimension offd(npro), t(npro), 291 t = pt5 * (rlam(:,1) - rlam(:,5)) / offd 292 t = sign(one, t) / ( abs(t) + sqrt(one + t**2) ) 296 Rcs = one / sqrt( one + t**2 ) 297 Rsn = t * Rcs 299 rlam(:,1) = rlam(:,1) - t * offd 300 rlam(:,5) = rlam(:,5) + t * offd 304 t = Rcs * Q(:,1,1) - Rsn * Q(:,1,5) 306 Q(:,1,1) = t 308 t = Rcs * Q(:,2,1) - Rsn * Q(:,2,5) [all …]
|
| H A D | secs.f | 1 real*8 function secs(t) 3 real t local
|
| /phasta/phSolver/common/ |
| H A D | blowerControl.f90 | 97 real*8 :: t local 134 t = mod(BC_t - blower(iBlower)%t0, t4) 136 if(t0 <= t .and. t < t1) then !rising edge 137 vmagref = (vmax - vmin)*(t - t0)/(t1 - t0) + vmin 138 elseif(t1 <= t .and. t < t2) then !high 140 elseif(t2 <= t .and. t < t3) then !falling edge 141 vmagref = (vmax - vmin)*(t - t3)/(t2 - t3) + vmin 142 elseif(t3 <= t .and. t < t4) then !low 154 t = (BC_t - blower(iBlower)%t0)/blower(iBlower)%t_cycle 155 vmagref = 0.5*(vmax - vmin)*sin(2*pi*t) + 0.5*(vmax + vmin)
|
| H A D | fillsparse.f | 87 integer f, g, r, s, t local 119 t = (g-1)*nflow 121 lhsK(t+f,k) = lhsK(t+f,k) + EGmass(e,r+f,s+g) 149 integer f, g, r, s, t local
|
| H A D | GaussLegendreSimplex.c | 9 double *r, double *s, double *t, double *J); 57 double *r, double *s, double *t, double *J) { in brickToTet() argument 63 *t = 0.5e0*(1.0e0+zeta)*rs1; in brickToTet()
|
| H A D | phIO.cc | 24 inline void printTime(const char* key, double t) { in printTime() argument 27 fprintf(stderr, "%s %f\n", key, t); in printTime()
|
| H A D | bctint.f | 303 real*8 BC(nshg,ndofBC), timel,t local 319 t=timel-nper*BCperiod(i) ! now time in periodic domain 323 if(BCtptr(i)%p(j,4).gt.t) then 326 wr=(t-BCtptr(i)%p(j-1,4))/(BCtptr(i)%p(j,4)-BCtptr(i)%p(j-1,4))
|
| /phasta/phSolver/testIC-DSIM/ |
| H A D | geom.sms | 164 2 24-1+2+t-d 165 2 24-t+3+6-e 172 0 34+m+x-p-t
|
| H A D | geom.sdm | 163 2 24-1+2+t-d 164 2 24-t+3+6-e
|
| /phasta/phSolver/AMG/ |
| H A D | amgread.txt | 25 statements: 1) setup AMG by calling ramg_control. ( don't worry you'll 136 coarse level vector. We don't need to do this every iteration. This can 147 the math hasn't changed. But here if you combine matrix-matrix product 176 3) Eventually you want your own eigen calculator, can't use ARPACK 185 When you do sparse matrix-matrix product, don't attempt to do it in one 198 activated. Don't worry about that though. 204 nodes, though I don't think there will be substantial difference. 222 which won't be a big issue for the overall performance. It is a direct 225 don't worry about it. (slightly different value of G_bb only). 248 coarse first fine second (or other way around). In direct solve we can't
|
| H A D | ramg_coarse.f | 473 integer :: i,j,k,t local 483 t = heap(j) 485 heap(k) = t 519 integer i,j,k,t local 524 t = heap(i) 526 heap(j) = t 543 t = heap(i) 545 heap(j) = t
|
| H A D | ramg_tools.f | 716 integer i,t variable 722 read(264,*)t,iarray(i)
|
| /phasta/CMakeFiles/ |
| H A D | FindPETSC.cmake | 16 message(WARNING "pkg-config with --static didn't work")
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 92 c$$$ t131 = -8.D0*t57*t18-20.D0*t6*t7+56.D0*t77*t7-42.D0*t12*t8-24.D0*t 269 c t | i sigma i| i sigma| i| 533 cad d_{,t} + sign(phi)*( d_{,i}/|d_{,i}| )* d_{,i} = sign(phi) 537 cad d_{,t} + u * d_{,i} = sign(phi)
|
| H A D | e3source.fRHSonly | 92 c$$$ t131 = -8.D0*t57*t18-20.D0*t6*t7+56.D0*t77*t7-42.D0*t12*t8-24.D0*t 269 c t | i sigma i| i sigma| i| 519 cad d_{,t} + sign(phi)*( d_{,i}/|d_{,i}| )* d_{,i} = sign(phi) 523 cad d_{,t} + u * d_{,i} = sign(phi)
|
| /phasta/ |
| H A D | CMakeLists.txt | 94 option(PHASTA_SKIP_MPI_CHECK "Bypass CMake MPI code. Don't use this" OFF)
|
| /phasta/converterIO/ |
| H A D | README.md | 54 … count mesh (up to a few thousands is ok) and not on multiple million part meshes (don't do that!).
|
| /phasta/phSolver/ |
| H A D | CHANGES | 259 7) (IC) fixed a bug in the lhs for the scalar solver. It was not accounting for the phi_{,t} term …
|