19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
+
+
+
|
## Structure matrix [states,nonstates,inputs,outputs,zero_outputs]
## ###############################################################
## ## Version control history
## ###############################################################
## ## $Id$
## ## $Log$
## ## Revision 1.46 2002/08/20 15:51:17 gawthrop
## ## Update to work with ident DIY rep
## ##
## ## Revision 1.45 2002/05/22 09:15:03 gawthrop
## ## Non-repetitive components no longer use _1 in names
## ##
## ## Revision 1.44 2001/11/11 18:12:30 geraint
## ## Moved fflush(structure_file) out of loop.
## ##
## ## Revision 1.43 2001/11/11 08:32:00 geraint
|
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
|
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
|
-
+
-
+
+
|
which_indices = getindex(structure_change,1);
which_indices = which_indices(:,2)';
for which_index=which_indices
value = structure(which_index);
value_change=value-old_structure(which_index);
for k=1:value_change
fprintf(structure_file, ...
"%s\t%i\t%s\t%s_%s\t%i\n", ...
"%s\t%i\t%s\t%s_%s\t%i\t%s\t%s\n", ...
structure_name(which_index,:), value-k+1, ...
comp_name, full_name_repetition, comp_name, repetition);
comp_name, full_name_repetition, comp_name, \
repetition, cause2name(-comp_bonds(1)), full_name_repetition );
endfor;
endfor;
endif;
endif
endfor
fflush (structure_file);
|