/dports/math/cppad/CppAD-20210000.8/cppad_ipopt/example/ |
H A D | ode_fast.hpp | 59 ADVector y(Ny), a(Na); in eval_r() 68 a[j] = u[Ny + j]; in eval_r() 74 { ADVector r(Ny), F(Ny); in eval_r() 86 ADVector ym(Ny), G(Ny), Gm(Ny), r(Ny); in eval_r() 93 for(j = 0; j < Ny; j++) in eval_r() 95 y[j] = u[Ny + j]; in eval_r() 101 for(j = 0; j < Ny; j++) in eval_r() 121 return Ny; in range_size() 169 for(j = 0; j < Ny; j++) in index() 170 I[j] = 1 + Ny + i * Ny + j; in index() [all …]
|
H A D | ode_simple.hpp | 62 size_t ny_inx = (S_[Nz] + 1) * Ny; in eval_r() 74 ADVector y(Ny); in eval_r() 78 { for(j = 0; j < Ny; j++) in eval_r() 79 y[j] = x[Ny*S_[k+1] + j]; in eval_r() 84 for(j = 0; j < Ny; j++) in eval_r() 89 ADVector ym(Ny), G(Ny), Gm(Ny); in eval_r() 98 { size_t Index = (j + S_[k]) * Ny; in eval_r() 104 for(i = 0; i < Ny; i++) in eval_r() 105 y[i] = x[Ny + Index + i]; in eval_r() 109 for(i = 0; i < Ny; i++) in eval_r() [all …]
|
H A D | ode2.omh | 281 y^i & = & ( x_{Ny * i} , \ldots , x_{Ny * i + Ny - 1} ) 375 x_{Ny * S(Nz) + Ny } \; , \; 377 x_{Ny * S(Nz) + Ny + Na - 1} 386 Ny * S(k+1) + Ny - 1 \; , \; 387 Ny * S(Nz) + Ny \; , \; 389 Ny * S(Nz) + Ny + Na - 1 428 x_{Ny * S(Nz) + Ny } \; , \; 440 Ny * S(Nz) + Ny \; , \; 442 Ny * S(Nz) + Ny + Na - 1 510 Ny * S(Nz) + Ny \; , \; [all …]
|
H A D | ode_check.cpp | 42 NumberVector a(Na), y_0(Ny), y_1(Ny), y_2(Ny); in ode_check() 45 for(j = 0; j < Ny; j++) in ode_check() 47 y_1[j] = x[Ny + j]; in ode_check() 48 y_2[j] = x[2 * Ny + j]; in ode_check() 68 for(j = 0; j < Ny; j++) in ode_check() 76 for(j = 0; j < Ny; j++) in ode_check() 85 for(j = 0; j < Ny; j++) in ode_check() 92 NumberVector y_i(Ny); in ode_check() 94 { for(j = 0; j < Ny; j++) in ode_check() 95 y_i[j] = x[S[k+1] * Ny + j]; in ode_check()
|
/dports/science/lammps/lammps-stable_29Sep2021/tools/phonon/ |
H A D | interpolate.cpp | 11 Ny = ny; in Interpolate() 13 Npt = Nx*Ny*Nz; in Interpolate() 51 int im = (ii-1+Nx)%Nx, jm = (jj-1+Ny)%Ny, km = (kk-1+Nz)%Nz; in tricubic_init() 141 vidx[0] = (ix*Ny+iy)*Nz+iz; in tricubic() 142 vidx[1] = (ixp*Ny+iy)*Nz+iz; in tricubic() 143 vidx[2] = (ix*Ny+iyp)*Nz+iz; in tricubic() 200 q[1] *= double(Ny); in trilinear() 204 iy = int(q[1])%Ny; in trilinear() 207 iyp = (iy+1)%Ny; in trilinear() 297 int ip1 = p1*Ny*Nz, ip2 = p2*Ny*Nz; in reset_gamma() [all …]
|
/dports/devel/ispc/ispc-1.16.1/examples/cpu/stencil/ |
H A D | stencil.cpp | 56 for (int y = 0; y < Ny; ++y) in InitData() 66 int Nx = 256, Ny = 256, Nz = 256; in main() local 73 Ny *= scale; in main() 84 Aserial[0] = new float[Nx * Ny * Nz]; in main() 85 Aserial[1] = new float[Nx * Ny * Nz]; in main() 86 Aispc[0] = new float[Nx * Ny * Nz]; in main() 87 Aispc[1] = new float[Nx * Ny * Nz]; in main() 92 InitData(Nx, Ny, Nz, Aispc, vsq); in main() 109 InitData(Nx, Ny, Nz, Aispc, vsq); in main() 127 InitData(Nx, Ny, Nz, Aserial, vsq); in main() [all …]
|
/dports/devel/upp/upp/bazaar/plugin/gdal/frmts/grib/degrib18/degrib/ |
H A D | scan.c | 58 sInt4 Ny) in ScanIndex2XY() argument 64 x = row / Ny; in ScanIndex2XY() 66 y = (Ny - 1) - (row % Ny); in ScanIndex2XY() 68 y = row % Ny; in ScanIndex2XY() 82 y = (Ny - 1 - y); in ScanIndex2XY() 126 sInt4 Ny) in XY2ScanIndex() argument 137 y = (Ny - 1 - y); in XY2ScanIndex() 141 row = Ny - 1 - y + x * Ny; in XY2ScanIndex() 143 row = y + x * Ny; in XY2ScanIndex() 187 sInt4 Nx = 2, Ny = 3; in main() local [all …]
|
/dports/graphics/p5-Geo-GDAL/gdal-3.3.3/frmts/grib/degrib/degrib/ |
H A D | scan.c | 58 sInt4 Ny) in ScanIndex2XY() argument 64 x = row / Ny; in ScanIndex2XY() 66 y = (Ny - 1) - (row % Ny); in ScanIndex2XY() 68 y = row % Ny; in ScanIndex2XY() 82 y = (Ny - 1 - y); in ScanIndex2XY() 127 sInt4 Ny) in XY2ScanIndex() argument 138 y = (Ny - 1 - y); in XY2ScanIndex() 142 row = Ny - 1 - y + x * Ny; in XY2ScanIndex() 144 row = y + x * Ny; in XY2ScanIndex() 189 sInt4 Nx = 2, Ny = 3; in main() local [all …]
|
/dports/graphics/gdal/gdal-3.3.3/frmts/grib/degrib/degrib/ |
H A D | scan.c | 58 sInt4 Ny) in ScanIndex2XY() argument 64 x = row / Ny; in ScanIndex2XY() 66 y = (Ny - 1) - (row % Ny); in ScanIndex2XY() 68 y = row % Ny; in ScanIndex2XY() 82 y = (Ny - 1 - y); in ScanIndex2XY() 127 sInt4 Ny) in XY2ScanIndex() argument 138 y = (Ny - 1 - y); in XY2ScanIndex() 142 row = Ny - 1 - y + x * Ny; in XY2ScanIndex() 144 row = y + x * Ny; in XY2ScanIndex() 189 sInt4 Nx = 2, Ny = 3; in main() local [all …]
|
/dports/graphics/py-gdal/gdal-3.3.3/frmts/grib/degrib/degrib/ |
H A D | scan.c | 58 sInt4 Ny) in ScanIndex2XY() argument 64 x = row / Ny; in ScanIndex2XY() 66 y = (Ny - 1) - (row % Ny); in ScanIndex2XY() 68 y = row % Ny; in ScanIndex2XY() 82 y = (Ny - 1 - y); in ScanIndex2XY() 127 sInt4 Ny) in XY2ScanIndex() argument 138 y = (Ny - 1 - y); in XY2ScanIndex() 142 row = Ny - 1 - y + x * Ny; in XY2ScanIndex() 144 row = y + x * Ny; in XY2ScanIndex() 189 sInt4 Nx = 2, Ny = 3; in main() local [all …]
|
/dports/math/cloog/cloog-cloog-0.20.0/test/ |
H A D | pouchet.c | 2 if (Ny >= 2) { 3 for (c0=1;c0<=floord(Ny+4,2);c0++) { 4 for (c1=max(ceild(c0+1,2),c0-1);c1<=min(floord(2*c0+Ny,4),c0);c1++) { 5 if (c0 >= ceild(4*c1-Ny+1,2)) { 11 if (2*c0 == 4*c1-Ny) { 13 if (Ny%2 == 0) { 14 if ((2*c0+3*Ny)%4 == 0) { 15 S2(((2*c0-Ny)/4),((2*c0+Ny)/4),((2*c0-Ny)/2),(Ny-1),c2);
|
H A D | pouchet.good.c | 13 void test(int Ny) in test() argument 19 if (Ny >= 2) { in test() 20 for (c0=1;c0<=floord(Ny+4,2);c0++) { in test() 21 for (c1=max(ceild(c0+1,2),c0-1);c1<=min(floord(2*c0+Ny,4),c0);c1++) { in test() 22 if (c0 >= ceild(4*c1-Ny+1,2)) { in test() 28 if (2*c0 == 4*c1-Ny) { in test() 30 if (Ny%2 == 0) { in test() 31 if ((2*c0+3*Ny)%4 == 0) { in test() 32 S2((2*c0-Ny)/4,(2*c0+Ny)/4,(2*c0-Ny)/2,Ny-1,c2); in test()
|
/dports/misc/adios2/ADIOS2-2.7.1/testing/adios2/engine/bp/ |
H A D | TestBPWriteReadMultiblock.cpp | 848 {Ny / 2, static_cast<size_t>(mpiRank * Nx)}, {Ny - Ny / 2, Nx}); in TEST_F() 1240 {Ny / 2, static_cast<size_t>(mpiRank * Nx)}, {Ny - Ny / 2, Nx}); in TEST_F() 1741 EXPECT_EQ(i8Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1746 EXPECT_EQ(i16Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1751 EXPECT_EQ(i32Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1756 EXPECT_EQ(i64Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1761 EXPECT_EQ(u8Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1766 EXPECT_EQ(u16Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1771 EXPECT_EQ(u32Info[i].Count[0], Ny - Ny / 2); in TEST_F() 1776 EXPECT_EQ(u64Info[i].Count[0], Ny - Ny / 2); in TEST_F() [all …]
|
/dports/science/qmcpack/qmcpack-3.11.0/src/einspline/ |
H A D | multi_bspline_create.c | 156 int Nx, Ny; in create_multi_UBspline_2d_s() local 168 Ny = My+3; in create_multi_UBspline_2d_s() 170 Ny = My+2; in create_multi_UBspline_2d_s() 207 int Nx, Ny; in set_multi_UBspline_2d_s() local 214 Ny = My+3; in set_multi_UBspline_2d_s() 216 Ny = My+2; in set_multi_UBspline_2d_s() 268 Ny = My+3; in create_multi_UBspline_3d_s() 270 Ny = My+2; in create_multi_UBspline_3d_s() 325 Ny = My+3; in set_multi_UBspline_3d_s() 544 int Nx, Ny; in create_multi_UBspline_2d_c() local [all …]
|
H A D | nubspline_create.c | 256 int Mx, My, Nx, Ny; in create_NUBspline_2d_s() local 265 spline->x_stride = Ny; in create_NUBspline_2d_s() 282 int doffset = ix*Ny; in create_NUBspline_2d_s() 283 int coffset = ix*Ny; in create_NUBspline_2d_s() 583 int Mx, My, Nx, Ny; in create_NUBspline_2d_d() local 592 spline->x_stride = Ny; in create_NUBspline_2d_d() 609 int doffset = ix*Ny; in create_NUBspline_2d_d() 610 int coffset = ix*Ny; in create_NUBspline_2d_d() 751 int Mx, My, Nx, Ny; in create_NUBspline_2d_c() local 760 spline->x_stride = Ny; in create_NUBspline_2d_c() [all …]
|
H A D | bspline_create.c | 369 int Nx, Ny; in create_UBspline_2d_s() local 413 int Nx, Ny; in recompute_UBspline_2d_s() local 647 int Nx, Ny; in create_UBspline_2d_c() local 711 int Nx, Ny; in recompute_UBspline_2d_c() local 718 Ny = My+3; in recompute_UBspline_2d_c() 720 Ny = My+2; in recompute_UBspline_2d_c() 1234 int Nx, Ny; in create_UBspline_2d_d() local 1280 int Nx, Ny; in recompute_UBspline_2d_d() local 1288 Ny = My+3; in recompute_UBspline_2d_d() 1545 int Nx, Ny; in create_UBspline_2d_z() local [all …]
|
H A D | multi_nubspline_create.c | 149 int Mx, My, Nx, Ny; in create_multi_NUBspline_2d_s() local 164 spline->x_stride = Ny*N; in create_multi_NUBspline_2d_s() 170 sizeof(float)*Nx*Ny*N); in create_multi_NUBspline_2d_s() 180 int Mx, My, Nx, Ny; in set_multi_NUBspline_2d_s() local 196 coefs+coffset, Ny*ys); in set_multi_NUBspline_2d_s() 409 int Mx, My, Nx, Ny; in create_multi_NUBspline_2d_c() local 444 int Mx, My, Nx, Ny; in set_multi_NUBspline_2d_c() local 705 int Mx, My, Nx, Ny; in create_multi_NUBspline_2d_d() local 736 int Mx, My, Nx, Ny; in set_multi_NUBspline_2d_d() local 989 int Mx, My, Nx, Ny; in create_multi_NUBspline_2d_z() local [all …]
|
H A D | bspline_create_cuda.cu | 49 int Ny = spline->y_grid.num+3; in create_UBspline_3d_c_cuda() local 54 cuda_spline->stride.x = Ny*N; in create_UBspline_3d_c_cuda() 72 for (int iy=0; iy<Ny; iy++) { in create_UBspline_3d_c_cuda() 86 cuda_spline->stride.x = 2*Ny*N; in create_UBspline_3d_c_cuda() 115 int Ny = spline->y_grid.num+3; in create_UBspline_3d_c_cuda_conv() local 119 cuda_spline->stride.x = Ny*N; in create_UBspline_3d_c_cuda_conv() 137 for (int iy=0; iy<Ny; iy++) in create_UBspline_3d_c_cuda_conv() 190 cuda_spline->stride.x = Ny*N; in create_UBspline_3d_s_cuda() 209 for (int iy=0; iy<Ny; iy++) in create_UBspline_3d_s_cuda() 250 cuda_spline->stride.x = Ny*N; in create_UBspline_3d_s_cuda_conv() [all …]
|
/dports/graphics/krita/krita-4.4.8/libs/image/3rdparty/einspline/ |
H A D | multi_bspline_create.cpp | 168 int Nx, Ny; in create_multi_UBspline_2d_s() local 180 Ny = My+3; in create_multi_UBspline_2d_s() 182 Ny = My+2; in create_multi_UBspline_2d_s() 219 int Nx, Ny; in set_multi_UBspline_2d_s() local 226 Ny = My+3; in set_multi_UBspline_2d_s() 228 Ny = My+2; in set_multi_UBspline_2d_s() 280 Ny = My+3; in create_multi_UBspline_3d_s() 282 Ny = My+2; in create_multi_UBspline_3d_s() 335 Ny = My+3; in set_multi_UBspline_3d_s() 473 int Nx, Ny; in create_multi_UBspline_2d_c() local [all …]
|
H A D | bspline_create.cpp | 389 int Nx, Ny; in create_UBspline_2d_s() local 433 int Nx, Ny; in recompute_UBspline_2d_s() local 666 int Nx, Ny; in create_UBspline_2d_c() local 730 int Nx, Ny; in recompute_UBspline_2d_c() local 737 Ny = My+3; in recompute_UBspline_2d_c() 739 Ny = My+2; in recompute_UBspline_2d_c() 1266 int Nx, Ny; in create_UBspline_2d_d() local 1312 int Nx, Ny; in recompute_UBspline_2d_d() local 1320 Ny = My+3; in recompute_UBspline_2d_d() 1567 int Nx, Ny; in create_UBspline_2d_z() local [all …]
|
H A D | nubspline_create.cpp | 287 int My, Nx, Ny; in create_NUBspline_2d_s() local 296 spline->x_stride = Ny; in create_NUBspline_2d_s() 313 int doffset = ix*Ny; in create_NUBspline_2d_s() 314 int coffset = ix*Ny; in create_NUBspline_2d_s() 340 int My, Mz, Nx, Ny, Nz; in create_NUBspline_3d_s() local 631 int My, Nx, Ny; in create_NUBspline_2d_d() local 640 spline->x_stride = Ny; in create_NUBspline_2d_d() 657 int doffset = ix*Ny; in create_NUBspline_2d_d() 658 int coffset = ix*Ny; in create_NUBspline_2d_d() 802 int My, Nx, Ny; in create_NUBspline_2d_c() local [all …]
|
/dports/science/py-pyscf/pyscf-2.0.1/examples/pbc/ |
H A D | 40-customize_kpts_hamiltonian.py | 22 Ny=2 variable 27 Nsite=Nx*Ny 36 [0.0 ,Ny*1.0,0.0], 50 def gen_H_tb(t,Nx,Ny,kvec): argument 51 H = np.zeros((Nx,Ny,Nx,Ny),dtype=np.complex) 53 for j in range(Ny): 64 if j == Ny-1: 70 H[i,j,i,Ny-1] += np.exp(-1j*np.dot(np.array(kvec),np.array([0,-Ny]))) 73 return -t*H.reshape(Nx*Ny,Nx*Ny) 77 def get_H_tb_array(kpts,Nx,Ny,t): argument [all …]
|
/dports/science/dakota/dakota-6.13.0-release-public.src-UI/packages/external/trilinos/packages/rol/example/PDE-OPT/ginzburg-landau/ |
H A D | plotresults.m | 2 function plotresults(Nx,Ny,x0,x1,y0,y1) 4 Ny = Ny+1; variable 44 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(sqrt(stateR.^2 + stateI.^2),N… 50 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(stateR,Nx,Ny)); 56 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(stateI,Nx,Ny)); 129 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(sqrt(stateR.^2 + stateI.^2),N… 135 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(stateR,Nx,Ny)); 141 contour(reshape(nodes(:,1),Nx,Ny), reshape(nodes(:,2),Nx,Ny), reshape(stateI,Nx,Ny));
|
/dports/science/py-GPy/GPy-1.10.0/GPy/likelihoods/ |
H A D | binomial.py | 69 Ny = N-y 73 t2[Ny>0] = Ny[Ny>0]*np.log(1.-inv_link_f[Ny>0]) 95 Ny = N-y 99 t2[Ny>0] = (Ny[Ny>0])/(1.-inv_link_f[Ny>0]) 126 Ny = N-y 130 t2[Ny>0] = -(Ny[Ny>0])/np.square(1.-inv_link_f[Ny>0]) 158 Ny = N-y 162 t2[Ny>0] = - 2*(Ny[Ny>0])/(1.-inv_link_f[Ny>0])**3
|
/dports/devel/llvm-cheri/llvm-project-37c49ff00e3eadce5d8703fdc4497f28458c64a8/polly/lib/External/isl/test_inputs/codegen/cloog/ |
H A D | pouchet.c | 1 for (int c0 = 1; c0 <= floord(Ny, 2) + 2; c0 += 1) 2 for (int c1 = max(c0 - 1, c0 / 2 + 1); c1 <= min(c0, (Ny + 2 * c0) / 4); c1 += 1) { 3 if (Ny + 2 * c0 >= 4 * c1 + 1) { 10 S2((-Ny + 2 * c0) / 4, (Ny + 2 * c0) / 4, (-Ny / 2) + c0, Ny - 1, c2);
|