Non-determinantal positive-definiteness conditions for $4 times 4$ Hermitian matrices with certain null...












2












$begingroup$


I have classes of $4 times 4$ Hermitian matrices $D$, the two $2 times 2$ diagonal blocks of which are themselves diagonal. That is, the (1,2),(2,1),(3,4) and (4,3) entries of $D$ are zero. (The other off-diagonal entries of $D$ can be real, complex or quaternionic.) I want necessary and sufficient (maybe minimal) conditions that $D$ be positive definite, that do not involve the full determinant of $D$. (In the quaternionic case, this would be the "Moore determinant", but the real and complex cases are the ones of immediate interest to me.)



As some background, this pertains to a quantum-information-theoretic problem, in which I am also interested in the "partial transpose" $D^{PT}$ of $D$, obtained, in general, by transposing in place the four $2 times 2$ blocks of $D$. The determinants of both $D$ and $D^{PT}$ are quite cumbersome, while the difference of the two determinants simplifies greatly, which is the basic motivation for my question. (I want to enforce the positive-definite nature of $D^{PT}$, subject to $D$ itself being positive-definite. Hopefully, the $2 times 2$ and $3 times 3$ minors of these matrices are more computationally amenable than their determinants. This pertains to the problem of determining the probability that two quantum bits ["qubits"] are disentangled/separable.)










share|cite|improve this question











$endgroup$












  • $begingroup$
    One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
    $endgroup$
    – Hanno
    Jan 15 at 22:34
















2












$begingroup$


I have classes of $4 times 4$ Hermitian matrices $D$, the two $2 times 2$ diagonal blocks of which are themselves diagonal. That is, the (1,2),(2,1),(3,4) and (4,3) entries of $D$ are zero. (The other off-diagonal entries of $D$ can be real, complex or quaternionic.) I want necessary and sufficient (maybe minimal) conditions that $D$ be positive definite, that do not involve the full determinant of $D$. (In the quaternionic case, this would be the "Moore determinant", but the real and complex cases are the ones of immediate interest to me.)



As some background, this pertains to a quantum-information-theoretic problem, in which I am also interested in the "partial transpose" $D^{PT}$ of $D$, obtained, in general, by transposing in place the four $2 times 2$ blocks of $D$. The determinants of both $D$ and $D^{PT}$ are quite cumbersome, while the difference of the two determinants simplifies greatly, which is the basic motivation for my question. (I want to enforce the positive-definite nature of $D^{PT}$, subject to $D$ itself being positive-definite. Hopefully, the $2 times 2$ and $3 times 3$ minors of these matrices are more computationally amenable than their determinants. This pertains to the problem of determining the probability that two quantum bits ["qubits"] are disentangled/separable.)










share|cite|improve this question











$endgroup$












  • $begingroup$
    One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
    $endgroup$
    – Hanno
    Jan 15 at 22:34














2












2








2





$begingroup$


I have classes of $4 times 4$ Hermitian matrices $D$, the two $2 times 2$ diagonal blocks of which are themselves diagonal. That is, the (1,2),(2,1),(3,4) and (4,3) entries of $D$ are zero. (The other off-diagonal entries of $D$ can be real, complex or quaternionic.) I want necessary and sufficient (maybe minimal) conditions that $D$ be positive definite, that do not involve the full determinant of $D$. (In the quaternionic case, this would be the "Moore determinant", but the real and complex cases are the ones of immediate interest to me.)



As some background, this pertains to a quantum-information-theoretic problem, in which I am also interested in the "partial transpose" $D^{PT}$ of $D$, obtained, in general, by transposing in place the four $2 times 2$ blocks of $D$. The determinants of both $D$ and $D^{PT}$ are quite cumbersome, while the difference of the two determinants simplifies greatly, which is the basic motivation for my question. (I want to enforce the positive-definite nature of $D^{PT}$, subject to $D$ itself being positive-definite. Hopefully, the $2 times 2$ and $3 times 3$ minors of these matrices are more computationally amenable than their determinants. This pertains to the problem of determining the probability that two quantum bits ["qubits"] are disentangled/separable.)










