/dports/games/powder-toy/The-Powder-Toy-95.0/src/simulation/elements/ |
H A D | DMG.cpp | 53 int r, rr, rx, ry, nxi, nxj, t, dist; in update() local 68 for (nxi=-rad; nxi<=rad; nxi++) in update() 69 if (x+nxi>=0 && y+nxj>=0 && x+nxi<XRES && y+nxj<YRES && (nxi || nxj)) in update() 74 rr = pmap[y+nxj][x+nxi]; in update() 77 angle = atan2((float)nxj, nxi); in update() 82 sim->vx[(y+nxj)/CELL][(x+nxi)/CELL] += fx; in update() 83 sim->vy[(y+nxj)/CELL][(x+nxi)/CELL] += fy; in update() 84 sim->pv[(y+nxj)/CELL][(x+nxi)/CELL] += 1.0f; in update() 89 sim->part_change_type(ID(rr), x+nxi, y+nxj, PT_BRMT); in update() 91 sim->part_change_type(ID(rr), x+nxi, y+nxj, PT_BGLA); in update() [all …]
|
H A D | BOMB.cpp | 66 int nxi, nxj; in update() local 69 for (nxi=-rad; nxi<=rad; nxi++) in update() 70 if ((pow((float)nxi,2))/(pow((float)rad,2))+(pow((float)nxj,2))/(pow((float)rad,2))<=1) in update() 72 int ynxj = y + nxj, xnxi = x + nxi; in update() 93 for (nxi=-(rad+1); nxi<=(rad+1); nxi++) in update() 94 …if ((pow((float)nxi,2))/(pow((float)(rad+1),2))+(pow((float)nxj,2))/(pow((float)(rad+1),2))<=1 && … in update() 96 nb = sim->create_part(-3, x+nxi, y+nxj, PT_EMBR); in update()
|
H A D | FRAY.cpp | 56 int r, nxx, nyy, len, nxi, nyi, rx, ry; in update() local 65 for (nxx = 0, nyy = 0, nxi = rx*-1, nyi = ry*-1, len = 0; ; nyy+=nyi, nxx+=nxi, len++) { in update() 66 if (!(x+nxi+nxx<XRES && y+nyi+nyy<YRES && x+nxi+nxx >= 0 && y+nyi+nyy >= 0) || len>curlen) { in update() 69 r = pmap[y+nyi+nyy][x+nxi+nxx]; in update() 71 r = sim->photons[y+nyi+nyy][x+nxi+nxx]; in update() 73 parts[ID(r)].vx += nxi*((parts[i].temp-273.15)/10.0f); in update()
|
H A D | CRAY.cpp | 53 int nxx, nyy, docontinue, nxi, nyi; in update() local 91 …for (docontinue = 1, nxi = rx*-1, nyi = ry*-1, nxx = spacesRemaining*nxi, nyy = spacesRemaining*ny… in update() 93 if (!(x+nxi+nxx<XRES && y+nyi+nyy<YRES && x+nxi+nxx >= 0 && y+nyi+nyy >= 0)) { in update() 96 r = pmap[y+nyi+nyy][x+nxi+nxx]; in update() 97 …if (!sim->IsWallBlocking(x+nxi+nxx, y+nyi+nyy, TYP(parts[i].ctype)) && (!sim->pmap[y+nyi+nyy][x+nx… in update() 98 … int nr = sim->create_part(-1, x+nxi+nxx, y+nyi+nyy, TYP(parts[i].ctype), ID(parts[i].ctype)); in update()
|
/dports/math/R-cran-gss/gss/src/ |
H A D | reg.f | 10 nn = nnull + nxi 80 23030 if(.not.(i.le.nxi))goto 23032 87 call dcopy (nxi, q, nxi+1, wk, 1) 89 23033 if(.not.(i.le.nxi))goto 23035 91 23036 if(.not.(j.le.nxi))goto 23038 101 23039 if(.not.(i.le.nxi))goto 23041 103 23042 if(.not.(j.le.nxi))goto 23044 113 call dcopy (nxi, wk, 1, q, nxi+1) 114 call dsyev ('n', 'u', nxi, sr, nobs, mu, wk, 3*nxi, info) 118 trc = trc + dlog (mu(nxi-i+1)) [all …]
|
H A D | hzdnewton10.f | 2 subroutine hzdnewton10 (cd, nxis, q, nxi, rs, nt, nobs, cntsum, cn 4 integer nxis, nxi, nt, nobs, cntsum, cnt(*), maxiter, jpvt(*), inf 6 double precision cd(*), q(nxi,*), rs(nt,*), intrs(*), rho(*), prec 43 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 44 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 61 if(j.le.nxi)then 71 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 116 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) 117 lkhdnew = lkhdnew + ddot (nxi, cdnew, 1, wk, 1) / 2.d0 226 integer nxis, nxi, nt, jpvt(*) [all …]
|
H A D | cdennewton10.f | 2 subroutine cdennewton10 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, 4 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 5 double precision cd(*), q(nxi,*), rs(nobs,*), intrs(*), prec, mchp 22 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 41 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 42 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 66 if(j.le.nxi)then 76 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 121 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) 122 lkhdnew = lkhdnew + ddot (nxi, cdnew, 1, wk, 1) / 2.d0 [all …]
|
H A D | dnewton10.f | 2 subroutine dnewton10 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, int 4 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 5 double precision cd(*), q(nxi,*), rs(nobs,*), intrs(*), prec, mchp 22 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 47 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 48 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 71 if(j.le.nxi)then 81 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 132 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) 133 lkhdnew = lkhdnew + ddot (nxi, cdnew, 1, wk, 1) / 2.d0 [all …]
|
H A D | hzdnewton.f | 2 subroutine hzdnewton (cd, nxis, q, nxi, rs, nt, nobs, cntsum, cnt, 4 integer nxis, nxi, nt, nobs, cntsum, cnt(*), nqd, nx, maxiter, jpv 6 double precision cd(*), q(nxi,*), rs(nxis,*), qdrs(nqd,nxis,*), qd 21 call hzdnewton1 (cd, nxis, q, nxi, rs, nt, nobs, cntsum, cnt, qdrs 81 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 82 lkhd = ddot (nxi, cd, 1, wk, 1) / 2.d0 - fitmean 115 23034 if(.not.(i.le.nxi))goto 23036 117 23037 if(.not.(j.le.nxi))goto 23039 126 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 192 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) [all …]
|
H A D | dnewton.f | 2 subroutine dnewton (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, qdrs, 4 integer nxis, nxi, nobs, cntsum, cnt(*), nqd, nt, maxiter, jpvt(*) 6 double precision cd(*), q(nxi,*), rs(nxis,*), qdrs(nqd,*), bwt(*), 23 call dnewton1 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, qdrs, nqd, 108 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 109 lkhd = ddot (nxi, cd, 1, wk, 1) / 2.d0 - fitmean + norm 148 23048 if(.not.(i.le.nxi))goto 23050 150 23051 if(.not.(j.le.nxi))goto 23053 159 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 241 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) [all …]
|
/dports/math/arpack-ng/arpack-ng-3.8.0/EXAMPLES/BAND/ |
H A D | znbdr2.f | 131 nxi = 10 132 n = nxi*nxi 196 kl = nxi 197 ku = nxi 203 h = one / dcmplx (nxi+1) 218 do 50 i = 1, nxi 219 lo = (i-1)*nxi 220 do 40 j = lo+1, lo+nxi-1 233 do 80 i = 1, nxi-1 234 lo = (i-1)*nxi [all …]
|
H A D | cnbdr2.f | 131 nxi = 10 132 n = nxi*nxi 196 kl = nxi 197 ku = nxi 203 h = one / cmplx(nxi+1) 218 do 50 i = 1, nxi 219 lo = (i-1)*nxi 220 do 40 j = lo+1, lo+nxi-1 233 do 80 i = 1, nxi-1 234 lo = (i-1)*nxi [all …]
|
/dports/math/R-cran-gss/gss/src/ratfor/ |
H A D | reg.r | 10 mchpr, v(nnull+nxi,*), mu(*), wk(*) 17 nn = nnull + nxi 65 for (i=1;i<=nxi;i=i+1) { 70 call dcopy (nxi, q, nxi+1, wk, 1) 71 for (i=1;i<=nxi;i=i+1) { 72 for (j=i;j<=nxi;j=j+1) 76 for (i=1;i<=nxi;i=i+1) { 77 for (j=i;j<=nxi;j=j+1) { 83 call dcopy (nxi, wk, 1, q, nxi+1) 85 call dsyev ('n', 'u', nxi, sr, nobs, mu, wk, 3*nxi, info) [all …]
|
H A D | hzdnewton10.r | 6 subroutine hzdnewton10 (cd, nxis, q, nxi, rs, nt, nobs, cntsum, cnt, 9 integer nxis, nxi, nt, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 37 integer nxis, nxi, nt, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 53 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 54 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 67 if (j<=nxi) v(i,j) = v(i,j) + q(i,j) 71 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 101 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) 102 lkhdnew = lkhdnew + ddot (nxi, cdnew, 1, wk, 1) / 2.d0 189 integer nxis, nxi, nt, jpvt(*) [all …]
|
H A D | dnewton10.r | 6 subroutine dnewton10 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, 9 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 21 call dnewton101 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, intrs, 33 subroutine dnewton101 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, 37 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 57 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 58 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 73 if (j<=nxi) v(i,j) = v(i,j) + q(i,j) 77 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 110 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) [all …]
|
H A D | cdennewton10.r | 6 subroutine cdennewton10 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, intrs, 9 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 21 call cdennewton101 (cd, nxis, q, nxi, rs, nobs, cntsum, cnt, intrs, 37 integer nxis, nxi, nobs, cntsum, cnt(*), maxiter, jpvt(*), info 53 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 54 lkhd = lkhd + ddot (nxi, cd, 1, wk, 1) / 2.d0 69 if (j<=nxi) v(i,j) = v(i,j) + q(i,j) 73 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 103 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) 104 lkhdnew = lkhdnew + ddot (nxi, cdnew, 1, wk, 1) / 2.d0 [all …]
|
H A D | hzdnewton.r | 6 subroutine hzdnewton (cd, nxis, q, nxi, rs, nt, nobs, cntsum, cnt, qdrs, nqd, 9 integer nxis, nxi, nt, nobs, cntsum, cnt(*), nqd, nx, maxiter, jpvt(*), info 44 integer nxis, nxi, nt, nobs, cntsum, cnt(*), nqd, nx, maxiter, jpvt(*), info 45 double precision cd(*), q(nxi,*), rs(nxis,*), qdrs(nqd,nxis,*), qdwt(nqd,*), 75 call dsymv ('u', nxi, 1.d0, q, nxi, cd, 1, 0.d0, wk, 1) 76 lkhd = ddot (nxi, cd, 1, wk, 1) / 2.d0 - fitmean 97 for (i=1;i<=nxi;i=i+1) { 98 for (j=i;j<=nxi;j=j+1) v(i,j) = v(i,j) + q(i,j) 101 call dsymv ('u', nxi, -1.d0, q, nxi, cd, 1, 1.d0, mu, 1) 144 call dsymv ('u', nxi, 1.d0, q, nxi, cdnew, 1, 0.d0, wk, 1) [all …]
|
/dports/science/afni/afni-AFNI_21.3.16/src/ |
H A D | 1dfft.c | 17 int ii , jj , nx,nfft=0,ny , nopt,nby2 , ignore=0,nxi , use=0 ; in main() local 156 nxi = nx - ignore ; in main() 158 if( use > nxi ){ in main() 164 nxi = use ; in main() 165 } else if( nxi < UMIN ){ in main() 186 if( nx != nxi ){ in main() 203 if( nfft < nxi ) nfft = nxi ; in main() 221 THD_quadratic_detrend( nxi , far+ignore , &f0,&f1,&f2 ) ; in main() 225 f0,f1,f2 , f0+0.5*nxi*f1+0.25*nxi*nxi*f2 , in main() 226 f0+(nxi-1)*f1+(nxi-1)*(nxi-1)*f2 ) ; in main() [all …]
|
H A D | mri_intedge.c | 38 int xm,xp,ym,yp,zm,zp , nxi,nyi,nzi , nxyi,nxyzi ; in mri_interior_edgeize() local 50 nxi = inim->nx ; if( nxi < 32 ) RETURN( -2 ) ; in mri_interior_edgeize() 53 nxyi = nxi*nyi ; nxyzi = nxyi*nzi ; in mri_interior_edgeize() 58 if( xm > 1 || ym > 1 || zm > 1 || xp < nxi-2 || yp < nyi-2 || zp <= nzi-2 ){ in mri_interior_edgeize() 59 int xc = nxi - (xm + ((nxi-1)-xp)) ; in mri_interior_edgeize() 64 cim = mri_zeropad_3D( -xm , xp-(nxi-1) , in mri_interior_edgeize() 89 … ININFO_message(" mri_interior_edgeize: input=%dx%dx%d cropped=%dx%dx%d",nxi,nyi,nzi,nx,ny,nz) ; in mri_interior_edgeize() 211 pp = ii + jj*nxi + kk*nxyi ; /* 1D in inar */ in mri_interior_edgeize()
|
/dports/cad/tochnog/Sources/Tochnog-Latest-jan-2014/ |
H A D | node.cc | 72 nxi = neta = nzeta = 1; in nod_nod() 74 nxi = 2; in nod_nod() 77 nxi = 3; in nod_nod() 80 nxi = 4; in nod_nod() 83 nxi = 2; in nod_nod() 87 nxi = 3; in nod_nod() 91 nxi = 4; in nod_nod() 95 nxi = 2; in nod_nod() 100 nxi = 3; in nod_nod() 105 nxi = 4; in nod_nod() [all …]
|
/dports/math/R-cran-gss/gss/R/ |
H A D | sshzd.R | 279 q.wk0[1:nxi,1:nxi] <- 10^(lambda[1])*q.wk 280 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 334 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 335 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 457 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 458 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 568 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 569 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 640 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 641 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- [all …]
|
H A D | ssllrm.R | 243 q.wk0[1:nxi,1:nxi] <- 10^(lambda[1])*q.wk 244 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 298 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 299 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 443 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 444 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 578 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 579 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 653 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 654 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- [all …]
|
H A D | gssanova.R | 139 nxi <- ncol(r) functionVar 162 q.wk[1:nxi,1:nxi] <- 10^(la.wk[1]+theta)*q 163 q.wk[(nxi+1):nxiz,(nxi+1):nxiz] <- 256 q.wk[1:nxi,1:nxi] <- 10^theta*q 257 q.wk[(nxi+1):nxiz,(nxi+1):nxiz] <- 274 nxi <- ncol(r) functionVar 316 q.wk[1:nxi,1:nxi] <- 10^nlambda*qq.wk 317 q.wk[(nxi+1):nxiz,(nxi+1):nxiz] <- 415 q.wk[1:nxi,1:nxi] <- qq.wk 416 q.wk[(nxi+1):nxiz,(nxi+1):nxiz] <- [all …]
|
H A D | sshzd1.R | 295 q.wk0[1:nxi,1:nxi] <- 10^(lambda[1])*q.wk 296 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 347 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 348 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 456 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 457 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 562 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 563 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- 631 q.wk0[1:nxi,1:nxi] <- 10^(lambda)*q.wk 632 q.wk0[(nxi+1):nxiz,(nxi+1):nxiz] <- [all …]
|
/dports/science/rdkit/rdkit-Release_2021_03_5/Code/Geometry/ |
H A D | testGrid.cpp | 243 unsigned int nxi, nyi, nzi; in testUniformGridIndexing() local 244 grd.getGridIndices(idx, nxi, nyi, nzi); in testUniformGridIndexing() 245 TEST_ASSERT(nxi == xi); in testUniformGridIndexing() 252 unsigned int nxi, nyi, nzi; in testUniformGridIndexing() local 254 TEST_ASSERT(nxi == xi); in testUniformGridIndexing() 261 unsigned int nxi, nyi, nzi; in testUniformGridIndexing() local 263 TEST_ASSERT(nxi == xi); in testUniformGridIndexing() 270 unsigned int nxi, nyi, nzi; in testUniformGridIndexing() local 272 TEST_ASSERT(nxi == xi); in testUniformGridIndexing() 279 unsigned int nxi, nyi, nzi; in testUniformGridIndexing() local [all …]
|