/dports/math/hpipm/hpipm-0.1.1/experimental/andrea/casadi_vde_interface/ |
H A D | main.c | 6 const int NX = 4; in main() local 10 double in[NX + NU + NX*NX + NX*NU]; in main() 11 double out[NX + NX*NX + NX*NU]; in main() 13 for (int i = 0; i < NX + NU + NX*NX + NX*NU; i++) in[i] = i; in main() 16 for (int i = 0; i < NX; i++) printf("x[%i] = %f\n", i, in[i]); in main() 18 for (int i = 0; i < NX*NX; i++) printf("dx/dx[%i] = %f\n", i, in[i+NX]); in main() 20 for (int i = 0; i < NX*NU; i++) printf("dx/du[%i] = %f\n", i, in[i+NX*NX]); in main() 22 for (int i = 0; i < NU; i++) printf("u[%i] = %f\n", i, in[i+NX*NX+NX*NU]); in main() 28 for (int i = 0; i < NX; i++) printf("x[%i] = %f\n", i, out[i]); in main() 30 for (int i = 0; i < NX*NX; i++) printf("dx/dx[%i] = %f\n", i, out[i+NX]); in main() [all …]
|
/dports/math/adept/adept-2.0.8/benchmark/ |
H A D | advection_schemes_K.h | 22 real q[NX], real jacobian[NX*NX]) { in lax_wendroff_K() argument 24 real flux_K[NX-1][NX]; // Flux Jacobian (dflux/dq_init) in lax_wendroff_K() 26 real q_K[NX][NX]; in lax_wendroff_K() 30 for (int i=0; i<NX; i++) { in lax_wendroff_K() 50 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in lax_wendroff_K() 69 void toon_K(int nt, real c, const real q_init[NX], real q[NX], in toon_K() argument 70 real jacobian[NX*NX]) { in toon_K() argument 72 real flux_K[NX-1][NX]; in toon_K() 73 real q_K[NX][NX]; in toon_K() 75 for (int i=0; i<NX; i++) { in toon_K() [all …]
|
H A D | advection_schemes_AD.h | 19 void lax_wendroff_AD(int nt, real c, const real q_init[NX], real q[NX], in lax_wendroff_AD() argument 20 const real q_AD_const[NX], real q_init_AD[NX]) { in lax_wendroff_AD() argument 22 real flux[NX-1]; in lax_wendroff_AD() 30 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in lax_wendroff_AD() 33 real q_AD[NX]; in lax_wendroff_AD() 67 const real q_AD_const[NX], real q_init_AD[NX]) { in toon_AD() argument 69 real flux[NX-1]; in toon_AD() 81 q += NX; in toon_AD() 83 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in toon_AD() 88 real q_AD[NX]; in toon_AD() [all …]
|
H A D | differentiator.h | 213 jac.resize(NX*NX); 296 jac.resize(NX*NX); 418 jac_.resize(NX*NX); 438 jac = myalloc2(NX,NX); 450 MINDEC(rc,fov_reverse(1, NX, NX, NX, I, jac)); 453 int rc = fov_forward(1, NX, NX, NX, &x[0], I, result, jac); 459 ::jacobian(1, NX, NX, &x[0], jac); 545 jac.resize(NX*NX); 546 jac_transpose_.resize(NX*NX); 579 Real (&jac_transpose2)[NX][NX] [all …]
|
H A D | advection_schemes.h | 25 void lax_wendroff(int nt, Real c, const aReal q_init[NX], aReal q[NX]) { in lax_wendroff() argument 31 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in lax_wendroff() 37 void toon(int nt, Real c, const aReal q_init[NX], aReal q[NX]) { in toon() argument 69 aReal q[NX]) { 73 my_vector Q(NX); 74 my_vector F(NX-1); 84 Q(0) = Q(NX-2); 85 Q(NX-1) = Q(1); 94 my_vector Q(NX); 106 Q(0) = Q(NX-2); [all …]
|
/dports/science/siesta/siesta-4.1.5/Util/Denchar/Src/ |
H A D | ranger.f | 220 . X0(NX), XA(NX,NA), XIJ(NX,*) 424 allocate(CELMSH(NX*NX)) 425 call memory('A','D',NX*NX,'ranger') 428 allocate(RCELL(NX*NX)) 429 call memory('A','D',NX*NX,'ranger') 432 allocate(RMCELL(NX*NX)) 433 call memory('A','D',NX*NX,'ranger') 435 allocate(INX(NX)) 445 allocate(DMX(NX)) 447 allocate(DX(NX)) [all …]
|
H A D | neighb.f | 95 INTEGER NX 96 PARAMETER ( NX = 3 ) 100 real(dp) CELL(NX,NX), RANGE, R2IJ(*), 101 . XA(NX,NA), XIJ(NX,*) 107 real(dp) CELAST(NX,NX), RGLAST, X0(NX) 112 DATA X0 / NX*0.D0 / 119 DO 20 IX = 1,NX 120 DO 10 JX = 1,NX 130 DO 40 IX = 1,NX 131 DO 30 JX = 1,NX [all …]
|
/dports/science/siesta/siesta-4.1.5/Util/STM/ol-stm/Src/ |
H A D | ranger.f | 220 . X0(NX), XA(NX,NA), XIJ(NX,*) 424 allocate(CELMSH(NX*NX)) 425 call memory('A','D',NX*NX,'ranger') 428 allocate(RCELL(NX*NX)) 429 call memory('A','D',NX*NX,'ranger') 432 allocate(RMCELL(NX*NX)) 433 call memory('A','D',NX*NX,'ranger') 435 allocate(INX(NX)) 445 allocate(DMX(NX)) 447 allocate(DX(NX)) [all …]
|
H A D | neighb.f | 94 INTEGER NX 95 PARAMETER ( NX = 3 ) 99 DOUBLE PRECISION CELL(NX,NX), RANGE, R2IJ(*), 100 . XA(NX,NA), XIJ(NX,*) 106 DOUBLE PRECISION CELAST(NX,NX), RGLAST, X0(NX) 111 DATA X0 / NX*0.D0 / 118 DO 20 IX = 1,NX 119 DO 10 JX = 1,NX 129 DO 40 IX = 1,NX 130 DO 30 JX = 1,NX [all …]
|
/dports/misc/libemos/libemos-4.5.9-Source/interpolation/ |
H A D | intocn.F | 297 NX = NX*2 315 RLVAL = PSEC4(NX/2+(J-1)*NX/2) 331 RLVAL = PSEC4(NX/2+(J-1)*NX/2) 335 TARRAY(NX-1+(J-1)*NX) = PSEC4(NX/2+(J-1)*NX/2) 374 NX = NX*2 399 RLVAL = PSEC4(NX/2+(NY-1)*NX/2) 403 TARRAY(NX-1+(NY-1)*NX) = PSEC4(NX/2+(NY-1)*NX/2) 415 RUVAL = PSEC4(NX/2+(J-2)*NX/2) 416 RDVAL = PSEC4(NX/2+(J )*NX/2) 431 RLVAL = PSEC4(NX/2+(1-1)*NX/2) [all …]
|
H A D | intocnu.F | 312 NX = NX*2 330 RLVAL = INOCEAN(NX/2+(J-1)*NX/2) 346 RLVAL = INOCEAN(NX/2+(J-1)*NX/2) 350 TARRAY(NX-1+(J-1)*NX) = INOCEAN(NX/2+(J-1)*NX/2) 389 NX = NX*2 414 RLVAL = INOCEAN(NX/2+(NY-1)*NX/2) 418 TARRAY(NX-1+(NY-1)*NX) = INOCEAN(NX/2+(NY-1)*NX/2) 430 RUVAL = INOCEAN(NX/2+(J-2)*NX/2) 431 RDVAL = INOCEAN(NX/2+(J )*NX/2) 446 RLVAL = INOCEAN(NX/2+(1-1)*NX/2) [all …]
|
/dports/math/adept/adept-2.0.8/test/ |
H A D | test_thread_safe.cpp | 32 #define NX 128 macro 38 toon(int nt, double c, const adouble q_init[NX], adouble q[NX]) { in toon() argument 45 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in toon() 62 adouble q[NX]; // Final values in compute() 65 for (int j = 0; j < NX; j++) { in compute() 80 Real jac_for[NX*NX]; // Where Jacobian will be stored from forward computation in compute() 81 Real jac_rev[NX*NX]; // Where Jacobian will be stored from reverse computation in compute() 87 for (int j = 0; j < NX*NX; j++) { in compute() 93 rmsd = sqrt(rmsd / (NX*NX)); in compute() 148 double q_init_save[NX]; in main() [all …]
|
H A D | test_checkpoint.cpp | 23 #define NX 100 macro 29 toon(int nt, double c, const adouble q_init[NX], adouble q[NX]) { in toon() argument 36 q[0] = q[NX-2]; q[NX-1] = q[1]; // Treat boundary conditions in toon() 59 double q_init_save[NX]; in main() 63 for (int i = 0; i < NX; i++) { in main() 90 double dJ_dq[NX]; in main() 93 for (int i = 0; i < NX; i++) { in main() 104 for (int i = 0; i < NX; i++) { in main() 148 adouble q_save[nblocks][NX]; in main() 151 adouble q[NX]; in main() [all …]
|
/dports/science/pnetcdf/parallel-netcdf-1.8.1/test/testcases/ |
H A D | ncmpi_vars_null_stride.c | 35 #define NX 2 macro 95 count[1] = NX; in main() 110 count[1] = NX; in main() 124 count[1] = NX; in main() 135 count[1] = NX; in main() 164 if (buf[i*nprocs*NX+j*NX+k] != j+10) { in main() 166 __LINE__,i*nprocs*NX+j*NX+k, j+10, buf[i*nprocs*NX+j*NX+k]); in main() 187 if (buf[i*nprocs*NX+j*NX+k] != j+10) { in main() 189 __LINE__,i*nprocs*NX+j*NX+k, j+10, buf[i*nprocs*NX+j*NX+k]); in main() 233 if (buf[i*nprocs*NX+j*NX+k] != j+10) { in main() [all …]
|
/dports/math/cloog/cloog-cloog-0.20.0/test/bigs/ |
H A D | applu.N.w2p21 | 9 # NX NY NZ 1 13 NX NY NZ # parameter names 21 # K J I M NX NY NZ 1 35 # K J I NX NY NZ 1 47 # K J I NX NY NZ 1 59 # K J I NX NY NZ 1 71 # K J I NX NY NZ 1 83 # K J I NX NY NZ 1 95 # K J I NX NY NZ 1 107 # K J I NX NY NZ 1 [all …]
|
H A D | applu.N.w2p12 | 9 # NZ NY NX 1 13 NZ NY NX # parameter names 35 # K NZ NY NX 1 43 # K J NZ NY NX 1 53 # K J I NZ NY NX 1 79 # K J I NZ NY NX 1 467 # K NZ NY NX 1 475 # K I NZ NY NX 1 899 # J NZ NY NX 1 907 # J I NZ NY NX 1 [all …]
|
/dports/science/siesta/siesta-4.1.5/Src/ |
H A D | mneighb.f | 30 integer :: INX(NX), I1NX(NX), I2NX(NX), 31 & J1NX(NX), J2NX(NX), I1EMX(NX), 32 & I2EMX(NX), IMX(NX), I1MX(NX), 33 & I2MX(NX), NEMX(NX), NMX(NX), NNX(NX) 34 real(dp) :: DMX(NX), DX(NX), DX0M(NX), 35 & CELMSH(NX*NX), RCELL(NX*NX), 36 & RMCELL(NX*NX) 147 do IX = 1,NX 158 do IX = 1,NX 403 & CELL(*), range, X0(NX), XA(NX,NA) [all …]
|
/dports/science/silo/silo-4.10.2-bsd/tests/ |
H A D | group_test.c | 77 #define NX 30 macro 742 d[k * NX * NY + j * NX + i] = dist; in build_block_rect3d() 864 d[n_z * NX * NY + n_y * NX + n_x]; in build_block_rect3d() 866 p[n_z * NX * NY + n_y * NX + n_x]; in build_block_rect3d() 1017 #undef NX in build_block_ucd3d() 1018 #define NX 30 in build_block_ucd3d() macro 1569 #undef NX in build_block_curv3d() 1570 #define NX 30 in build_block_curv3d() macro 1780 matlist[k * NX * NY + j * NX + i] = 1; in build_block_curv3d() 1784 matlist[k * NX * NY + j * NX + i] = 2; in build_block_curv3d() [all …]
|
H A D | multispec.c | 79 #define NX 30 macro 143 int speclist[NX*NY],speclist2[NX*NY]; in build_dbfile() 234 dims[0] = NX + 1; in build_dbfile() 258 zdims[0] = NX; in build_dbfile() 260 for (i = 0; i < NX; i++) in build_dbfile() 264 xave = (x[(j) * (NX + 1) + i] + x[(j) * (NX + 1) + i + 1] + in build_dbfile() 265 x[(j + 1) * (NX + 1) + i + 1] + x[(j + 1) * (NX + 1) + i]) / 4.; in build_dbfile() 266 yave = (y[(j) * (NX + 1) + i] + y[(j) * (NX + 1) + i + 1] + in build_dbfile() 267 y[(j + 1) * (NX + 1) + i + 1] + y[(j + 1) * (NX + 1) + i]) / 4.; in build_dbfile() 270 d[j * NX + i] = dist*((float)i/(float)NX); in build_dbfile() [all …]
|
H A D | multi_test.c | 60 #define NX 30 macro 2459 d[k * NX * NY + j * NX + i] = dist; in build_block_rect3d() 2580 d[n_z * NX * NY + n_y * NX + n_x]; in build_block_rect3d() 2582 p[n_z * NX * NY + n_y * NX + n_x]; in build_block_rect3d() 2765 float xstrip[NX + NY + NZ], ystrip[NX + NY + NZ], zstrip[NX + NY + NZ]; in build_block_ucd3d() 2865 d[i * NX * NY + j * NX + k] = dist; in build_block_ucd3d() 3403 d[k * (NX) * (NY) + j * (NX) + i] = dist; in build_block_curv3d() 3456 matlist[k * NX * NY + j * NX + i] = 1; in build_block_curv3d() 3460 matlist[k * NX * NY + j * NX + i] = 2; in build_block_curv3d() 3464 matlist[k * NX * NY + j * NX + i] = 3; in build_block_curv3d() [all …]
|
H A D | simple.c | 66 #undef NX 67 #define NX 10 macro 107 float val[NX * NY * NZ]; in main() 108 long long lval[NX * NY * NZ]; in main() 110 float val2[NX * NY * NZ]; in main() 145 dims[2] = NX; in main() 150 val[i + j * NX + k * NX * NY] = i + j * NX + k * NX * NY; in main() 151 lval[i + j * NX + k * NX * NY] = ((long long) 1 << 35) + i + j * NX + k * NX * NY; in main() 225 if (val2[i + j * NX2 + k * NX2 * NY2] != val[i + j * NX + k * NX * NY]) in main() 235 k2 = NX * NY * NZ; in main() [all …]
|
/dports/science/dcl/dcl-5.4.8/demo/grph2/g2pack/ |
H A D | g2pk02.f | 8 REAL UX(NX), UY(NY) 9 REAL UYW(NX), UXW(NY) 10 REAL CX(NX,NY), CY(NX,NY) 11 REAL Z(NX,NY) 12 REAL TERRAIN(NX) 19 DO 10 I=1,NX 29 DO 20 I=1,NX 51 DO 30 I=1,NX 59 CALL UETONE(Z, NX, NX, NY) 64 DO 40 I=1,NX [all …]
|
H A D | g2pk01.f | 6 PARAMETER(NX=15,NY=15) 8 REAL UX(NX), UY(NY) 10 REAL CX(NX,NY), CY(NX,NY) 11 REAL Z(NX,NY) 16 DO 10 I=1,NX 24 DO 20 I=1,NX 55 DO 30 I=1,NX 61 CALL UWSGXA(UX,NX) 63 CALL UETONE(Z, NX, NX, NY) 68 DO 40 I=1,NX [all …]
|
/dports/science/dcl/dcl-5.4.8/demo/grph1/scpack/ |
H A D | scpkt7.f | 6 PARAMETER( NX=37 ) 10 REAL ALON(NX), ALAT(NX), R(NX), A(NX) 13 DO 10 I=1,NX 14 ALON(I) = XMIN + (XMAX-XMIN)*(I-1)/(NX-1) 15 ALAT(I) = YMIN + (YMAX-YMIN)*(I-1)/(NX-1) 37 DO 100 I=1, NX 38 DO 110 J=1, NX 41 CALL SCPLU(NX, R, ALON, A) 44 DO 200 I=1, NX 45 DO 210 J=1, NX [all …]
|
/dports/graphics/pgplot/pgplot/examples/ |
H A D | pgdemo7.f | 47 NX = KX 53 DO 10 I=1,NX 58 DELTAX = SIZE/(NX+NY) 66 MX = NX+1 96 DO 80 I=2,NX 133 DO 90 I=2,NX 149 DO 100 I=2,NX 219 170 KI = NX 247 10 IF(LEFT.LE.NX .OR. MOD(LEFT-1,NX).EQ.0 .OR. 248 1 RIGHT.LE.NX .OR. MOD(RIGHT-1,NX).EQ.0) GOTO 100 [all …]
|