share|cite|improve this question











$endgroup$




I have classes of $4 times 4$ Hermitian matrices $D$, the two $2 times 2$ diagonal blocks of which are themselves diagonal. That is, the (1,2),(2,1),(3,4) and (4,3) entries of $D$ are zero. (The other off-diagonal entries of $D$ can be real, complex or quaternionic.) I want necessary and sufficient (maybe minimal) conditions that $D$ be positive definite, that do not involve the full determinant of $D$. (In the quaternionic case, this would be the "Moore determinant", but the real and complex cases are the ones of immediate interest to me.)



As some background, this pertains to a quantum-information-theoretic problem, in which I am also interested in the "partial transpose" $D^{PT}$ of $D$, obtained, in general, by transposing in place the four $2 times 2$ blocks of $D$. The determinants of both $D$ and $D^{PT}$ are quite cumbersome, while the difference of the two determinants simplifies greatly, which is the basic motivation for my question. (I want to enforce the positive-definite nature of $D^{PT}$, subject to $D$ itself being positive-definite. Hopefully, the $2 times 2$ and $3 times 3$ minors of these matrices are more computationally amenable than their determinants. This pertains to the problem of determining the probability that two quantum bits ["qubits"] are disentangled/separable.)







linear-algebra matrices complex-numbers determinant quaternions






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited May 27 '17 at 16:00







Paul B. Slater

















asked May 27 '17 at 15:53









Paul B. SlaterPaul B. Slater

22319




22319












  • $begingroup$
    One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
    $endgroup$
    – Hanno
    Jan 15 at 22:34


















  • $begingroup$
    One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
    $endgroup$
    – Hanno
    Jan 15 at 22:34
















$begingroup$
One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
$endgroup$
– Hanno
Jan 15 at 22:34




$begingroup$
One way to characterise positive matrices is by the condition that all eigenvalues are positive. So I do not understand your phrase "... conditions that $D$ be positive definite, that do not involve the full determinant of $D$.", because the determinant alone, which equals the product of all eigenvalues, is not sufficient to detect or disproof positivity.
$endgroup$
– Hanno
Jan 15 at 22:34










1 Answer
1






active

oldest

votes


















0












$begingroup$

Let's go and derive first an equivalent condition on a related $2times 2$ matrix which guarantees positive-definiteness of $,D$, valid in either the complex or the real case.



(1) Fixing notation:
$$D:=;begin{pmatrix} D_1& M\ M^*& D_2end{pmatrix}
:=;begin{pmatrix} begin{pmatrix}d_1& 0\ 0& d_2end{pmatrix}& begin{pmatrix}m_{13}& m_{14}\ m_{23}& m_{24}end{pmatrix}\[1ex]
M^*& begin{pmatrix}d_3& 0\ 0& d_4end{pmatrix}end{pmatrix}$$

There are no a priori restrictions on the entries of $M$, and their indices stick to the indexation as of $D$.



It is henceforth assumed that $,d_1,d_2,d_3,d_4>0$.

This feature is a necessary condition to $,D,$ being positive-definite, because
$,d_i =langle e_i| De_irangle>0$, with $e_i$ being a standard unit vector.
Using it right away one has the equivalence
$$0<D;iff;
0,<:begin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
Dbegin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
;=; begin{pmatrix}mathbb{1}& D_1^{-1/2}MD_2^{-1/2}\
D_2^{-1/2}M^*D_1^{-1/2}& mathbb{1}end{pmatrix}$$

since positive-definiteness is preserved under congruences
$Amapsto T^*AT$ with invertible $T$. Now use that
$$0 < begin{pmatrix}mathbb{1}& C\ C^*& mathbb{1}end{pmatrix}
;iff;|C|<1tag{1}$$

