312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
|
sed "s/$space\($non_space\)$spaces=$spaces\($fun_name\)(\([a-zA-Z0-9,]*\))/\2(\1,\3)/" |\
sed "s/$space\($non_space\)$spaces=$spaces\($fun_name\)/\2(\1)/" |\
awk '
function printvar(Name,N) {
if (N<1) return;
width = 10;
kk = 0;
printf("VAR\n");
for (k=1;k<=N;k++) {
printf("mtt%s%i", Name, k);
if (k==N)
printf(" : REAL;\n")
else
printf(", ");
kk++;
if (kk==width){printf("\n"); kk=0}
}
}
BEGIN{
comment_regexp = "{"
doing_header = 0
doing_globals = 0
Nt = 99
inc ="$I"
}
{
if ($1=="global")
doing_globals = 1
else{
if (doing_globals==1){
|
<
|
|
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
|
sed "s/$space\($non_space\)$spaces=$spaces\($fun_name\)(\([a-zA-Z0-9,]*\))/\2(\1,\3)/" |\
sed "s/$space\($non_space\)$spaces=$spaces\($fun_name\)/\2(\1)/" |\
awk '
function printvar(Name,N) {
if (N<1) return;
width = 10;
kk = 0;
for (k=1;k<=N;k++) {
printf("mtt%s%i", Name, k);
if (k==N)
printf(" : REAL;\n")
else
printf(", ");
kk++;
if (kk==width){printf("\n"); kk=0}
}
}
BEGIN{
comment_regexp = "{"
doing_header = 0
doing_globals = 0
Nt = 9
inc ="$I"
}
{
if ($1=="global")
doing_globals = 1
else{
if (doing_globals==1){
|
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
|
printf("{%s $MTTPATH/trans/p/sign.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_euler.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_solve.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_implicit.p}\n",inc)
# printf("{%s $MTTPATH/trans/p/zero_matrix.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_input.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_state.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_getargs.p}\n",inc)
printf("{%s %s_simpar.p}\n",inc,Sys)
printf("{%s %s_numpar.p}\n",inc,Sys)
printf("{%s %s_state.p}\n",inc,Sys)
printf("{%s %s_input.p}\n",inc,Sys)
if (Method=="euler") {
printf("{%s %s_ode.p}\n",inc,Sys)
printf("{%s %s_odeo.p}\n",inc,Sys)
}
if (Method=="implicit") {
printf("{%s %s_cse.p}\n",inc,Sys)
printf("{%s %s_cseo.p}\n",inc,Sys)
printf("\n{%s %s_smx.p}\n",inc,Sys);
}
printf("{%s %s_switchopen.p}\n\n",inc,Sys)
for (k=1;k<=j;k++) printf("%s\n", comment[k])
printf("\n")
printf("\nBEGIN{%s}\n", Sys_rep)
}
else{
for (k=1;k<=j;k++) printf("%s\n", comment[k])
printf("\n")
printf("VAR mmtti : INTEGER;");
printf("\nBEGIN{%s}\n", Sys_rep)
if (rep=="smxa"){
printf("FOR mmtti:= 1 to %i DO\n\tmtta[mmtti] := 0;\n", Nx*Nx);
}
if (rep=="smxax"){
printf("FOR mmtti:= 1 to %i DO\n\tmttax[mmtti] := 0;\n", Nx);
}
}
doing_header = 0;
}
if (match($1,comment_regexp)>0){
if (doing_header==1)
comment[++j] = $0
|
<
>
|
>
>
>
|
|
|
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
|
printf("{%s $MTTPATH/trans/p/sign.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_euler.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_solve.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_implicit.p}\n",inc)
# printf("{%s $MTTPATH/trans/p/zero_matrix.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_input.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_state.p}\n",inc)
printf("{%s %s_simpar.p}\n",inc,Sys)
printf("{%s %s_numpar.p}\n",inc,Sys)
printf("{%s $MTTPATH/trans/p/mtt_getargs.p}\n",inc)
printf("{%s %s_state.p}\n",inc,Sys)
printf("{%s %s_input.p}\n",inc,Sys)
if (Method=="euler") {
printf("{%s %s_ode.p}\n",inc,Sys)
printf("{%s %s_odeo.p}\n",inc,Sys)
}
if (Method=="implicit") {
printf("{%s %s_cse.p}\n",inc,Sys)
printf("{%s %s_cseo.p}\n",inc,Sys)
printf("\n{%s %s_smx.p}\n",inc,Sys);
}
printf("{%s %s_switchopen.p}\n\n",inc,Sys)
for (k=1;k<=j;k++) printf("%s\n", comment[k])
printf("\n")
printf("\nBEGIN{%s}\n", Sys_rep)
}
else{
for (k=1;k<=j;k++) printf("%s\n", comment[k])
printf("\n")
printf("VAR mtt_i : INTEGER;\n");
if ( (rep=="smxa")||(rep=="smxax") ) {
printvar("_t",Nt);
}
printf("\nBEGIN{%s}\n", Sys_rep)
if (rep=="smxa"){
printf("FOR mtt_i:= 1 to %i DO\n\tmtta[mtt_i] := 0;\n", Nx*Nx);
}
if (rep=="smxax"){
printf("FOR mtt_i:= 1 to %i DO\n\tmttax[mtt_i] := 0;\n", Nx);
}
}
doing_header = 0;
}
if (match($1,comment_regexp)>0){
if (doing_header==1)
comment[++j] = $0
|