ADDED   mttroot/mtt/bin/trans/m/abg2cbg.m
Index: mttroot/mtt/bin/trans/m/abg2cbg.m
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/m/abg2cbg.m
@@ -0,0 +1,135 @@
+function [bonds,status] = abg2cbg(system_name,bonds,infofile)
+% [bonds,status] = abg2cbg(system_name,bonds,infofile)
+% 
+%     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
+%     %%%%% Model Transformation Tools %%%%%
+%     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% 
+% Matlab function  abg2cbg.m
+% Acausal bond graph to causal bond graph: mfile format
+
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% %% Version control history
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% %% $Id$
+% %% $Log$
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+
+pc = '%';
+if nargin<3
+  infofile = 'stdout';
+end;
+
+% Evaluate the system function to get the bonds
+fun_name = [system_name, '_abg']
+if exist(fun_name)~=2
+  mtt_info(['m-file ', fun_name, ' does not exist'], infofile);
+  bonds = [];
+  status = [];
+else
+  eval(['[bonds,components]=', fun_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);
+
+  % Set initial status
+  status = -ones(n_components,1);
+  total = 2*n_bonds;
+  done = sum(sum(abs(bonds)))/total*100;
+  mtt_info(sprintf('Initial causality is %3.0f%s complete.', done, pc), infofile);
+
+  old_done = inf;
+
+  while done~=old_done
+    disp(sprintf('Causality is %3.0f%s complete.', done, pc));
+    old_done = done;
+
+    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;
+      eval([ '[comp_type,name,cr,arg] = ', system_name, '_cmp(i);' ]);
+      % 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;
+      
+      % Invoke  the appropriate causality procedure
+      cause_name = [comp_type, '_cause'];
+      if exist(cause_name)~=2 % Try a compound component -- need doing
+       mtt_info(sprintf('Component %s is unknown', comp_type), infofile);
+% $$$ 	sys_name = [comp_type, '_abg'];
+% $$$ 	agb2cbg(sys_name,comp_bonds,infofile);
+      else % its a simple component
+	cause_name
+	eval([ '[comp_bonds,status(i)] = ', cause_name, '(comp_bonds);' ]);
+        % Convert from component orientated to arrow orientated causality
+	bonds(bond_list,:) = comp_bonds.*direction; 
+      end;
+    end;
+
+    done = sum(sum(abs(bonds)))/total*100;
+%  mtt_info(sprintf('Causality is %3.0f%s complete.', done, pc), infofile);
+
+  end;
+end;
+
+final_done =  (sum(status==zeros(n_components,1))/n_components)*100;;
+
+mtt_info(sprintf('Final causality is %3.0f%s complete.', final_done, pc),
+infofile);
+
+% List overcausal bonds
+[over_causal_bonds,n] = getindex(status,1)
+if n>0
+  for i=over_causal_bonds'
+    eval([ '[comp_type,name] = ', system_name, '_cmp(i);' ]);
+    mtt_info(sprintf('Component %s (%s) is overcausal', name, comp_type), ...
+      infofile);
+  end;
+end;
+
+% List undercausal bonds
+[under_causal_bonds,n] = getindex(status,-1)
+if n>0
+  for i=under_causal_bonds'
+    eval([ '[comp_type,name] = ', system_name, '_cmp(i);' ]);
+    mtt_info(sprintf('Component %s (%s) is undercausal', name, comp_type), ...
+      infofile);
+  end;
+end;
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+