function structure = SS_eqn(name,Bond_number,Bonds,Direction,cr,args, ... structure,filenum); disp("SS_eqn"); name,Bond_number,Bonds,Direction ## 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 ## SS_eqn - equations for SS component ## ## ####################################### ## ###### Model Transformation Tools ##### ## ###################################### ## ## Matlab function SS_eqn ## structure = SS_eqn(name,bond_number,bonds,direction,cr,args, ... ## structure,filenum); ## Copyright (c) P.J. Gawthrop, 1996. ## ############################################################### ## ### Version control history ## ############################################################### ## ### $Id$ ## ### $Log$ ## ### Revision 1.27 2005/03/21 10:34:47 gawthrop ## ### Tidied. ## ### ## ### Revision 1.26 2004/08/29 13:11:21 geraint ## ### Writes MTTyz vector in the same format as MTTx and MTTy: MTTyz($d,1) ## ## ## ### Revision 1.25 2003/01/07 18:09:21 geraint ## ### More informative error message. ## ## ## ### Revision 1.24 1999/10/22 02:20:32 peterg ## ### Changed error handleing of number of attributes >2 ## ## ## ### Revision 1.23 1998/12/03 16:46:16 peterg ## ### Deblanked attributes so that zero attribute works. ## ## ## ### Revision 1.22 1998/07/28 19:05:46 peterg ## ### Fixed a few bugs. ## ## ## ### Revision 1.21 1998/07/28 14:21:31 peterg ## ### Vector SS ports. ## ## ## ### Revision 1.20 1998/07/08 14:42:52 peterg ## ### Removed the annoying info message - it causes problems with big ## ### systems ## ## ## ### Revision 1.19 1998/07/08 11:33:54 peterg ## ### Replace mtt_info by mtt_error when appropriate ## ## ## ### Revision 1.18 1998/07/08 11:30:45 peterg ## ### Removed second (fileID) argument from mtt_info ## ## ## ### Revision 1.17 1998/07/04 07:15:44 peterg ## ### Back under RCS ## ## ## ### Revision 1.16 1998/04/11 19:07:16 peterg ## ### Now do named ports as ordinary ports iff at top level. ## ### --- not yet complete, need to pass necesssary info though to this ## ### function ## ## ## ### Revision 1.15 1997/12/16 19:16:07 peterg ## ### Added unknown input to the effort part. ## ## ## ### Revision 1.14 1997/12/16 18:25:19 peterg ## ### Added unknown_input attribure to flow -- effort still needs doing ## ## ## ### Revision 1.13 1997/09/18 13:15:15 peterg ## ### Fixed incorrect error message flagging inappropriate flow outputs ## ### -- used to give the effort rather than the flow in the error message. ## ## ## ### Revision 1.12 1997/08/26 07:51:30 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.11 1997/05/09 08:21:07 peterg ## ### Explicit computation of port number -- avoids str2num ## ## ## ### Revision 1.10 1997/03/22 17:13:03 peterg ## ### Fixed bug for port nos. > 1 digit! ## ## ## ### Revision 1.9 1997/03/22 15:50:59 peterg ## ### Changed #1.0f to %d format. ## ## ## ### Revision 1.8 1996/12/10 16:52:29 peterg ## ### Detect null string using strcmp, not length. ## ### Put filnum argument to mtt_info. ## ## ## ### Revision 1.7 1996/12/07 17:17:40 peterg ## ### Added some ; ## ## ## ### Revision 1.6 1996/12/05 09:49:09 peterg ## ### Explicit computation of port number from string. ## ## ## ### Revision 1.5 1996/12/04 21:27:53 peterg ## ### Replaced str2num by sprintf ## ## ## ### Revision 1.4 1996/08/18 20:06:21 peter ## ### Included zero outputs. ## ## ## ### Revision 1.3 1996/08/09 14:08:04 peter ## ### Empty effort and flow attributes replaced by 'external'. ## ## ## ### Revision 1.2 1996/08/08 16:38:19 peter ## ### Changed to filenumer type of fprintf ## ## ## ### Revision 1.1 1996/08/08 16:19:08 peter ## ### Initial revision ## ## ## ############################################################### ## The string 'cr' contains the effort attributes; ## 'arg' contains the flow attributes. ## external indicates an external variable ## internal indicates an internal variable ## anything else is the name of a constant. ## Blank is taken to be external ## If its a numbered port: ## cr contains 'MTT_port' ## arg contains port number if (strcmp(cr,"SS")) # Then its the standard file a = char(strsplit(args,",")) [N,M]=size(a); if (N~=2) # Must have 2 arguments mtt_error(sprintf("SS should have 2 args not %i", N)); endif effort_attribute = deblank(a(1,:)); flow_attribute = deblank(a(2,:)); else # Old style file effort_attribute = cr; flow_attribute = args; endif ## Default attributes if strcmp(effort_attribute,"") effort_attribute = "external"; endif if strcmp(flow_attribute,"") flow_attribute = "external"; endif inputs = structure(3); outputs = structure(4); zero_outputs = structure(5); unknown_inputs = structure(6); n_bonds = length(Bond_number) # Multi port? for i=1:n_bonds # Loop over all the bonds bond_number = Bond_number(i); bonds = Bonds(i,:); direction = Direction(i,:); Named_Port = name(1)=="["; if Named_Port sname = name(2:length(name)); #Strip the [] else sname = name; endif if Named_Port&&~at_top_level # It's a named port ## Effort if bonds(1,1)==-1 # Source local_u_index = local_u_index + 1 fprintf(filenum, "%s := %s_MTTu%d;\n", ... varname(sname, bond_number,1), sname, local_u_index); else # Sensor local_y_index = local_y_index + 1 fprintf(filenum, "%s_MTTy%d := %s;\n", ... sname, local_y_index, varname(sname, bond_number,1)); endif ## Flow if bonds(1,2)==1 # Source local_u_index = local_u_index + 1 fprintf(filenum, "%s := %s_MTTu%d;\n", ... varname(sname, bond_number,-1), sname, local_u_index); else # Sensor local_y_index = local_y_index + 1 fprintf(filenum, "%s_MTTy%d := %s;\n", ... sname, local_y_index, varname(sname, bond_number,-1)); endif else # Now do SS which are not ports ## Effort if strcmp(effort_attribute, "external") if bonds(1,1)==-1 # Source inputs = inputs+1; fprintf(filenum, "%s := MTTu(%d,1);\n", ... varname(sname, bond_number,1),inputs); else # Sensor outputs = outputs+1; fprintf(filenum, "MTTy(%d,1) := %s;\n", ... outputs, varname(sname, bond_number,1)); endif elseif strcmp(effort_attribute, "unknown") # Unknown input unknown_inputs = unknown_inputs + 1; fprintf(filenum, "%s := MTTUi%d;\n", ... varname(sname, bond_number,1), unknown_inputs); elseif strcmp(effort_attribute, "internal") ## Do nothing else if bonds(1,1)==-1 # Named or unknown source fprintf(filenum, "%s := %s;\n", ... varname(sname, bond_number,1), effort_attribute); else # Sensor if strcmp(effort_attribute, "zero") #Zero output zero_outputs = zero_outputs + 1; fprintf(filenum, "MTTyz(%d,1) := %s;\n", ... zero_outputs, varname(sname, bond_number,1)); else mtt_error([effort_attribute, " not appropriate for an output (", name, ")"]); endif endif endif ## Flow if strcmp(flow_attribute, "external") if bonds(1,2)==1 # Source inputs = inputs+1 fprintf(filenum, "%s := MTTu(%d,1);\n", varname(sname, bond_number,-1),inputs); else # Sensor outputs = outputs+1 fprintf(filenum, "MTTy(%d,1) := %s;\n", outputs, ... varname(sname, bond_number,-1)); endif elseif strcmp(flow_attribute, "unknown") # Unknown input unknown_inputs = unknown_inputs + 1 fprintf(filenum, "%s := MTTUi%d;\n", ... varname(sname, bond_number,-1), unknown_inputs); elseif strcmp(flow_attribute, "internal") ## Do nothing else # Named constant if bonds(1,2)==1 # Source fprintf(filenum, "%s := %s;\n", ... varname(sname, bond_number,-1), flow_attribute); else # Sensor if strcmp(flow_attribute, "zero") #Zero output zero_outputs = zero_outputs + 1 fprintf(filenum, "MTTyz(%d,1) := %s;\n", ... zero_outputs, varname(sname, bond_number,-1)); else mtt_error([flow_attribute, " not appropriate for an output (", name, ")"]); endif endif endif endif structure(3) = inputs; structure(4) = outputs; structure(5) = zero_outputs; structure(6) = unknown_inputs; endfor endfunction