Changeset 41 for distools/pe_mtimes.m

Timestamp:

02/17/13 14:56:55 (12 years ago)

Author:

bduin

Message:

 

File:

• distools/pe_mtimes.m (modified) (2 diffs)

distools/pe_mtimes.m

@@ -3 +3 @@
 %  C = PE_MTIMES(A,B)
 %
-% A and/or B should be datasets. Their PE signature SIG is retrieved,
+% A should be a PE dataset. Its PE signature SIG is retrieved,
 % C = A*J*B, where J is a diagonal matrix with 1's, followed by -1's.
 % J = diag ([ONES(SIG(1),1);  -ONES(SIG(2),1)])
@@ -9 +9 @@
 function c = pe_mtimes(a,b)
 
-siga = getsig(a);
-sigb = getsig(b); 
-
-if isdataset(a) & isdataset(b)
-  if any(siga~=sigb)
-    error('PE signatures should be equal')
-  end
-  J = diag([ones(1,siga(1))  -ones(1,siga(2))]);
-  c = a*J*b;
-elseif isdataset(a)
-  J = diag([ones(1,siga(1))  -ones(1,siga(2))]);
-  c = a*J*b;
-elseif isdataset(b) % problem! make doubles
-  J = diag([ones(1,sigb(1))  -ones(1,sigb(2))]);
-  c = a*J*(+b);
-else % just doubles, no PE space
-  c = a*b;
-end
+isdataset(a);
+sig = getsig(a); 
+J = diag([ones(1,sig(1))  -ones(1,sig(2))]);
+c = a*J*b;