Nathaniel: /* A complicated operator on 4 local dimensions */

2018-06-26T13:43:01Z

A complicated operator on 4 local dimensions

Rmmateri at 16:35, 13 June 2018

2018-06-13T16:35:45Z

Nathaniel: Created page with "{{Function |name=IsBlockPositive |desc=Determines whether or not an operator is block positive |req=[http://cvxr.com/cvx/ cvx] |rel=IsPPT
IsPSD
IsSeparable..."

2014-09-23T18:51:51Z

Created page with "{{Function |name=IsBlockPositive |desc=Determines whether or not an operator is block positive |req=[http://cvxr.com/cvx/ cvx] |rel=IsPPT<br />IsPSD<br />IsSeparable..."

New page

{{Function
|name=IsBlockPositive
|desc=Determines whether or not an operator is block positive
|req=[http://cvxr.com/cvx/ cvx]
|rel=[[IsPPT]]<br />[[IsPSD]]<br />[[IsSeparable]]
|cat=[[List of functions#Entanglement_and_separability|Entanglement and separability]]
|upd=September 23, 2014
|v=0.50}}
<tt>'''IsBlockPositive'''</tt> is a [[List of functions|function]] that determines whether or not a bipartite operator is block positive (roughly speaking, it determines whether the operator is an entanglement witness or not). A value of 0 indicates that the operator is not block positive (and hence is not an entanglement witness), a value of 1 indicates that the operator is block positive (and hence is either an entanglement witness or is positive semidefinite), and a value of -1 indicates that the script was unable to determine whether or not the operator is block positive.

==Syntax==
* <tt>IBP = IsBlockPositive(X)</tt>
* <tt>IBP = IsBlockPositive(X,K)</tt>
* <tt>IBP = IsBlockPositive(X,K,DIM)</tt>
* <tt>IBP = IsBlockPositive(X,K,DIM,STR)</tt>
* <tt>IBP = IsBlockPositive(X,K,DIM,STR,TOL)</tt>
* <tt>[IBP,WIT] = IsBlockPositive(X,K,DIM,STR,TOL)</tt>

==Argument descriptions==
* <tt>X</tt>: A bipartite Hermitian operator.
* <tt>K</tt> (optional, default 1): A positive integer indicating that the script should determine whether or not <tt>X</tt> is <tt>K</tt>-block positive (i.e., whether or not it remains nonnegative under left and right multiplication by vectors with [[SchmidtRank|Schmidt rank]] <tt><= K</tt>).
* <tt>DIM</tt> (optional, by default has both subsystems of equal dimension): A 1-by-2 vector containing the dimensions of the two subsystems that <tt>X</tt> acts on.
* <tt>STR</tt> (optional, default 2): An integer that determines how hard the script should work to determine block positivity before giving up (<tt>STR = -1</tt> means that the script won't stop working until it finds an answer). Other valid values are <tt>0, 1, 2, 3, ...</tt>. In practice, if <tt>STR >= 4</tt> then most computers will run out of memory and/or the sun will explode before computation completes.
* <tt>TOL</tt> (optional, default <tt>eps^(3/8)</tt>): The numerical tolerance used throughout the script.

==Examples==
===The [[SwapOperator|swap operator]] is block positive===
The [[SwapOperator|swap operator]] always block positive, since it is the Choi matrix of the transpose map (and it is easy to see that the transpose map is positive). We can verify this in small dimensions as follows:
<syntaxhighlight>
>> for d = 2:10
res(d) = IsBlockPositive(SwapOperator(d));
end
>> res(2:10)

ans =

1 1 1 1 1 1 1 1 1
</syntaxhighlight>

However, the swap operator is ''not'' 2-block positive, which we can verify as follows:
<syntaxhighlight>
>> IsBlockPositive(SwapOperator(5),2)

ans =

0
</syntaxhighlight>

===The [[ChoiMap|Choi map]] is positive===
The [[ChoiMap|Choi map]] is a well-known example of a positive linear map that acts on 3-by-3 matrices. The following code verifies that this map is indeed positive, but not 2-positive:
<syntaxhighlight>
>> IsBlockPositive(ChoiMap())

ans =

1

>> IsBlockPositive(ChoiMap(),2)

ans =

0
</syntaxhighlight>

===A complicated operator on 4 local dimensions===
A somewhat messy positive linear map was introduced in Lemma 3 of <ref name="BCJRWY14">S. Bandyopadhyay, A. Cosentino, N. Johnston, V. Russo, J. Watrous, and N. Yu. ''Limitations on separable measurements by convex optimization''. E-print: [http://arxiv.org/abs/1408.6981 arXiv:1408.6981] [quant-ph], 2014.</ref>. The following code verifies that this map is indeed positive in the $t = 1$ case.
<syntaxhighlight>
>> b0 = Bell(0);
>> b3 = Bell(3);
>> v0 = kron(b0,b0);
>> Y = (kron(eye(4),b0*b0')/2 + kron(b3*b3',PartialTranspose(b3*b3',1)))/3 - v0*v0'/4;
>> Y = Swap(Y,[2,2,2,2],[2,3]); % this is the Choi matrix of the linear map given in Lemma 3
>> IsBlockPositive(Y)

ans =

1
</syntaxhighlight>

{{SourceCode|name=IsBlockPositive}}

==References==
<references />

@@ Line 70: / Line 70: @@
 >> v0 = kron(b0,b0);
 >> Y = (kron(eye(4),b0*b0')/2 + kron(b3*b3',PartialTranspose(b3*b3',1)))/3 - v0*v0'/4;
->> Y = Swap(Y,[2,2,2,2],[2,3]); % this is the Choi matrix of the linear map given in Lemma 3
+>> Y = Swap(Y,[2,3],[2,2,2,2]); % this is the Choi matrix of the linear map given in Lemma 3
 >> IsBlockPositive(Y)

@@ Line 2: / Line 2: @@
 |name=IsBlockPositive
 |desc=Determines whether or not an operator is block positive
-|req=[http://cvxr.com/cvx/ cvx]
+|req=[http://cvxr.com/cvx/ CVX]
 |rel=[[IsPPT]]<br />[[IsPSD]]<br />[[IsSeparable]]
 |cat=[[List of functions#Entanglement_and_separability|Entanglement&nbsp;and&nbsp;separability]]

IsBlockPositive - Revision history

Nathaniel: /* A complicated operator on 4 local dimensions */

Rmmateri at 16:35, 13 June 2018

Nathaniel: Created page with "{{Function |name=IsBlockPositive |desc=Determines whether or not an operator is block positive |req=[http://cvxr.com/cvx/ cvx] |rel=IsPPTIsPSDIsSeparable..."

Nathaniel: Created page with "{{Function |name=IsBlockPositive |desc=Determines whether or not an operator is block positive |req=[http://cvxr.com/cvx/ cvx] |rel=IsPPT
IsPSD
IsSeparable..."