Trace: compute <X1,Z1> + ... <Xp,Zp> % % SDPNAL: % Copyright (c) 2008 by % Xinyuan Zhao, Defeng Sun, and Kim-Chuan Toh %**********************************************************************
0001 % Trace: compute <X1,Z1> + ... <Xp,Zp> 0002 %% 0003 %% SDPNAL: 0004 %% Copyright (c) 2008 by 0005 %% Xinyuan Zhao, Defeng Sun, and Kim-Chuan Toh 0006 %%********************************************************************** 0007 0008 function trXZ = blktrace(blk,X,Z) 0009 0010 trXZ = 0; 0011 for p = 1:size(blk,1) 0012 pblk = blk(p,:); 0013 if strcmp(pblk{1},'s') 0014 trXZ = trXZ + sum(sum(X{p}.*Z{p})); 0015 else 0016 trXZ = trXZ + X{p}'*Z{p}; 0017 end 0018 end 0019 %%********************************************************************** 0020 0021 0022 0023