function structure = cbg2ese(system_name, system_type, system_cr, ...
system_args, full_name, full_name_repetition, ...
repetition,...
structure, structure_file,infofilenum)
% Set up globals to count the component inputs and outputs. This relies on
% the named SS (the ports) being in the correct order. Using globals here
% avoids changing the common argument list for all _eqn files for something
% which is only used for named SS components.
global local_u_index
global local_y_index
global at_top_level
%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %%%%% Model Transformation Tools %%%%%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Matlab function cbg2ese.m
% Acausal bond graph to causal bond graph: mfile format
% Structure matrix [states,nonstates,inputs,outputs,zero_outputs]
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.24 1998/07/27 10:26:02 peterg
% %% No change - but fixed bug in alias_args
% %%
% %% Revision 1.23 1998/07/08 12:33:51 peterg
% %% Reinstated the infofilenum parameter.
% %%
% %% Revision 1.22 1998/07/04 07:10:27 peterg
% %% Don't evaluate alias if no constitutive relationship or args and write
% %% message.
% %%
% %% Revision 1.21 1998/07/03 18:58:58 peterg
% %% Put arg alias stuff within function alias_args
% %% Called recursively to handle arithmetic expressions
% %%
% %% Revision 1.20 1998/07/03 14:39:09 peterg
% %% Added info messages a bit busy now!
% %%
% %% Revision 1.19 1998/04/12 11:58:19 peterg
% %% Rename port components by changing name_r to [name_r
% %%
% %% Revision 1.18 1998/04/11 18:59:16 peterg
% %% at_top_level now global - passed to SS components
% %%
% %% Revision 1.17 1998/04/04 10:47:31 peterg
% %% Uses (coerced) components from _cbg file - _abg not now used here.
% %%
% %% Revision 1.16 1998/03/06 09:38:58 peterg
% %% Now writes correct structure file for multiport C & I
% %%
% %% Revision 1.15 1997/12/16 18:24:33 peterg
% %% Added unknown_input to structure list
% %%
% %% Revision 1.14 1997/08/26 07:51:10 peterg
% %% Now counts the local input and outputs by order of appearence rather
% %% than by port number - it therfore handles ports with bicausality correctely.
% %%
% %% Revision 1.13 1997/05/19 16:45:56 peterg
% %% Fixed ISW bug -- deleted spurious ISW_eqn.m file
% %%
% %% Revision 1.12 1997/04/15 09:17:26 peterg
% %% Added the structure file - contains details of states etc.
% %%
% %% Revision 1.11 1996/12/07 18:20:11 peterg
% %% Replaces null argument by a default and tells user.
% %%
% %% Revision 1.10 1996/12/07 17:37:07 peterg
% %% Now handles numbered SS ports appearing at top level.
% %%
% %% Revision 1.9 1996/12/04 21:49:47 peterg
% %% Compares full-name with empty string (instead of testing for zero
% %% length)
% %%
% %% Revision 1.8 1996/08/30 16:36:08 peter
% %% More info written to ese files.
% %%
% %% Revision 1.7 1996/08/30 11:23:13 peter
% %% Argument substitution implemented.
% %%
% %% Revision 1.6 1996/08/27 08:04:52 peterg
% %% Handles complex components and repetative components.
% %%
% %% Revision 1.5 1996/08/24 15:02:23 peter
% %% Writes `END;' to keep reduce happy.
% %%
% %% Revision 1.4 1996/08/19 09:03:41 peter
% %% Handles repeating components.
% %%
% %% Revision 1.3 1996/08/18 20:08:02 peter
% %% Included additional structure: structure(5) = zero_outputs.
% %%
% %% Revision 1.2 1996/08/08 18:08:11 peter
% %% Sorted out file naming sceme
% %%
% %% Revision 1.1 1996/08/08 15:53:23 peter
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
system_name, system_type, full_name,
pc = "%";
% Set up the names corresponding to the structure matrix.
structure_name = [
"state ",
"nonstate ",
"input ",
"output ",
"zero_output ",
"unknown_input"];
% Are we at the top level of the heirarchy?
at_top_level = strcmp(full_name, "");
% Create the (full) system name
if at_top_level
full_name = system_name;
full_name_repetition = system_name;
system_type = system_name;
else
full_name = [full_name, "_", system_name];
full_name_repetition = [full_name_repetition, ...
"_", system_name, "_", num2str(repetition)];
end;
% $$$ full_name_type = [full_name, "_", system_type];
% $$$ cbg_name = [full_name_type, "_cbg"];
cbg_name = [full_name, "_cbg"];
% Return if cbg file doesn't exist
if exist(cbg_name)~=2
disp([cbg_name, " does not exist"]);
return
end;
% Setup files
ese_name = [full_name_repetition, "_ese.r"];
ese_file = fopen(ese_name, "w");
fprintf(ese_file, "\n%s%s Equation file for system %s (file %s)\n", ...
pc, pc, full_name_repetition, ese_name);
fprintf(ese_file, "%s%s Generated by MTT\n\n", pc, pc);
% Evaluate the system function to get the bonds
eval(["[bonds,status,system_type,components] = ", cbg_name, ";"]);
abg_name = [system_type, "_abg"];
cmp_name = [system_type, "_cmp"];
alias_name = [system_type, "_alias"];
% No longer needed - cbg now has components:
% eval(["[junk,components]=", abg_name, ";"]);
% Find number of bonds
[n_bonds,columns] = size(bonds);
if (columns ~= 2)&(n_bonds>0)
error("Incorrect bonds matrix: must have 2 columns");
end;
% Find number of components
[n_components,columns] = size(components);
n_components = n_components
% Set up the first dummy bond number - needed for repetative components
next_bond = max(max(abs(components)))+1;
% Set up the counters for the labelled SS. These are, by definition,
% encountered first and so the counters will not be messed up by subsystems.
local_u_index = 0;
local_y_index = 0;
if (length(system_args)==0)
mtt_info(sprintf("No arguments given so no argument aliasing done for system %s(%s)",\
system_name,system_type), infofilenum);
else
eval([ "alias = ", alias_name ";"]);
end;
if (length(system_cr)==0)
mtt_info(sprintf("No cr given so no cr aliasing done for system %s(%s)",\
system_name,system_type), infofilenum);
else
eval([ "alias = ", alias_name ";"]);
end;
for i = 1:n_components
comp = nozeros(components(i,:));
bond_list = abs(comp);
direction = sign(comp)'*[1 1];
% Convert from arrow orientated to component orientated causality
comp_bonds = bonds(bond_list,:).*direction;
disp("---- Component ---");
% Get the component details
eval([ "[comp_type,comp_name,cr,args,repetitions] = ", cmp_name, "(i)" ...
]);
% Alias the args list -- if not at top level
message = sprintf("\tfor component %s (%s) within %s",\
comp_name,comp_type,full_name);
if (length(system_args)>0)
args = alias_args(args,alias,";",message,infofilenum)
end;
if (length(system_cr)>0)
cr = alias_args(cr,alias,";",message,infofilenum)
end;
% Substitute positional ($1 etc) arguments
cr = subs_arg(cr,system_cr, ...
"lin",full_name,comp_type,comp_name,infofilenum);
args = subs_arg(args,system_args, ...
"1",full_name,comp_type,comp_name,infofilenum);
% change name of 0 and 1 components -- matlab doesn't like numbers here
if strcmp(comp_type,"0")
comp_type = "zero";
end;
if strcmp(comp_type,"1")
comp_type = "one";
end;
ports = length(bond_list);
if repetitions>1
port_pairs = ports/2;
if round(port_pairs)~=port_pairs;
mtt_info(["Repeated component ", comp_name, ...
" has an odd number of ports - ignoring repetitions"], infofilenum);
repetitions = 1;
end;
end;
if repetitions>1
odd_bonds = bond_list(1:2:ports-1);
even_bonds = bond_list(2:2:ports);
end;
for k = 1:repetitions
if repetitions>1
if k==1
bond_list(1:2:ports-1) = odd_bonds;
else
bond_list(1:2:ports-1) = bond_list(2:2:ports);
end;
if k==repetitions
bond_list(2:2:ports) = even_bonds;
else
new_bonds = [next_bond:next_bond+port_pairs-1];
next_bond = next_bond+port_pairs;
bond_list(2:2:ports) = new_bonds;
end;
end;
% Invoke the appropriate equation-generating procedure
name_r = full_name_repetition;
eqn_name = [comp_type, "_eqn"]
if exist(eqn_name)~=2 % Try a compound component
disp("---PUSH---"); bond_list, comp_name, comp_type
structure = cbg2ese(comp_name, comp_type, cr, args, ...
full_name, full_name_repetition, ...
k, structure, structure_file, infofilenum);
% Link up the bonds
fprintf(ese_file, ...
"\n\t%s Equations linking up subsystem %s (%s)\n\n", ...
pc, comp_name, comp_type);
name_comp_name = sprintf("%s_%s_%d", ...
full_name_repetition, comp_name, k);
printf("\n\t%s Equations linking up subsystem %s (%s)\n\n",\
pc, comp_name, comp_type);
bond_list
u_index = 0; y_index = 0; % Count the inputs and outputs
for port_number=1:length(bond_list)
% Effort
if comp_bonds(port_number,1)==1 % Source
u_index = u_index + 1;
fprintf(ese_file, "%s_MTTu%d := %s;\n", ...
name_comp_name, u_index, varname(name_r, ...
bond_list(port_number),1));
else % Sensor
y_index = y_index + 1;
fprintf(ese_file, "%s := %s_MTTy%d;\n", ...
varname(name_r, ...
bond_list(port_number),1), name_comp_name, y_index);
end;
% Flow
if comp_bonds(port_number,2)==-1 % Source
u_index = u_index + 1;
fprintf(ese_file, "%s_MTTu%d := %s;\n", ...
name_comp_name, u_index, varname(name_r, ...
bond_list(port_number),-1));
else % Sensor
y_index = y_index + 1;
fprintf(ese_file, "%s := %s_MTTy%d;\n", ...
varname(name_r, ...
bond_list(port_number),-1), name_comp_name, y_index);
end;
end;
disp("---POP---");
else % its a simple component
fprintf(ese_file, "\n\t%s Equations for component %s (%s), repetition %d\n\n", ...
pc, comp_name, comp_type,k);
% Is it a named port? (Name begins with [)
Named_Port = (comp_name(1)=="[");
if Named_Port %Add [ to start of name
name_r = ["[" name_r];
end;
# % Take port SS to be ordinary SS at top level
# if at_top_level & strcmp(comp_type, "SS")
# effort_attribute = cr;
# flow_attribute = args;
# if strcmp(effort_attribute, "MTT_port") % Its a numbered port
# effort_attribute = "external";
# flow_attribute = "external";
# cr = effort_attribute;
# args = flow_attribute;
# end;
# end;
% Save the current structure
old_structure = structure;
% Generate the simple component equations
eval(["structure = ", ...
eqn_name, ...
"(name_r,bond_list,comp_bonds, ...
direction,cr,args,structure,ese_file);" ]);
% If structure has changed, write info to structure file.
structure_change = structure-old_structure;
% The following line is to avoid probs with multiport C or I
% it needs changing to handle multi ports correctly
structure_changes = sum(structure_change);
structure_change = structure_change>zeros(size(structure_change));
if structure_changes>0
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%1.0f\t%s\t%s\t%1.0f\n", ...
structure_name(which_index,:), value-k+1, ...
comp_name, full_name, repetition);
end;
end;
end;
end;
end;
end;
% Close the files
fclose(ese_file);