where $|cdot|$ denotes the operator norm. It equals the largest singular value.
The little sister of $(1)$ within $M_2(mathbb C)$ is the equivalence
$;0<left(begin{smallmatrix}1& z\ overline z& 1end{smallmatrix}right)$
iff $|z|<1$ (and she usually shows up in the proof(*) of the general case).




Thus the equivalence criterion for positive-definiteness of $,D,$ reads
$,left|D_1^{-1/2}MD_2^{-1/2}right|,stackrel{!}{<}, 1$

put another way: $,D_1^{-1/2}MD_2^{-1/2} = left(frac{m_{ij}}{sqrt{d_id_j}}right)_{i=1,2;j=3,4}$ is a strict contraction.




(2) An alternative, more hands-on approach would be based on the
Cholesky decomposition. It is the "$Rightarrow$" in
$$A>0:iff,exists! text{ lower-triangular matrix }Ltext{ with },LL^*=A,,$$



and many (to most) implementations of the Cholesky algorithm accept as input a matrix without definiteness properties, and return $L$, or some educated error message if the input was not positive (enough). Cholesky still works for positive-semidefinite input and yields a lower-triangular $L$ but its uniqueness is lost.



As $n=4$ this "positivity criterion" is checked within a blink of the eye $:ddotsmile$



* e.g. Proposition I.3.5 in R. Bhatia's book "Matrix Analysis"
(Springer, GTM 169, 1997)






share|cite|improve this answer









