2 [1] 7 5 C 2.1 (kikuchi-fem-mac ) input.dat (cat input.dat type input.dat ))
|
|
|
- あつの あかさか
- 6 years ago
- Views:
Transcription
1 , , (1) (2) (3) 2 Ω Poisson u = f (in Ω), u = 0 (on Γ 1 ), u n = 0 (on Γ 2) f Γ 1, Γ 2 Γ := Ω n Γ Ω (1980) Fortran : kikuchi-fem-mac.tar.gz 1 ( fem, 6701) tar xzf kikuchi-fem-mac.tar.gz kikuchi-fem-v2 cd kikuchi-fem-mac ls naive.c, band.c, input.dat, make-input.c Windows Shift JIS Linux Mac (nkf make euc, make utf8 ) 1 1
2 2 [1] 7 5 C 2.1 (kikuchi-fem-mac ) input.dat (cat input.dat type input.dat )) : 2 (input.dat ) 1. 1 (nnode) nelmt Dirichlet Γ 1 (nbc) 2. ( ) (x i, y i ) (i = 0, 1,, nnode 1) 3. ( ) (0 nelmt 1 ) ( ) 4. Γ 1 2
3 make-input.c /* * make-input.c */ #include <stdio.h> #include <math.h> int main(void) { int i, j, n; int elmt1, elmt11, elmt12, elmt13, elmt2, elmt21, elmt22, elmt23; double h; /* n */ scanf("%d", &n); h = 1.0 / n; /*,, */ printf("%d %d %d\n", (n + 1) * (n + 1), 2 * n * n, 2 * n + 1); /* */ for (i = 0; i <= n; i++) for (j = 0; j <= n; j++) printf("%f %f\n", i * h, j * h); /* */ for (j = 0; j < n; j++) { for (i = 0; i < n; i++) { /* I */ elmt1 = 2 * (i + n * j); elmt11 = i + (n + 1) * j; elmt12 = elmt11 + (n + 1); elmt13 = elmt12 + 1; /* I */ elmt2 = elmt1 + 1; elmt21 = elmt11; elmt22 = elmt21 + (n + 2); elmt23 = elmt21 + 1; printf("%d %d %d ", elmt11, elmt12, elmt13); printf("%d %d %d\n", elmt21, elmt22, elmt23); /* */ for (j = 0; j <= n; j++) printf("%d ", j); for (i = 1; i <= n; i++) printf("%d ", (n + 1) * i); printf("\n"); return 0; ccmg (glsc -d ) 3
4 input4.dat mathpc00% ccmg make-input.c mathpc00%./make-input ( ) mathpc00%./make-input > input4.dat 4 mathpc00% cat input4.dat ( ) 2.2 disp-glsc disp-glsc GLSC ccmg GLSC mathpc00% ccmg disp-glsc.c 4
5 input4.dat mathpc00%./disp-glsc input4.dat input4.dat mathpc00% cat input4.dat./disp-glsc make-input disp-glsc mathpc00%./make-input./disp-glsc 4 2 DISP DISP PostScript mathpc00% g out -i DISP DISP.i00 PostScript PostScript DISP.i00 mathpc00% gv DISP.i00 & mathpc00% gsview32 DISP.i00 & (PostScript ) PostScript DISP.i00 mathpc00% lpr DISP.i00 L A TEX \includegraphics[angle=90,width=10cm]{disp.i00 (90 10 cm ) GLSCWIN glsc -d GLSC (6701 glscd ) 2 GLSC g sleep() 5
6 mathpc00% glsc -d disp-glsc mathpc00%./disp-glsc < input4.dat ( < ) mathpc00% cat input4.dat./disp-glsc mathpc00%./make-input./disp-glsc 4 3 DISP000.emf L A TEX mathpc00% convert DISP000.emf DISP000.eps (ImageMagick ) disp-glsc input4.dat disp-glsc mathpc00% cat input4.dat GLSCWIN Windows 6
7 mathpc00% cat input4.dat./disp-glsc basic data nnode= 25 nelmt= 32 nbc= 9 nband= 0 x,y-coordinates of nodes ( x,y ) i x(i) y(i) i x(i) y(i) nodes of elements ( )
8 nodes with zero Dirichlet data ( ) mathpc00% ( 2 ) 2: input4.dat 2.3 naive.c [1] FORTRAN C naive.c 1 /* naive.c -- Poisson */ 2 3 /* 4 * 5 * 6 * 7 Fortran C 7 * 8 * 9 * 1, 2 10 * f 11 * 12 * - u=f in 13 * 14 * u=0 on 1, u/ n=0 on 2 15 * 16 * u=u(x,y) n 17 */ #include <stdio.h> 20 #include <stdlib.h> 21 #include <math.h> 22 8
9 23 typedef struct { 24 int node[3]; 25 threenodes; typedef int *ivector; 28 typedef double *vector; 29 typedef vector *matrix; double f(double, double); 32 ivector new_ivector(int); 33 vector new_vector(int); 34 matrix new_matrix(int, int); 35 void errorexit(char *); 36 void input0(file *, int *, int *, int *); 37 void input(file *, int, int, int, threenodes *, vector, vector, ivector); 38 void assem(int, int, int, matrix, vector, 39 threenodes *, vector, vector, ivector); 40 void solve(matrix, vector, int); 41 void ecm(int, threenodes *, vector, vector, double[3][3], double[3]); 42 void output(int, vector); /* main program 45 * the finite element method for the Poisson equation 46 */ int main(int argc, char **argv) 49 { 50 /* 51 * nnode 52 * nelmt 53 * nbc (Dirichlet ) 54 * x[], y[] 55 * ielmt[][3] 56 * ibc[] 57 * am[][] 58 * fm[] 59 * 60 * : 0 61 * 0,1,...,nnode-1 62 * 0,1,...,nelmt-1 63 * 64 */ 65 int nnode, nelmt, nbc; 66 vector x, y, fm; 67 matrix am; 68 ivector ibc; 69 threenodes *ielmt; 70 FILE *fp; if (argc == 2) 73 fp = fopen(argv[1], "r"); 74 else 75 fp = stdin; /* nnode, nelmt, nbc */ 78 input0(fp, &nnode, &nelmt, &nbc); /* */ 9
10 81 if ((ielmt = malloc(sizeof(threenodes) * nelmt)) == NULL) 82 errorexit(" threenodes "); 83 if ((ibc = new_ivector(nnode)) == NULL) 84 errorexit(" ibc "); 85 if ((x = new_vector(nnode)) == NULL) 86 errorexit(" x x "); 87 if ((y = new_vector(nnode)) == NULL) 88 errorexit(" y y "); 89 if ((fm = new_vector(nnode)) == NULL) 90 errorexit(" fm "); 91 if ((am = new_matrix(nnode, nnode)) == NULL) 92 errorexit(" am "); /* */ 95 input(fp, nnode, nelmt, nbc, ielmt, x, y, ibc); /* 1 */ 98 assem(nnode, nelmt, nbc, am, fm, ielmt, x, y, ibc); 99 /* 1 */ 100 solve(am, fm, nnode); 101 /* */ 102 output(nnode, fm); 103 return 0; /*,, Dirichlet */ 107 void input0(file *fp, int *nnode, int *nelmt, int *nbc) 108 { 109 char buf[bufsiz]; 110 fgets(buf, sizeof(buf), fp); 111 sscanf(buf, "%d %d %d", nnode, nelmt, nbc); /*, 115, 116 Dirichlet */ 117 void input(file *fp, int nnode, int nelmt, int nbc, 118 threenodes *ielmt, vector x, vector y, ivector ibc) 119 { 120 int i, j; 121 /* input */ 122 for (i = 0; i < nnode; i++) 123 fscanf(fp, "%lf %lf", &(x[i]), &(y[i])); 124 for (i = 0; i < nelmt; i++) 125 for (j = 0; j < 3; j++) 126 fscanf(fp, "%d", &(ielmt[i].node[j])); 127 for (i = 0; i < nbc; i++) 128 fscanf(fp, "%d", &ibc[i]); /* output of input data */ 131 printf("basic data\n nnode=%4d nelmt=%4d nbc=%4d\n", 132 nnode, nelmt, nbc); 133 printf("x,y-coordinates of nodes ( x,y )\n"); 134 printf(" i x(i) y(i) i x(i) y(i)\n"); 135 for (i = 0; i < nnode; i++) { 136 printf("%4d %8.4f %8.4f", i, x[i], y[i]); 137 if (i % 2 == 1) printf("\n");
11 139 printf("\nnodes of elements ( )\n" 140 " i i1 i2 i3\n i i1 i2 i3\n"); 141 for (i = 0; i < nelmt; i++) { 142 printf("%4d %4d %4d %4d ", 143 i, ielmt[i].node[0], ielmt[i].node[1], ielmt[i].node[2]); 144 if (i % 2 == 1) 145 printf("\n"); if (nbc > 0) { 148 printf("\nnodes with zero Dirichlet data ( )\n"); 149 for (i = 0; i < nbc; i++) 150 printf("%5d", ibc[i]); 151 if (i % 8 == 7) printf("\n"); printf("\n"); /* " " */ 157 void assem(int nnode, int nelmt, int nbc, 158 matrix am, vector fm, threenodes *ielmt, 159 vector x, vector y, ivector ibc) 160 { 161 int i, j, ie, ii, jj; 162 /* the direct stiffness method */ 163 double ae[3][3], fe[3]; 164 /* initial clear */ 165 for (i = 0; i < nnode; i++) { 166 fm[i] = 0.0; 167 for (j = 0; j < nnode; j++) 168 am[i][j] = 0.0; /* assemblage of total matrix and vector; */ 171 for (ie = 0; ie < nelmt; ie++) { 172 ecm(ie, ielmt, x, y, ae, fe); 173 for (i = 0; i < 3; i++) { 174 ii = ielmt[ie].node[i]; 175 fm[ii] += fe[i]; 176 for (j = 0; j < 3; j++) { 177 jj = ielmt[ie].node[j]; 178 am[ii][jj] += ae[i][j]; /* the homogeneous Dirichlet condition */ 183 for (i = 0; i < nbc; i++) { 184 ii = ibc[i]; 185 fm[ii] = 0.0; 186 for (j = 0; j < nnode; j++) 187 am[ii][j] = am[j][ii] = 0.0; 188 am[ii][ii] = 1.0; #ifdef DEBUG 191 for (i = 0; i < nnode; i++) { 192 for (j = 0; j < nnode; j++) 193 printf("%5g ", am[i][j]); 194 printf(" %5g\n", fm[i]); #endif 11
12 /*, */ 200 void ecm(int ie, threenodes *ielmt, vector x, vector y, 201 double ae[3][3], double fe[3]) 202 { 203 int i, j, k; 204 double d, s; 205 /* element coefficient matrix and free vector */ 206 double xe[3], ye[3], b[3], c[3]; 207 /* */ 208 for (i = 0; i < 3; i++) { 209 j = ielmt[ie].node[i]; 210 xe[i] = x[j]; 211 ye[i] = y[j]; /* */ 214 d = xe[0] * (ye[1] - ye[2]) + xe[1] * (ye[2] - ye[0]) xe[2] * (ye[0] - ye[1]); 216 s = fabs(d) / 2.0; 217 /* calculation of element coefficient matrix */ 218 for (i = 0; i < 3; i++) { 219 /* j, k */ 220 j = (i == 2)? 0 : i + 1; 221 k = (i == 0)? 2 : i - 1; 222 b[i] = (ye[j] - ye[k]) / d; 223 c[i] = (xe[k] - xe[j]) / d; for (i = 0; i < 3; i++) 226 for (j = 0; j < 3; j++) 227 ae[i][j] = s * (b[j] * b[i] + c[j] * c[i]); 228 /* calculation of element free vector */ 229 for (i = 0; i < 3; i++) 230 fe[i] = s * f(xe[i], ye[i]) / 3.0; /* Gauss 1 */ 234 void solve(matrix a, vector f, int m) 235 { 236 int m1, i, j, k; 237 double aa; 238 /* forward elimination; */ 239 m1 = m - 1; 240 for (i = 0; i < m1; i++) { 241 for (j = i + 1; j < m; j++) { 242 aa = a[j][i] / a[i][i]; 243 f[j] -= aa * f[i]; 244 for (k = i + 1; k < m; k++) 245 a[j][k] -= aa * a[i][k]; /* backward substitution */ 249 f[m-1] /= a[m-1][m-1]; 250 for (i = m - 2; i >= 0; i--) { 251 for (j = i + 1; j < m; j++) 252 f[i] -= a[i][j] * f[j]; 253 f[i] /= a[i][i];
13 void output(int nnode, vector fm) 258 { 259 int i; 260 /* output of approximate nodal values of u */ 261 printf("nodal values of u ( u )\n"); 262 for (i = 0; i < 3; i++) 263 printf(" i u "); 264 for (i = 0; i < nnode; i++) { 265 if (i % 3 == 0) printf("\n"); 266 printf("%4d %11.3e", i, fm[i]); printf("\n"); double f(double x, double y) 272 { 273 return 1.0; vector new_vector(int n) 277 { 278 return malloc(sizeof(double) * n); ivector new_ivector(int n) 282 { 283 return malloc(sizeof(int) * n); void del_vector(vector a) 287 { 288 free(a); void del_ivector(ivector a) 292 { 293 free(a); matrix new_matrix(int m, int n) 297 { 298 int i; 299 matrix a; 300 if ((a = malloc(m * sizeof(vector))) == NULL) { 301 return NULL; for (i = 0; i < m; i++) { 304 if ((a[i] = new_vector(n)) == NULL) { 305 while (--i >= 0) free(a[i]); 306 free(a); 307 return NULL; return a;
14 313 void errorexit(char *msg) 314 { 315 fprintf(stderr, msg); 316 exit(1); naive.c main() input() assem() A, f ( ) ecm() A e, f e solve() output() ( ) f() Poisson u = f f nnode nelmt nbc (Dirichlet ) x[], y[] ielmt[][3] ibc[] am[][] fm[] ecm() ecm() assem() 1 void ecm(int ie, threenodes *ielmt, vector x, vector y, double ae[3][3], double fe[3]) ie e = e ie A e, f e x[], y[] (P j = (x[j], y[j]) ) /* */ for (i = 0; i < 3; i++) { j = ielmt[ie].node[i]; xe[i] = x[j]; ye[i] = y[j]; e = e ie i N i j (xe[i], ye[i]) 14
15 /* */ d = xe[0] * (ye[1] - ye[2]) + xe[1] * (ye[2] - ye[0]) + xe[2] * (ye[0] - ye[1]); s = fabs(d) / 2.0; e = e ie s ( ) e = 1 2 det x 1 x 0 x 2 x 0 = 1 y 1 y 0 y 2 y 0 2 [(x 1 x 0 )(y 2 y 0 ) (x 2 x 0 )(y 1 y 0 )]. /* calculation of element coefficient matrix */ for (i = 0; i < 3; i++) { /* j, k */ j = (i == 2)? 0 : i + 1; k = (i == 0)? 2 : i - 1; b[i] = (ye[j] - ye[k]) / d; c[i] = (xe[k] - xe[j]) / d; for (i = 0; i < 3; i++) for (j = 0; j < 3; j++) ae[i][j] = s * (b[j] * b[i] + c[j] * c[i]); L i (x, y) = a i + b i x + c i y a i := x jy k x k y j 2 e, b i := y j y k, c i := x k x j 2 e 2 e ((i, j, k) (0, 1, 2) ) (i, j, k) (0, 1, 2) (i, j, k) = (0, 1, 2), (1, 2, 0), (2, 0, 1) A e (i, j) L j, L i e ( L j, L i e = L j L i dxdy = e e b j c j ) ( b i c i ) dxdy = (b j b i + c j c i ) e /* calculation of element free vector */ for (i = 0; i < 3; i++) fe[i] = s * f(xe[i], ye[i]) / 3.0; f e f (e) i := (f, L i ) e f f 2 f(n j )L j j=0 15
16 1 f (e) 0 e 12 (2f(N 0) + f(n 1 ) + f(n 2 )), f (e) 1 e 12 (f(n 0) + 2f(N 1 ) + f(n 2 )), f (e) 2 e 12 (f(n 0) + f(n 1 ) + 2f(N 2 )) assem() void assem(int nnode, int nelmt, int nbc, matrix am, vector fm, threenodes *ielmt, vector x, vector y, ivector ibc) am, fm 0 ( = += 0 ) A = N e 1 k=0 A k, f = N e 1 k=0 f k /* assemblage of total matrix and vector; */ for (ie = 0; ie < nelmt; ie++) { ecm(ie, ielmt, x, y, ae, fe); for (i = 0; i < 3; i++) { ii = ielmt[ie].node[i]; fm[ii] += fe[i]; for (j = 0; j < 3; j++) { jj = ielmt[ie].node[j]; am[ii][jj] += ae[i][j]; ie for ie A ie, f ie am, fm am A, fm f A f (A P i Γ 1 i i )) (f P i Γ 1 i i ) 16
17 A (A P i Γ 1 i i i ) f (f a i g 1 (P i ) ) P i Γ 1 P i Γ 1 i A i, i e i = (0,..., 0, 1, 0,..., 0) f i 0 i u i = 0 ( Dirichlet u = 0 on Γ 1 ) /* the homogeneous Dirichlet condition */ for (i = 0; i < nbc; i++) { ii = ibc[i]; fm[ii] = 0.0; for (j = 0; j < nnode; j++) am[ii][j] = am[j][ii] = 0.0; am[ii][ii] = 1.0; 3 band [1] band.c [1] FORTRAN C band.c kikuchi-fem-mac 3.2 band.c m, i j > m = a ij = 0 m nband m am[nnode][nnode] am[nnode][nband] ecm() naive.c assem() solve() 17
18 3.2.1 assem() nband void assem(int nnode, int nelmt, int nbc, int nband, matrix am, vector fm, threenodes *ielmt, vector x, vector y, ivector ibc) A = (a ij ) ( ) a ij (j i) am[i][j-i] /* assemblage of total matrix and vector; */ for (ie = 0; ie < nelmt; ie++) { ecm(ie, ielmt, x, y, ae, fe); for (i = 0; i < 3; i++) { ii = ielmt[ie].node[i]; fm[ii] += fe[i]; for (j = 0; j < 3; j++) { jj = ielmt[ie].node[j]; if ((ij = jj - ii) >= 0) am[ii][ij] += ae[i][j]; /* the homogeneous Dirichlet condition */ for (i = 0; i < nbc; i++) { ii = ibc[i]; fm[ii] = 0.0; for (j = 0; j < nband; j++) { if ((ij = ii - j) >= 0) am[ij][j] = 0.0; /* a_{ii,jj = 0 */ am[ii][j] = 0.0; am[ii][0] = 1.0; 3.3 band A = (a ij ) m i j > m = a ij = 0 18
19 4 a ij 0 = i j m make-band-input.c mathpc00% ccmg make-band-input.c mathpc00%./make-band-input ( ) mathpc00% ccmg band.c mathpc00%./make-band-input./band 4 mathpc00% ccmg contour.c mathpc00%./contour < band.out 3: 4 4: 20 4 A 3 A 4 max i j a ij 0 19
20 4 (1) (2) (3) (1) naive band (2) u = 0 on Γ 1 u = g 1 on Γ 1 ( naive ) P i Γ 1 a i g 1 (P i ) (3) u n = 0 on Γ 2 u n = g 2 on Γ 2 ( naive ) [g 2, v] = g 2 v dσ Γ (2) naive.c P i Γ 1 i g 1 (P i ) input.data (g 1 0 ) assem() 4.2 (3) naive.c Γ 2 4 6, 7, 8, 5, 2 Γ 2 2 4, 6, 7, 7, 8, 8, 5, 5, 2 N n I j, J j (j = 0,..., N n 1) [g 2, v] = g 2 v ds = Γ 2 N n 1 j=0 P Ij P Jj g 2 v ds. 20
21 P Q g 2 (P ), g 2 (Q), v(p ), v(q) P Q φ(t) = (1 t) OP + t OQ (t [0, 1]) g 2 1 g 2 v ds = P Q P Q 1 0 g 2 (φ(t)) = (1 t)g 2 (P ) + tg 2 (Q), v(φ(t)) = (1 t) v(p ) + t v(q), ds = P Q dt. [(1 t)g 2 (P ) + tg 2 (Q)] [(1 t) v(p ) + t v(q)] dt 1 3 = ( v(p ) v(q))p Q ( ) 6 g 2 (P ) = ( v(p ) v(q)) P Q 1 g 2 (Q) 6 3 ( ) 2g 2 (P ) + g 2 (Q). g 2 (P ) + 2g 2 (Q) Γ 2 P I P J I J g 2 g 2 (P I ) g 2 (P J ) assem() A kikuchi-fem-mac.tar.gz README kikuchi-fem-v2.tar.gz (2010/6/8 ) 1. (?) Windows Shift JIS nkf SunOS, CentOS make euc euc qkc qkc -eu *.c OK Windows make sjis Shift JIS nkf ( ShiftJIS ) 21
22 2. (a) SunOS, CentOS make demo (b) Cygwin Makefile CFLAGS = -W -Wall -O -finput-charset=cp932 -fexec-charset=cp932 CCGLSC make demo (c) ( ) (i) make.bat (glsc -d ) (ii) make-demo.bat : Cygwin glscd Makefile #CFLAGS = -W -Wall -O CFLAGS = -W -Wall -O -finput-charset=cp932 -fexec-charset=cp932 #CCGLSC = ccmg #CCGLSC = glsc -d CCGLSC = glscd [1],, ( 1980, 1999). [2],, suurikeisantokuron/members/fem.pdf, [3], 1 1, linear-eq-1.pdf, [4],, fem,
第7章 有限要素法のプログラミング
April 3, 2019 1 / 34 7.1 ( ) 2 Poisson 2 / 34 7.2 femfp.c [1] main( ) input( ) assem( ) ecm( ) f( ) solve( ) gs { solve( ) output( ) 3 / 34 7.3 fopen() #include FILE *fopen(char *fname, char
1 5 13 4 1 41 1 411 1 412 2 413 3 414 3 415 4 42 6 43 LU 7 431 LU 10 432 11 433 LU 11 44 12 441 13 442 13 443 SOR ( ) 14 444 14 445 15 446 16 447 SOR 16 448 16 45 17 4 41 n x 1,, x n a 11 x 1 + a 1n x
QR
1 7 16 13 1 13.1 QR...................................... 2 13.1.1............................................ 2 13.1.2..................................... 3 13.1.3 QR........................................
Ｃ言語によるアルゴリズムとデータ構造
Algorithms and Data Structures in C 4 algorithm List - /* */ #include List - int main(void) { int a, b, c; int max; /* */ Ÿ 3Ÿ 2Ÿ 3 printf(""); printf(""); printf(""); scanf("%d", &a); scanf("%d",
新版明解Ｃ言語 実践編
2 List - "max.h" a, b max List - max "max.h" #define max(a, b) ((a) > (b)? (a) : (b)) max List -2 List -2 max #include "max.h" int x, y; printf("x"); printf("y"); scanf("%d", &x); scanf("%d", &y); printf("max(x,
double float
2015 3 13 1 2 2 3 2.1.......................... 3 2.2............................. 3 3 4 3.1............................... 4 3.2 double float......................... 5 3.3 main.......................
£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (2018) - Â裵²ó ¨¡ À©¸æ¹½Â¤¡§¾ò·ïʬ´ô ¨¡
(2018) 2018 5 17 0 0 if switch if if ( ) if ( 0) if ( ) if ( 0) if ( ) (0) if ( 0) if ( ) (0) ( ) ; if else if ( ) 1 else 2 if else ( 0) 1 if ( ) 1 else 2 if else ( 0) 1 if ( ) 1 else 2 (0) 2 if else
[1] #include<stdio.h> main() { printf("hello, world."); return 0; } (G1) int long int float ± ±
![[1] #include<stdio.h> main() { printf(hello, world.); return 0; } (G1) int long int float ± ±](/thumbs/99/143405735.jpg "[1] #include<stdio.h> main() { printf(hello, world.); return 0; } (G1) int long int float ± ±")
[1] #include printf("hello, world."); (G1) int -32768 32767 long int -2147483648 2147483647 float ±3.4 10 38 ±3.4 10 38 double ±1.7 10 308 ±1.7 10 308 char [2] #include int a, b, c, d,
comment.dvi
( ) (sample1.c) (sample1.c) 2 2 Nearest Neighbor 1 (2D-class1.dat) 2 (2D-class2.dat) (2D-test.dat) 3 Nearest Neighbor Nearest Neighbor ( 1) 2 1: NN 1 (sample1.c) /* -----------------------------------------------------------------
XMPによる並列化実装2
2 3 C Fortran Exercise 1 Exercise 2 Serial init.c init.f90 XMP xmp_init.c xmp_init.f90 Serial laplace.c laplace.f90 XMP xmp_laplace.c xmp_laplace.f90 #include int a[10]; program init integer
: CR (0x0d) LF (0x0a) line separator CR Mac LF UNIX CR+LF MS-DOS WINDOWS Japan Advanced Institute of Science and Technology
I117 8 1 School of Information Science, Japan Advanced Institute of Science and Technology : CR (0x0d) LF (0x0a) line separator CR Mac LF UNIX CR+LF MS-DOS WINDOWS Japan Advanced Institute of Science and
ex14.dvi
1,, 0, b (b b 2 b ) n k n = n j b j, (0 n j b 1), n =(n k n k 1...n 1 n 0 ) b, n j j j +1, 0,...,b 1 (digit). b b, n b 1 ñ, ñ = k (b 1 n j )b j b N, n b n, n = b N n, n =ñ+1 b N, n m n + m (mod b N ),
9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0)
E-mail: takio-kurita@aist.go.jp 1 ( ) CPU ( ) 2 1. a f(a) =(a 1.0) 2 (1) a ( ) 1(a) f(a) a (1) a f(a) a =2(a 1.0) (2) 2 0 a f(a) a =2(a 1.0) = 0 (3) 1 9 8 7 (x-1.0)*(x-1.0) 6 4 2.0*(x-1.0) 6 2 5 4 0 3-2
(2 Linux Mozilla [ ] [ ] [ ] [ ] URL 2 qkc, nkc ~/.cshrc (emacs 2 set path=($path /usr/meiji/pub/linux/bin tcsh b
![(2 Linux Mozilla [ ] [ ] [ ] [ ] URL 2 qkc, nkc ~/.cshrc (emacs 2 set path=($path /usr/meiji/pub/linux/bin tcsh b](/thumbs/92/108952818.jpg "(2 Linux Mozilla [ ] [ ] [ ] [ ] URL 2 qkc, nkc ~/.cshrc (emacs 2 set path=($path /usr/meiji/pub/linux/bin tcsh b")
II 5 (1 2005 5 26 http://www.math.meiji.ac.jp/~mk/syori2-2005/ UNIX (Linux Linux 1 : 2005 http://www.math.meiji.ac.jp/~mk/syori2-2005/jouhousyori2-2005-00/node2. html ( (Linux 1 2 ( ( http://www.meiji.ac.jp/mind/tool/internet-license/
橡Pro PDF
1 void main( ) char c; /* int c; */ int sum=0; while ((c = getchar())!= EOF) if(isdigit(c) ) sum += (c-'0'); printf("%d\n", sum); main()int i,sum=0; for(i=0;i
II 3 yacc (2) 2005 : Yacc 0 ~nakai/ipp2 1 C main main 1 NULL NULL for 2 (a) Yacc 2 (b) 2 3 y
II 3 yacc (2) 2005 : Yacc 0 ~nakai/ipp2 1 C 1 6 9 1 main main 1 NULL NULL 1 15 23 25 48 26 30 32 36 38 43 45 47 50 52 for 2 (a) 2 2 1 Yacc 2 (b) 2 3 yytext tmp2 ("") tmp2->next->word tmp2 yytext tmp2->next->word
新・明解Ｃ言語 ポインタ完全攻略
2 1-1 1-1 /* 1-1 */ 1 int n = 100; int *p = &n; printf(" n %d\n", n); /* n int */ printf("*&n %d\n", *&n); /* *&n int */ printf(" p %p\n", p); /* p int * */ printf("&*p %p\n", &*p); /* &*p int * */ printf("sizeof(n)
lexex.dvi
(2018, c ) http://istksckwanseiacjp/ ishiura/cpl/ 4 41 1 mini-c lexc,, 2 testlexc, lexc mini-c 1 ( ) mini-c ( ) (int, char, if, else, while, return 6 ) ( ) (+, -, *, /, %, &, =, ==,!=, >, >=,
1 4 2 EP) (EP) (EP)
2003 2004 2 27 1 1 4 2 EP) 5 3 6 3.1.............................. 6 3.2.............................. 6 3.3 (EP)............... 7 4 8 4.1 (EP).................... 8 4.1.1.................... 18 5 (EP)
Microsoft Word - C.....u.K...doc
C uwêííôöðöõ Ð C ÔÖÐÖÕ ÐÊÉÌÊ C ÔÖÐÖÕÊ C ÔÖÐÖÕÊ Ç Ê Æ ~ if eíè ~ for ÒÑÒ ÌÆÊÉÉÊ ~ switch ÉeÍÈ ~ while ÒÑÒ ÊÍÍÔÖÐÖÕÊ ~ 1 C ÔÖÐÖÕ ÐÊÉÌÊ uê~ ÏÒÏÑ Ð ÓÏÖ CUI Ô ÑÊ ÏÒÏÑ ÔÖÐÖÕÎ d ÈÍÉÇÊ ÆÒ Ö ÒÐÑÒ ÊÔÎÏÖÎ d ÉÇÍÊ
新・明解Ｃ言語 実践編
第 1 章 見 21 1-1 見えないエラー 見 List 1-1 "max2x1.h" a, b max2 List 1-1 chap01/max2x1.h max2 "max2x1.h" #define max2(a, b) ((a) > (b)? (a) : (b)) max2 List 1-2 List 1-2 chap01/max2x1test.c max2 #include
ohp03.dvi
19 3 ( ) 2019.4.20 CS 1 (comand line arguments) Unix./a.out aa bbb ccc ( ) C main void int main(int argc, char *argv[]) {... 2 (2) argc argv argc ( ) argv (C char ) ( 1) argc 4 argv NULL. / a. o u t \0
r07.dvi
19 7 ( ) 2019.4.20 1 1.1 (data structure ( (dynamic data structure 1 malloc C free C (garbage collection GC C GC(conservative GC 2 1.2 data next p 3 5 7 9 p 3 5 7 9 p 3 5 7 9 1 1: (single linked list 1
ohp07.dvi
19 7 ( ) 2019.4.20 1 (data structure) ( ) (dynamic data structure) 1 malloc C free 1 (static data structure) 2 (2) C (garbage collection GC) C GC(conservative GC) 2 2 conservative GC 3 data next p 3 5
II ( ) prog8-1.c s1542h017%./prog8-1 1 => 35 Hiroshi 2 => 23 Koji 3 => 67 Satoshi 4 => 87 Junko 5 => 64 Ichiro 6 => 89 Mari 7 => 73 D
II 8 2003 11 12 1 6 ( ) prog8-1.c s1542h017%./prog8-1 1 => 35 Hiroshi 2 => 23 Koji 3 => 67 Satoshi 4 => 87 Junko 5 => 64 Ichiro 6 => 89 Mari 7 => 73 Daisuke 8 =>. 73 Daisuke 35 Hiroshi 64 Ichiro 87 Junko
tuat1.dvi
( 1 ) http://ist.ksc.kwansei.ac.jp/ tutimura/ 2012 6 23 ( 1 ) 1 / 58 C ( 1 ) 2 / 58 2008 9 2002 2005 T E X ptetex3, ptexlive pt E X UTF-8 xdvi-jp 3 ( 1 ) 3 / 58 ( 1 ) 4 / 58 C,... ( 1 ) 5 / 58 6/23( )
kiso2-09.key
座席指定はありません 計算機基礎実習II 2018 のウェブページか 第9回 ら 以下の課題に自力で取り組んで下さい 計算機基礎実習II 第7回の復習課題(rev07) 第9回の基本課題(base09) 第8回試験の結果 中間試験に関するコメント コンパイルできない不完全なプログラムなど プログラミングに慣れていない あるいは複雑な問題は 要件 をバラして段階的にプログラムを作成する exam08-2.c
Gauss Strassen LU LU LU LU 22 5 Gauss LU
I, 208 3 23 3 2 4 3 6 3 6 32 6 32 6 322 7 323 Gauss 8 324 Strassen 9 325 0 326 0 327 LU 0 33 4 LU 2 4 2 4 2 42 3 43 4 42 8 43 0 LU 20 44 LU 22 5 Gauss LU 23 5 23 52 Gauss 23 53 LU 24 53 24 532 25 533 LU
C 2 / 21 1 y = x 1.1 lagrange.c 1 / Laglange / 2 #include <stdio.h> 3 #include <math.h> 4 int main() 5 { 6 float x[10], y[10]; 7 float xx, pn, p; 8 in
![C 2 / 21 1 y = x 1.1 lagrange.c 1 / Laglange / 2 #include <stdio.h> 3 #include <math.h> 4 int main() 5 { 6 float x[10], y[10]; 7 float xx, pn, p; 8 in](/thumbs/88/114913715.jpg "C 2 / 21 1 y = x 1.1 lagrange.c 1 / Laglange / 2 #include <stdio.h> 3 #include <math.h> 4 int main() 5 { 6 float x[10], y[10]; 7 float xx, pn, p; 8 in")
C 1 / 21 C 2005 A * 1 2 1.1......................................... 2 1.2 *.......................................... 3 2 4 2.1.............................................. 4 2.2..............................................
I. Backus-Naur BNF S + S S * S S x S +, *, x BNF S (parse tree) : * x + x x S * S x + S S S x x (1) * x x * x (2) * + x x x (3) + x * x + x x (4) * *
2015 2015 07 30 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF S + S S * S S x S +, *, x BNF S (parse tree) : * x + x x S * S x + S S S x x (1) * x x * x (2) * + x x x (3) +
OHP.dvi
0 7 4 0000 5.. 3. 4. 5. 0 0 00 Gauss PC 0 Gauss 3 Gauss Gauss 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 u [] u [3] u [4] u [4] P 0 = P 0 (),3,4 (,), (3,), (4,) 0,,,3,4 3 3 3 3 4 4 4 4 0 3 6 6 0 6 3 6 0 6
I. Backus-Naur BNF : N N 0 N N N N N N 0, 1 BNF N N 0 11 (parse tree) 11 (1) (2) (3) (4) II. 0(0 101)* (
2016 2016 07 28 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF : 11011 N N 0 N N 11 1001 N N N N 0, 1 BNF N N 0 11 (parse tree) 11 (1) 1100100 (2) 1111011 (3) 1110010 (4) 1001011
/* do-while */ #include <stdio.h> #include <math.h> int main(void) double val1, val2, arith_mean, geo_mean; printf( \n ); do printf( ); scanf( %lf, &v
1 http://www7.bpe.es.osaka-u.ac.jp/~kota/classes/jse.html kota@fbs.osaka-u.ac.jp /* do-while */ #include #include int main(void) double val1, val2, arith_mean, geo_mean; printf( \n );
G1. tateyama~$ gcc -c xxxxx.c ( ) xxxxx.o tateyama~$ gcc -o xxxxx.o yyyyy.o..... zzzzz.o Makefile make Makefile : xxxxx.o yyyyy.o... zzzzz.o ; gcc -o
G1. tateyama~$ gcc -c xxxxx.c ( ) xxxxx.o tateyama~$ gcc -o xxxxx.o yyyyy.o..... zzzzz.o Makefile make Makefile : xxxxx.o yyyyy.o... zzzzz.o ; gcc -o xxxxx.o yyyyy.o... zzzzz.o [1] [5] 1 [1] (matrix multi
ex01.dvi
,. 0. 0.0. C () /******************************* * $Id: ex_0_0.c,v.2 2006-04-0 3:37:00+09 naito Exp $ * * 0. 0.0 *******************************/ #include int main(int argc, char **argv) { double
1 return main() { main main C 1 戻り値の型 関数名 引数 関数ブロックをあらわす中括弧 main() 関数の定義 int main(void){ printf("hello World!!\n"); return 0; 戻り値 1: main() 2.2 C main
C 2007 5 29 C 1 11 2 2.1 main() 1 FORTRAN C main() main main() main() 1 return 1 1 return main() { main main C 1 戻り値の型 関数名 引数 関数ブロックをあらわす中括弧 main() 関数の定義 int main(void){ printf("hello World!!\n"); return
Original : Hello World! (0x0xbfab85e0) Copy : Hello World! (0x0x804a050) fgets mstrcpy malloc mstrcpy (main ) mstrcpy malloc free fgets stream 1 ( \n
2008 3 10 1 mstrcpy char *mstrcpy(const char *src); mstrcpy malloc (main free ) stdio.h fgets char *fgets(char *s, int size, FILE *stream); s size ( ) stream FILE ( man ) 40 ( ) %./a.out String : test
Microsoft Word - Cプログラミング演習(12)
第 12 回 (7/9) 4. いくつかのトピック (5)main 関数の引数を利用したファイル処理 main 関数は, 起動する環境から引数を受け取ることができる 例えば 次に示すように,main 関数に引数を用いたプログラムを作成する 01 /* sample */ 02 /* main 関数の引数 */ 03 #include 04 05 main(int argc, char
A 30 A A ( ) 2 C C (, machine language) C (C compiler) ( ) Mac Apple Xcode Clan
C 2017 9 29, 30 5 13 http://nalab.mind.meiji.ac.jp/~mk/labo/text/ 1 2 2 C 2 3 4 3.1 C................................... 4 3.2 Hello world........................................ 5 3.3 5...............................
r03.dvi
19 ( ) 019.4.0 CS 1 (comand line arguments) Unix./a.out aa bbb ccc ( ) C main void... argc argv argc ( ) argv (C char ) ( 1) argc 4 argv NULL. / a. o u t \0 a a \0 b b b \0 c c c \0 1: // argdemo1.c ---
A/B (2010/10/08) Ver kurino/2010/soft/soft.html A/B
A/B (2010/10/08) Ver. 1.0 kurino@math.cst.nihon-u.ac.jp http://edu-gw2.math.cst.nihon-u.ac.jp/ kurino/2010/soft/soft.html 2010 10 8 A/B 1 2010 10 8 2 1 1 1.1 OHP.................................... 1 1.2.......................................
:30 12:00 I. I VI II. III. IV. a d V. VI
2017 2017 08 03 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF X [ S ] a S S ; X X X, S [, a, ], ; BNF X (parse tree) (1) [a;a] (2) [[a]] (3) [a;[a]] (4) [[a];a] : [a] X 2 222222
ex01.dvi
,. 0. 0.0. C () /******************************* * $Id: ex_0_0.c,v.2 2006-04-0 3:37:00+09 naito Exp $ * * 0. 0.0 *******************************/ #include int main(int argc, char **argv) double
6 6.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave t T = N t N 44.1 khz t = 1 sec j t f j {f 0, f 1, f 2,, f N 1
6 6.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave t T = t 44.1 khz t = 1 sec 44100 j t f j {f 0, f 1, f 2,, f 1 6.2 T {f 0, f 1, f 2,, f 1 T ft) f j = fj t) j = 0, 1, 2,,
Gauss
15 1 LU LDL T 6 : 1g00p013-5 1 6 1.1....................................... 7 1.2.................................. 8 1.3.................................. 8 2 Gauss 9 2.1.....................................
£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (2018) - Â裶²ó ¨¡ À©¸æ¹½Â¤¡§·«¤êÊÖ¤· ¨¡
(2018) 2018 5 24 ( ) while ( ) do while ( ); for ( ; ; ) while int i = 0; while (i < 100) { printf("i = %3d\n", i); i++; while int i = 0; i while (i < 100) { printf("i = %3d\n", i); i++; while int i =
:30 12:00 I. I VI II. III. IV. a d V. VI
2018 2018 08 02 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF N N y N x N xy yx : yxxyxy N N x, y N (parse tree) (1) yxyyx (2) xyxyxy (3) yxxyxyy (4) yxxxyxxy N y N x N yx
実際の株価データを用いたオプション料の計算
2002 2 20 1 1 3 2 3 2.1 : : : : : : : : : : : : : : : : : : : : : : : : : : : : 5 2.1.1 : : : : : : : : : : : : : : : : : : : : 5 2.1.2 : : : : : : : : : : : : : : : : : : : : 6 2.2 : : : : : : : : : :
p = 1, 2, cos 2n + p)πj = cos 2nπj 2n + p)πj, sin = sin 2nπj 7.1) f j = a ) 0 + a p + a n+p cos 2nπj p=1 p=0 1 + ) b n+p p=0 sin 2nπj 1 2 a 0 +
7 7.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave T > 0 t 44.1 khz t = 1 44100 j t f j {f 0, f 1, f 2,, f 1 = T t 7.2 T {f 0, f 1, f 2,, f 1 T ft) f j = fj t) j = 0, 1,
(K&R 2.9) ~, &,, >>, << 2. (K&R 5.7) 3. (K&R 5.9) 4. (K&R 5.10) (argc argv atoi(), atof() ) 5. (K&R 7.5) (K&R 7.6) - FILE, stdin, stdout, std
3 3 1. 1 (K&R.9) ~, &,, >>,
file"a" file"b" fp = fopen("a", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose(fp); fp = fopen("b", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose
I117 9 2 School of Information Science, Japan Advanced Institute of Science and Technology file"a" file"b" fp = fopen("a", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose(fp); fp = fopen("b", "r"); while(fgets(line,
£Ã¥×¥í¥°¥é¥ß¥ó¥°(2018) - Âè11²ó – ½ÉÂꣲ¤Î²òÀ⡤±é½¬£² –
(2018) 11 2018 12 13 2 g v dv x dt = bv x, dv y dt = g bv y (1) b v 0 θ x(t) = v 0 cos θ ( 1 e bt) (2) b y(t) = 1 ( v 0 sin θ + g ) ( 1 e bt) g b b b t (3) 11 ( ) p14 2 1 y 4 t m y > 0 y < 0 t m1 h = 0001
joho07-1.ppt
0xbffffc5c 0xbffffc60 xxxxxxxx xxxxxxxx 00001010 00000000 00000000 00000000 01100011 00000000 00000000 00000000 xxxxxxxx x y 2 func1 func2 double func1(double y) { y = y + 5.0; return y; } double func2(double*
yacc.dvi
2017 c 8 Yacc Mini-C C/C++, yacc, Mini-C, run,, Mini-C 81 Yacc Yacc, 1, 2 ( ), while ::= "while" "(" ")" while yacc 1: st while : lex KW WHILE lex LPAREN expression lex RPAREN statement 2: 3: $$ = new
Informatics 2014
C 計算機の歴史 手回し計算機 新旧のソロバン バベッジの階差機関 スパコン ENIAC (1946) パソコン 大型汎用計算機 電卓 現在のコンピュータ Input Output Device Central Processing Unit I/O CPU Memory OS (Operating System) OS Windows 78, Vista, XP Windows Mac OS X
<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
Visual Basic でわかるやさしい有限要素法の基礎 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/092001 このサンプルページの内容は, 初版 1 刷発行当時のものです. URL http://www.morikita.co.jp/soft/92001/ horibe@mx.ibaraki.ac.jp
[ 1] 1 Hello World!! 1 #include <s t d i o. h> 2 3 int main ( ) { 4 5 p r i n t f ( H e l l o World!! \ n ) ; 6 7 return 0 ; 8 } 1:
![[ 1] 1 Hello World!! 1 #include <s t d i o. h> 2 3 int main ( ) { 4 5 p r i n t f ( H e l l o World!! \ n ) ; 6 7 return 0 ; 8 } 1:](/thumbs/97/132863283.jpg "[ 1] 1 Hello World!! 1 #include <s t d i o. h> 2 3 int main ( ) { 4 5 p r i n t f ( H e l l o World!! \ n ) ; 6 7 return 0 ; 8 } 1:")
005 9 7 1 1.1 1 Hello World!! 5 p r i n t f ( H e l l o World!! \ n ) ; 7 return 0 ; 8 } 1: 1 [ ] Hello World!! from Akita National College of Technology. 1 : 5 p r i n t f ( H e l l o World!! \ n ) ;
‚æ4›ñ
( ) ( ) ( ) A B C D E F G H I J K L M N O P Q R S T U V W X Y Z a b c d e f g h i j k l m n o p q r s t u v w x y z 0 1 2 3 4 5 6 7 8 9 (OUS) 9 26 1 / 28 ( ) ( ) ( ) A B C D Z a b c d z 0 1 2 9 (OUS) 9
2 1 list0415.c forfor #include int i, j; for (i = 1; i スライド 1
数値解析 平成 24 年度前期第 13 週 [7 月 11 日 ] 静岡大学創造科学技術大学院情報科学専攻工学部機械工学科計測情報講座 三浦憲二郎 講義アウトライン [7 月 11 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ラグランジュ補間 形状処理工学の基礎 点列からの曲線の生成 T.Kanai, U.Tokyo 関数近似 p.116 複雑な関数を簡単な関数で近似する関数近似 閉区間
sim98-8.dvi
8 12 12.1 12.2 @u @t = @2 u (1) @x 2 u(x; 0) = (x) u(0;t)=u(1;t)=0fort 0 1x, 1t N1x =1 x j = j1x, t n = n1t u(x j ;t n ) Uj n U n+1 j 1t 0 U n j =1t=(1x) 2 = U n j+1 0 2U n j + U n j01 (1x) 2 (2) U n+1
2 P.S.P.T. P.S.P.T. wiki 26
P.S.P.T. C 2011 4 10 2 P.S.P.T. P.S.P.T. wiki p.s.p.t.since1982@gmail.com http://www23.atwiki.jp/pspt 26 3 2 1 C 8 1.1 C................................................ 8 1.1.1...........................................
‚æ2›ñ C„¾„ê‡Ìš|
I 8 10 10 I ( 6 ) 10 10 1 / 23 1 C ( ) getchar(), gets(), scanf() ( ) putchar(), puts(), printf() 1 getchar(), putchar() 1 I ( 6 ) 10 10 2 / 23 1 (getchar 1 1) 1 #include 2 void main(void){ 3 int
スライド 1
数値解析 平成 29 年度前期第 14 週 [7 月 10 日 ] 静岡大学工学研究科機械工学専攻ロボット 計測情報分野創造科学技術大学院情報科学専攻 三浦憲二郎 期末試験 7 月 31 日 ( 月 ) 9 10 時限 A : 佐鳴会議室 B : 佐鳴ホール 講義アウトライン [7 月 10 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ( 復習 ) ラグランジュ補間 形状処理工学の基礎
I
I 998 208 0 9 Dirichlet ( ) 5 5 2 6 2 6 22 6 23 6 3 Dirichlet 7 3 7 32 7 4 Neumann 9 4 9 42 0 43 44 4 45 Neumann 4 5 6 5 Robin 6 52 7 53 Robin 22 6 23 6 23 62 24 2 2 25 2 Fourier 25 22 27 3 2 32 3 Target
BW BW
Induced Sorting BW 11T2042B 2015 3 23 1 1 1.1................................ 1 1.2................................... 1 2 BW 1 2.1..................................... 2 2.2 BW.................................
演習課題No12
演習課題 No.12 ( 課題は 3 題ある ) 課題 12-1 時間内提出 従来の C 言語には複素数を直接扱うデータ型はないので (*), 構造体で複素数 ( 英語で complex) を表すことにする. 複素数を表す構造体を以下のように定義する. struct complex float r; // 実部 ( 英語で real) float i; // 虚部 ( 英語で imaginary)
Taro-リストⅠ（公開版）.jtd
0. 目次 1. 再帰的なデータ構造によるリストの表現 1. 1 リストの作成と表示 1. 1. 1 リストの先頭に追加する方法 1. 1. 2 リストの末尾に追加する方法 1. 1. 3 昇順を保存してリストに追加する方法 1. 2 問題 問題 1 問題 2-1 - 1. 再帰的なデータ構造によるリストの表現 リストは データの一部に次のデータの記憶場所を示す情報 ( ポインタという ) を持つ構造をいう
file:///D|/C言語の擬似クラス.txt
愛知障害者職業能力開発校 システム設計科 修了研究発表会報告書 題名 : C 言語の擬似クラス あらまし : C 言語でクラスを作れるという噂の真偽を確かめるために思考錯誤した まえがき : VC++ や Java その他オブジェクト指向の言語にはクラスが存在して クラスはオブジェクトの設計図である 手法 : C++ のクラスを解析して C++ のクラスを作成して C 言語に翻訳する class struct
ｽﾗｲﾄﾞ ﾀｲﾄﾙなし
ファイル入出力 (2) これまでのおさらい ( 入出力 ) これまでの入出力は 入力 scanf 出力 printf キーボードと画面 ( 端末 ) scanf/printf は 書式つき入出力 フォーマットを指定する 標準入出力を対象とする 何もしなければ 標準入出力は キーボードと画面 ストリームという考え方 ストリーム (stream) = データの列 キーボードから打つ文字列 画面に出力される文字列
memo
数理情報工学演習第一 C プログラミング演習 ( 第 5 回 ) 2015/05/11 DEPARTMENT OF MATHEMATICAL INFORMATICS 1 今日の内容 : プロトタイプ宣言 ヘッダーファイル, プログラムの分割 課題 : 疎行列 2 プロトタイプ宣言 3 C 言語では, 関数や変数は使用する前 ( ソースの上のほう ) に定義されている必要がある. double sub(int
2008 ( 13 ) C LAPACK 2008 ( 13 )C LAPACK p. 1
2008 ( 13 ) C LAPACK LAPACK p. 1 Q & A Euler http://phase.hpcc.jp/phase/mppack/long.pdf KNOPPIX MT (Mersenne Twister) SFMT., ( ) ( ) ( ) ( ). LAPACK p. 2 C C, main Asir ( Asir ) ( ) (,,...), LAPACK p.
r08.dvi
19 8 ( ) 019.4.0 1 1.1 (linked list) ( ) next ( 1) (head) (tail) ( ) top head tail head data next 1: NULL nil ( ) NULL ( NULL ) ( 1 ) (double linked list ) ( ) 1 next 1 prev 1 head cur tail head cur prev
Informatics 2015
C 計算機の歴史 新旧のソロバン バベッジの階差機関 19C前半 手回し計算機 19C後半 20C後半 スパコン 1960年代 ENIAC (1946) 大型汎用計算機 1950年代 1980年代 電卓 1964 パソコン 1970年代 現在のコンピュータ Input Output Device Central Processing Unit I/O CPU Memory OS (Operating
untitled
II 4 Yacc Lex 2005 : 0 1 Yacc 20 Lex 1 20 traverse 1 %% 2 [0-9]+ { yylval.val = atoi((char*)yytext); return NUM; 3 "+" { return + ; 4 "*" { return * ; 5 "-" { return - ; 6 "/" { return / ; 7 [ \t] { /*
ohp08.dvi
19 8 ( ) 2019.4.20 1 (linked list) ( ) next ( 1) (head) (tail) ( ) top head tail head data next 1: 2 (2) NULL nil ( ) NULL ( NULL ) ( 1 ) (double linked list ) ( 2) 3 (3) head cur tail head cur prev data
A common.h include #include <stdio.h> #include <time.h> #define MAXN int A[MAXN], n; double start,end; void inputdata(
![A common.h include #include <stdio.h> #include <time.h> #define MAXN int A[MAXN], n; double start,end; void inputdata(](/thumbs/99/142484560.jpg "A common.h include #include <stdio.h> #include <time.h> #define MAXN int A[MAXN], n; double start,end; void inputdata(")
2 065762A 19 7 13 1 2 2.1 common.h include #include #include #define MAXN 1000000 int A[MAXN], n; double start,end; void inputdata(void) int i; // printf(" "); scanf("%d",&n); // printf("
memo
計数工学プログラミング演習 ( 第 3 回 ) 2016/04/26 DEPARTMENT OF MATHEMATICAL INFORMATICS 1 内容 ポインタ malloc 構造体 2 ポインタ あるメモリ領域 ( アドレス ) を代入できる変数 型は一致している必要がある 定義時には値は不定 ( 何も指していない ) 実際にはどこかのメモリを指しているので, #include
ex12.dvi
1 0. C, char., char, 0,. C, ("),., char str[]="abc" ; str abc.,, str 4. str 3. char str[10]="abc" ;, str 10, str 3., char s[]="abc", t[10] ;, t = s. ASCII, 0x00 0x7F, char., "abc" 3, 1. 1 8 256, 2., 2
A/B (2018/10/19) Ver kurino/2018/soft/soft.html A/B
A/B (2018/10/19) Ver. 1.0 kurino@math.cst.nihon-u.ac.jp http://edu-gw2.math.cst.nihon-u.ac.jp/ kurino/2018/soft/soft.html 2018 10 19 A/B 1 2018 10 19 2 1 1 1.1 OHP.................................... 1
プログラミング方法論 II 第 14,15 回 ( 担当 : 鈴木伸夫 ) 問題 17. x 座標と y 座標をメンバに持つ構造体 Point を作成せよ 但し座標 は double 型とする typedef struct{ (a) x; (b) y; } Point; 問題 18. 問題 17 の
プログラミング方法論 II 第 14,15 回 ( 担当 : 鈴木伸夫 ) 問題 17. x 座標と y 座標をメンバに持つ構造体 Point を作成せよ 但し座標 は double 型とする typedef struct{ (a) x; (b) y; Point; 問題 18. 問題 17 の Point を用いて 2 点の座標を入力するとその 2 点間の距 離を表示するプログラムを作成せよ 平方根は
Ｃ言語による数値計算プログラミング演習
5. 行列の固有値問題 n n 正方行列 A に対する n 個の固有値 λ i (i=1,,,n) と対応する固有ベクトル u i は次式を満たす Au = λ u i i i a11 a1 L a1 n u1i a1 a a n u i A =, ui = M O M M an 1 an L ann uni これらはまとめて, つぎのように書ける 5.1 ヤコビ法 = Λ, = [ u1 u u
2004 2005 2 2 1G01P038-0 1 2 1.1.............................. 2 1.2......................... 2 1.3......................... 3 2 4 2.1............................ 4 2.2....................... 4 2.3.......................
void hash1_init(int *array) int i; for (i = 0; i < HASHSIZE; i++) array[i] = EMPTY; /* i EMPTY */ void hash1_insert(int *array, int n) if (n < 0 n >=
![void hash1_init(int *array) int i; for (i = 0; i < HASHSIZE; i++) array[i] = EMPTY; /* i EMPTY */ void hash1_insert(int *array, int n) if (n < 0 n >=](/thumbs/89/97513163.jpg "void hash1_init(int *array) int i; for (i = 0; i < HASHSIZE; i++) array[i] = EMPTY; /* i EMPTY */ void hash1_insert(int *array, int n) if (n < 0 n >=")
II 14 2018 7 26 : : proen@mm.ics.saitama-u.ac.jp 14,, 8 2 12:00 1 O(1) n O(n) O(log n) O(1) 32 : 1G int 4 250 M 2.5 int 21 2 0 100 0 100 #include #define HASHSIZE 100 /* */ #define NOTFOUND 0
Informatics 2010.key
http://math.sci.hiroshima-u.ac.jp/ ~ryo/lectures/informatics2010/ 1 2 C ATM etc. etc. (Personal Computer) 3 4 Input Output Device Central Processing Unit I/O CPU Memory 5 6 (CPU),,... etc. C, Java, Fortran...
スライド 1
数値解析 2019 年度前期第 13 週 [7 月 11 日 ] 静岡大学創造科学技術大学院情報科学専攻工学部機械工学科計測情報講座 三浦憲二郎 講義アウトライン [7 月 11 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ラグランジュ補間 T.Kanai, U.Tokyo 関数近似 p.116 複雑な関数を簡単な関数で近似する 関数近似 閉区間 [a,b] で定義された関数 f(x)
1 1.1 C 2 1 double a[ ][ ]; 1 3x x3 ( ) malloc() malloc 2 #include <stdio.h> #include
![1 1.1 C 2 1 double a[ ][ ]; 1 3x x3 ( ) malloc() malloc 2 #include <stdio.h> #include](/thumbs/98/138365120.jpg "1 1.1 C 2 1 double a[ ][ ]; 1 3x x3 ( ) malloc() malloc 2 #include <stdio.h> #include")
1 1.1 C 2 1 double a[ ][ ]; 1 3x3 0 1 3x3 ( ) 0.240 0.143 0.339 0.191 0.341 0.477 0.412 0.003 0.921 1.2 malloc() malloc 2 #include #include #include enum LENGTH = 10 ; int
超初心者用
3 1999 10 13 1. 2. hello.c printf( Hello, world! n ); cc hello.c a.out./a.out Hello, world printf( Hello, world! n ); 2 Hello, world printf n printf 3. ( ) int num; num = 100; num 100 100 num int num num
1.ppt
/* * Program name: hello.c */ #include int main() { printf( hello, world\n ); return 0; /* * Program name: Hello.java */ import java.io.*; class Hello { public static void main(string[] arg)
P02.ppt
int If 2 1 ,,, 3 a + b ab a - b ab a * b ab a / b ab a % b ab a + b 4 2 list0201.c /, % /*/ int vx, vy; puts(""); printf("vx"); scanf("%d", &vx); printf("vy"); scanf("%d", &vy); printf("vx + vy = %d\n",
C¥×¥í¥°¥é¥ß¥ó¥° ÆþÌç
C (3) if else switch AND && OR (NOT)! 1 BMI BMI BMI = 10 4 [kg]) ( [cm]) 2 bmi1.c Input your height[cm]: 173.2 Enter Input your weight[kg]: 60.3 Enter Your BMI is 20.1. 10 4 = 10000.0 1 BMI BMI BMI = 10
273? C
TSG Theoretical Science Group 273? C 2-1.................................. 2 -1. Windows Mac Mac UNIX CUI bash >_ Finder TSG No.273? 2 3 pwd ls cd ( ) change directory 3 TSG No.273? cd hoge cd hoge cd....../
17 1. strucr Potter ( ) Harry Potter and the Philosopher s Stone 1997 English Title : Harry Potter and the Philosopher s Stone : : 1997 #include<stdio
17 1. strucr Potter ( ) Harry Potter and the Philosopher s Stone 1997 English Title : Harry Potter and the Philosopher s Stone : : 1997 #include typedef struct Potter{ Potter; int main(void){
8 / 0 1 i++ i 1 i-- i C !!! C 2
C 2006 5 2 printf() 1 [1] 5 8 C 5 ( ) 6 (auto) (static) 7 (=) 1 8 / 0 1 i++ i 1 i-- i 1 2 2.1 C 4 5 3 13!!! C 2 2.2 C ( ) 4 1 HTML はじめ mkdir work 作業用ディレクトリーの作成 emacs hoge.c& エディターによりソースプログラム作成 gcc -o fuga
Microsoft PowerPoint - kougi9.ppt
C プログラミング演習 第 9 回ポインタとリンクドリストデータ構造 1 今まで説明してきた変数 #include "stdafx.h" #include int _tmain(int argc, _TCHAR* argv[]) { double x; double y; char buf[256]; int i; double start_x; double step_x; FILE*
(300, 150) 120 getchar() HgBox(x, y, w, h) (x, y), w, h #include <stdio.h> #include <handy.h> int main(void) { int i; double w, h; } HgO
Handy Graphic for Handy Graphic Version 0.5 2008-06-09 1 Handy Graphic Handy Graphic C Handy Graphic Handy Graphic Mac OS X Handy Graphic HgDisplayer Handy Graphic HgDisplayer 2 Handy Graphic 1 Handy Graphic
3.1 stdio.h iostream List.2 using namespace std C printf ( ) %d %f %s %d C++ cout cout List.2 Hello World! cout << float a = 1.2f; int b = 3; cout <<
C++ C C++ 1 C++ C++ C C++ C C++? C C++ C *.c *.cpp C cpp VC C++ 2 C++ C++ C++ [1], C++,,1999 [2],,,2001 [3], ( )( ),,2001 [4] B.W. /D.M.,, C,,1989 C Web [5], http://kumei.ne.jp/c_lang/ 3 Hello World Hello
2000 7 17 iii,,.,,,,.,.,,. =,,.,,.,.,,,,,,, fortran,.,,,,.,,, iv (i),., 10,, 10 1. (ii),, fortran, fortran, C,.. (i),., 10 23.. (ii),.,, fortran Pasal,, for, if then else. C., fortran C.,.,,.,.,,.,,.,.,,,,.
untitled
Q 8 1 8.1 (C++) C++ cin cout 5 C++ 16 6 p.63 8.3 #include 7 showbase noshowbase showpoint noshowpoint 8.3 uppercase 16 nouppercase 16 setfill(int) setprecision(int) setw(int) setbase(int) dec
C
C 1 2 1.1........................... 2 1.2........................ 2 1.3 make................................................ 3 1.4....................................... 5 1.4.1 strip................................................