print_stiffness_load.m

matlab实例

%***********************************************************************
%  function print_stiff_load
%
%  print out non-zero stiffness matrix and load vector members
%***********************************************************************
function print_stiffness_load(kglobal,load_vector,ndof)

%fprintf ('\n<<<<---- Internal Results ---->>>>\n\n')
if (ndof <= 8)   % print in matrix format
  fprintf ('STIFFNESS MATRIX\nr/c%9i%9i%9i%9i%9i%9i%9i%9i\n',...
           1,2,3,4,5,6,7,8)
  for i = 1:ndof
    fprintf ('%3i',i)
    for j = 1:ndof
      fprintf('%9.2g',kglobal(i,j))
    end
    fprintf('\n')
  end

  fprintf ('\nLOAD VECTOR\nRow     Value\n')
  for i = 1:ndof
    fprintf ('%3i %14.5g\n',i,load_vector(i))
  end
else               
  fprintf ('\n Stiffness and force terms only printed for small problems\n')
%   fprintf ('NON-ZERO STIFFNESS MATRIX ELEMENTS\n  #   row  col    value\n')
%   count = 0;
%   for i = 1:ndof
%     for j = 1:ndof
%       if (kglobal(i,j) ~= 0.0) % only print non-zero terms
%         count = count + 1;
%         fprintf ('%4i%4i%4i%14.5g\n',count,i,j,kglobal(i,j))
%       end
%     end
%   end
% 
%   fprintf ('NON-ZERO LOAD VECTOR ELEMENTS\n  #   row    value\n')
% 
%   count = 0;
%   for i = 1:ndof
%     if (load_vector(i) ~= 0.0) % only print non-zero terms
%       count = count + 1;
%       fprintf ('%3i%5i%14.5g\n',count,i,load_vector(i))
%     end
%   end
end