$endgroup$













    Your Answer





    StackExchange.ifUsing("editor", function () {
    return StackExchange.using("mathjaxEditing", function () {
    StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
    StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
    });
    });
    }, "mathjax-editing");

    StackExchange.ready(function() {
    var channelOptions = {
    tags: "".split(" "),
    id: "69"
    };
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function() {
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled) {
    StackExchange.using("snippets", function() {
    createEditor();
    });
    }
    else {
    createEditor();
    }
    });

    function createEditor() {
    StackExchange.prepareEditor({
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: true,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: 10,
    bindNavPrevention: true,
    postfix: "",
    imageUploader: {
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    },
    noCode: true, onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    });


    }
    });














    draft saved

    draft discarded


















    StackExchange.ready(
    function () {
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2299000%2fnon-determinantal-positive-definiteness-conditions-for-4-times-4-hermitian-ma%23new-answer', 'question_page');
    }
    );

    Post as a guest















    Required, but never shown

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    0












    $begingroup$

    Let's go and derive first an equivalent condition on a related $2times 2$ matrix which guarantees positive-definiteness of $,D$, valid in either the complex or the real case.



    (1) Fixing notation:
    $$D:=;begin{pmatrix} D_1& M\ M^*& D_2end{pmatrix}
    :=;begin{pmatrix} begin{pmatrix}d_1& 0\ 0& d_2end{pmatrix}& begin{pmatrix}m_{13}& m_{14}\ m_{23}& m_{24}end{pmatrix}\[1ex]
    M^*& begin{pmatrix}d_3& 0\ 0& d_4end{pmatrix}end{pmatrix}$$

    There are no a priori restrictions on the entries of $M$, and their indices stick to the indexation as of $D$.



    It is henceforth assumed that $,d_1,d_2,d_3,d_4>0$.

    This feature is a necessary condition to $,D,$ being positive-definite, because
    $,d_i =langle e_i| De_irangle>0$, with $e_i$ being a standard unit vector.
    Using it right away one has the equivalence
    $$0<D;iff;
    0,<:begin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
    Dbegin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
    ;=; begin{pmatrix}mathbb{1}& D_1^{-1/2}MD_2^{-1/2}\
    D_2^{-1/2}M^*D_1^{-1/2}& mathbb{1}end{pmatrix}$$

    since positive-definiteness is preserved under congruences
    $Amapsto T^*AT$ with invertible $T$. Now use that
    $$0 < begin{pmatrix}mathbb{1}& C\ C^*& mathbb{1}end{pmatrix}
    ;iff;|C|<1tag{1}$$

    where $|cdot|$ denotes the operator norm. It equals the largest singular value.
    The little sister of $(1)$ within $M_2(mathbb C)$ is the equivalence
    $;0<left(begin{smallmatrix}1& z\ overline z& 1end{smallmatrix}right)$
    iff $|z|<1$ (and she usually shows up in the proof(*) of the general case).




    Thus the equivalence criterion for positive-definiteness of $,D,$ reads
    $,left|D_1^{-1/2}MD_2^{-1/2}right|,stackrel{!}{<}, 1$

    put another way: $,D_1^{-1/2}MD_2^{-1/2} = left(frac{m_{ij}}{sqrt{d_id_j}}right)_{i=1,2;j=3,4}$ is a strict contraction.




    (2) An alternative, more hands-on approach would be based on the
    Cholesky decomposition. It is the "$Rightarrow$" in
    $$A>0:iff,exists! text{ lower-triangular matrix }Ltext{ with },LL^*=A,,$$



    and many (to most) implementations of the Cholesky algorithm accept as input a matrix without definiteness properties, and return $L$, or some educated error message if the input was not positive (enough). Cholesky still works for positive-semidefinite input and yields a lower-triangular $L$ but its uniqueness is lost.



    As $n=4$ this "positivity criterion" is checked within a blink of the eye $:ddotsmile$



    * e.g. Proposition I.3.5 in R. Bhatia's book "Matrix Analysis"
    (Springer, GTM 169, 1997)






    share|cite|improve this answer









    $endgroup$


















      0












      $begingroup$

      Let's go and derive first an equivalent condition on a related $2times 2$ matrix which guarantees positive-definiteness of $,D$, valid in either the complex or the real case.



      (1) Fixing notation:
      $$D:=;begin{pmatrix} D_1& M\ M^*& D_2end{pmatrix}
      :=;begin{pmatrix} begin{pmatrix}d_1& 0\ 0& d_2end{pmatrix}& begin{pmatrix}m_{13}& m_{14}\ m_{23}& m_{24}end{pmatrix}\[1ex]
      M^*& begin{pmatrix}d_3& 0\ 0& d_4end{pmatrix}end{pmatrix}$$

      There are no a priori restrictions on the entries of $M$, and their indices stick to the indexation as of $D$.



      It is henceforth assumed that $,d_1,d_2,d_3,d_4>0$.

      This feature is a necessary condition to $,D,$ being positive-definite, because
      $,d_i =langle e_i| De_irangle>0$, with $e_i$ being a standard unit vector.
      Using it right away one has the equivalence
      $$0<D;iff;
      0,<:begin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
      Dbegin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
      ;=; begin{pmatrix}mathbb{1}& D_1^{-1/2}MD_2^{-1/2}\
      D_2^{-1/2}M^*D_1^{-1/2}& mathbb{1}end{pmatrix}$$

      since positive-definiteness is preserved under congruences
      $Amapsto T^*AT$ with invertible $T$. Now use that
      $$0 < begin{pmatrix}mathbb{1}& C\ C^*& mathbb{1}end{pmatrix}
      ;iff;|C|<1tag{1}$$

      where $|cdot|$ denotes the operator norm. It equals the largest singular value.
      The little sister of $(1)$ within $M_2(mathbb C)$ is the equivalence
      $;0<left(begin{smallmatrix}1& z\ overline z& 1end{smallmatrix}right)$
      iff $|z|<1$ (and she usually shows up in the proof(*) of the general case).




      Thus the equivalence criterion for positive-definiteness of $,D,$ reads
      $,left|D_1^{-1/2}MD_2^{-1/2}right|,stackrel{!}{<}, 1$

      put another way: $,D_1^{-1/2}MD_2^{-1/2} = left(frac{m_{ij}}{sqrt{d_id_j}}right)_{i=1,2;j=3,4}$ is a strict contraction.




      (2) An alternative, more hands-on approach would be based on the
      Cholesky decomposition. It is the "$Rightarrow$" in
      $$A>0:iff,exists! text{ lower-triangular matrix }Ltext{ with },LL^*=A,,$$



      and many (to most) implementations of the Cholesky algorithm accept as input a matrix without definiteness properties, and return $L$, or some educated error message if the input was not positive (enough). Cholesky still works for positive-semidefinite input and yields a lower-triangular $L$ but its uniqueness is lost.



      As $n=4$ this "positivity criterion" is checked within a blink of the eye $:ddotsmile$



      * e.g. Proposition I.3.5 in R. Bhatia's book "Matrix Analysis"
      (Springer, GTM 169, 1997)






      share|cite|improve this answer









      $endgroup$
















        0












        0








        0





        $begingroup$

        Let's go and derive first an equivalent condition on a related $2times 2$ matrix which guarantees positive-definiteness of $,D$, valid in either the complex or the real case.



        (1) Fixing notation:
        $$D:=;begin{pmatrix} D_1& M\ M^*& D_2end{pmatrix}
        :=;begin{pmatrix} begin{pmatrix}d_1& 0\ 0& d_2end{pmatrix}& begin{pmatrix}m_{13}& m_{14}\ m_{23}& m_{24}end{pmatrix}\[1ex]
        M^*& begin{pmatrix}d_3& 0\ 0& d_4end{pmatrix}end{pmatrix}$$

        There are no a priori restrictions on the entries of $M$, and their indices stick to the indexation as of $D$.



        It is henceforth assumed that $,d_1,d_2,d_3,d_4>0$.

        This feature is a necessary condition to $,D,$ being positive-definite, because
        $,d_i =langle e_i| De_irangle>0$, with $e_i$ being a standard unit vector.
        Using it right away one has the equivalence
        $$0<D;iff;
        0,<:begin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
        Dbegin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
        ;=; begin{pmatrix}mathbb{1}& D_1^{-1/2}MD_2^{-1/2}\
        D_2^{-1/2}M^*D_1^{-1/2}& mathbb{1}end{pmatrix}$$

        since positive-definiteness is preserved under congruences
        $Amapsto T^*AT$ with invertible $T$. Now use that
        $$0 < begin{pmatrix}mathbb{1}& C\ C^*& mathbb{1}end{pmatrix}
        ;iff;|C|<1tag{1}$$

        where $|cdot|$ denotes the operator norm. It equals the largest singular value.
        The little sister of $(1)$ within $M_2(mathbb C)$ is the equivalence
        $;0<left(begin{smallmatrix}1& z\ overline z& 1end{smallmatrix}right)$
        iff $|z|<1$ (and she usually shows up in the proof(*) of the general case).




        Thus the equivalence criterion for positive-definiteness of $,D,$ reads
        $,left|D_1^{-1/2}MD_2^{-1/2}right|,stackrel{!}{<}, 1$

        put another way: $,D_1^{-1/2}MD_2^{-1/2} = left(frac{m_{ij}}{sqrt{d_id_j}}right)_{i=1,2;j=3,4}$ is a strict contraction.




        (2) An alternative, more hands-on approach would be based on the
        Cholesky decomposition. It is the "$Rightarrow$" in
        $$A>0:iff,exists! text{ lower-triangular matrix }Ltext{ with },LL^*=A,,$$



        and many (to most) implementations of the Cholesky algorithm accept as input a matrix without definiteness properties, and return $L$, or some educated error message if the input was not positive (enough). Cholesky still works for positive-semidefinite input and yields a lower-triangular $L$ but its uniqueness is lost.



        As $n=4$ this "positivity criterion" is checked within a blink of the eye $:ddotsmile$



        * e.g. Proposition I.3.5 in R. Bhatia's book "Matrix Analysis"
        (Springer, GTM 169, 1997)






        share|cite|improve this answer









        $endgroup$



        Let's go and derive first an equivalent condition on a related $2times 2$ matrix which guarantees positive-definiteness of $,D$, valid in either the complex or the real case.



        (1) Fixing notation:
        $$D:=;begin{pmatrix} D_1& M\ M^*& D_2end{pmatrix}
        :=;begin{pmatrix} begin{pmatrix}d_1& 0\ 0& d_2end{pmatrix}& begin{pmatrix}m_{13}& m_{14}\ m_{23}& m_{24}end{pmatrix}\[1ex]
        M^*& begin{pmatrix}d_3& 0\ 0& d_4end{pmatrix}end{pmatrix}$$

        There are no a priori restrictions on the entries of $M$, and their indices stick to the indexation as of $D$.



        It is henceforth assumed that $,d_1,d_2,d_3,d_4>0$.

        This feature is a necessary condition to $,D,$ being positive-definite, because
        $,d_i =langle e_i| De_irangle>0$, with $e_i$ being a standard unit vector.
        Using it right away one has the equivalence
        $$0<D;iff;
        0,<:begin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
        Dbegin{pmatrix} D_1^{-1/2}& mathbb{0}\ mathbb{0}& D_2^{-1/2}end{pmatrix}
        ;=; begin{pmatrix}mathbb{1}& D_1^{-1/2}MD_2^{-1/2}\
        D_2^{-1/2}M^*D_1^{-1/2}& mathbb{1}end{pmatrix}$$

        since positive-definiteness is preserved under congruences
        $Amapsto T^*AT$ with invertible $T$. Now use that
        $$0 < begin{pmatrix}mathbb{1}& C\ C^*& mathbb{1}end{pmatrix}
        ;iff;|C|<1tag{1}$$

        where $|cdot|$ denotes the operator norm. It equals the largest singular value.
        The little sister of $(1)$ within $M_2(mathbb C)$ is the equivalence
        $;0<left(begin{smallmatrix}1& z\ overline z& 1end{smallmatrix}right)$
        iff $|z|<1$ (and she usually shows up in the proof(*) of the general case).




        Thus the equivalence criterion for positive-definiteness of $,D,$ reads
        $,left|D_1^{-1/2}MD_2^{-1/2}right|,stackrel{!}{<}, 1$

        put another way: $,D_1^{-1/2}MD_2^{-1/2} = left(frac{m_{ij}}{sqrt{d_id_j}}right)_{i=1,2;j=3,4}$ is a strict contraction.




        (2) An alternative, more hands-on approach would be based on the
        Cholesky decomposition. It is the "$Rightarrow$" in
        $$A>0:iff,exists! text{ lower-triangular matrix }Ltext{ with },LL^*=A,,$$



        and many (to most) implementations of the Cholesky algorithm accept as input a matrix without definiteness properties, and return $L$, or some educated error message if the input was not positive (enough). Cholesky still works for positive-semidefinite input and yields a lower-triangular $L$ but its uniqueness is lost.



        As $n=4$ this "positivity criterion" is checked within a blink of the eye $:ddotsmile$



        * e.g. Proposition I.3.5 in R. Bhatia's book "Matrix Analysis"
        (Springer, GTM 169, 1997)







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered Jan 16 at 15:45









        HannoHanno

        2,240528




        2,240528






























            draft saved

            draft discarded




















































            Thanks for contributing an answer to Mathematics Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid



            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.


            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2299000%2fnon-determinantal-positive-definiteness-conditions-for-4-times-4-hermitian-ma%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            Popular posts from this blog

            MongoDB - Not Authorized To Execute Command

            Npm cannot find a required file even through it is in the searched directory

            in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith