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;