Weak Jacobian of Proximal Operator












5














Given a convex function $g(x):mathbb{R}^nrightarrow mathbb{R}$, the proximal operator of $g$ is defined as



$P_g(x)=underset{u}{argmin}quad frac{1}{2}||x-u||_2^2+g(u)$.



Since $g(x)$ is convex, the proximal is a singleton, i.e., there is a unique minimizer $u$. Thus, the proximal operator is a vector function $P_g:mathbb{R}^nrightarrowmathbb{R}^n$.



I'm trying to find an expression for the Jacobian (or weak Jacobian) of $P_g$. In the case where $g$ is differentiable, I believe I can find the Jacobian using the implicit function theorem. However, I'm interested in the general case where $g$ is convex but not differentiable.



As first try, I tried to look on the Moreau Envelope of $g$:
$M_g(x)=underset{u}{min}quad frac{1}{2}||x-u||_2^2+g(u)$.



This function is differentiable and its gradient is given by



$nabla M_g(x) = x-P_g(x)$ .



Therefore, if I could compute the Hessian $H_g$ of $M_g$, I'll get



$H_g=I-J_g$



where $J_g$ is the desired Jacobian.



However, I'm not sure under which conditions the Hessian exists (even in the weak sense) and what the expression for it? Moreover, maybe there is another way to approach this.










share|cite|improve this question







New contributor




Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

























    5














    Given a convex function $g(x):mathbb{R}^nrightarrow mathbb{R}$, the proximal operator of $g$ is defined as



    $P_g(x)=underset{u}{argmin}quad frac{1}{2}||x-u||_2^2+g(u)$.



    Since $g(x)$ is convex, the proximal is a singleton, i.e., there is a unique minimizer $u$. Thus, the proximal operator is a vector function $P_g:mathbb{R}^nrightarrowmathbb{R}^n$.



    I'm trying to find an expression for the Jacobian (or weak Jacobian) of $P_g$. In the case where $g$ is differentiable, I believe I can find the Jacobian using the implicit function theorem. However, I'm interested in the general case where $g$ is convex but not differentiable.



    As first try, I tried to look on the Moreau Envelope of $g$:
    $M_g(x)=underset{u}{min}quad frac{1}{2}||x-u||_2^2+g(u)$.



    This function is differentiable and its gradient is given by



    $nabla M_g(x) = x-P_g(x)$ .



    Therefore, if I could compute the Hessian $H_g$ of $M_g$, I'll get



    $H_g=I-J_g$



    where $J_g$ is the desired Jacobian.



    However, I'm not sure under which conditions the Hessian exists (even in the weak sense) and what the expression for it? Moreover, maybe there is another way to approach this.










    share|cite|improve this question







    New contributor




    Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
    Check out our Code of Conduct.























      5












      5








      5


      4





      Given a convex function $g(x):mathbb{R}^nrightarrow mathbb{R}$, the proximal operator of $g$ is defined as



      $P_g(x)=underset{u}{argmin}quad frac{1}{2}||x-u||_2^2+g(u)$.



      Since $g(x)$ is convex, the proximal is a singleton, i.e., there is a unique minimizer $u$. Thus, the proximal operator is a vector function $P_g:mathbb{R}^nrightarrowmathbb{R}^n$.



      I'm trying to find an expression for the Jacobian (or weak Jacobian) of $P_g$. In the case where $g$ is differentiable, I believe I can find the Jacobian using the implicit function theorem. However, I'm interested in the general case where $g$ is convex but not differentiable.



      As first try, I tried to look on the Moreau Envelope of $g$:
      $M_g(x)=underset{u}{min}quad frac{1}{2}||x-u||_2^2+g(u)$.



      This function is differentiable and its gradient is given by



      $nabla M_g(x) = x-P_g(x)$ .



      Therefore, if I could compute the Hessian $H_g$ of $M_g$, I'll get



      $H_g=I-J_g$



      where $J_g$ is the desired Jacobian.



      However, I'm not sure under which conditions the Hessian exists (even in the weak sense) and what the expression for it? Moreover, maybe there is another way to approach this.










      share|cite|improve this question







      New contributor




      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.











      Given a convex function $g(x):mathbb{R}^nrightarrow mathbb{R}$, the proximal operator of $g$ is defined as



      $P_g(x)=underset{u}{argmin}quad frac{1}{2}||x-u||_2^2+g(u)$.



      Since $g(x)$ is convex, the proximal is a singleton, i.e., there is a unique minimizer $u$. Thus, the proximal operator is a vector function $P_g:mathbb{R}^nrightarrowmathbb{R}^n$.



      I'm trying to find an expression for the Jacobian (or weak Jacobian) of $P_g$. In the case where $g$ is differentiable, I believe I can find the Jacobian using the implicit function theorem. However, I'm interested in the general case where $g$ is convex but not differentiable.



      As first try, I tried to look on the Moreau Envelope of $g$:
      $M_g(x)=underset{u}{min}quad frac{1}{2}||x-u||_2^2+g(u)$.



      This function is differentiable and its gradient is given by



      $nabla M_g(x) = x-P_g(x)$ .



      Therefore, if I could compute the Hessian $H_g$ of $M_g$, I'll get



      $H_g=I-J_g$



      where $J_g$ is the desired Jacobian.



      However, I'm not sure under which conditions the Hessian exists (even in the weak sense) and what the expression for it? Moreover, maybe there is another way to approach this.







      optimization convex-analysis convex-optimization weak-derivatives






      share|cite|improve this question







      New contributor




      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.











      share|cite|improve this question







      New contributor




      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.









      share|cite|improve this question




      share|cite|improve this question






      New contributor




      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.









      asked Dec 26 at 15:54









      Regev Cohen

      262




      262




      New contributor




      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.





      New contributor





      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.






      Regev Cohen is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.



























          active

          oldest

          votes











          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
          });


          }
          });






          Regev Cohen is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3053058%2fweak-jacobian-of-proximal-operator%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown






























          active

          oldest

          votes













          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          Regev Cohen is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          Regev Cohen is a new contributor. Be nice, and check out our Code of Conduct.













          Regev Cohen is a new contributor. Be nice, and check out our Code of Conduct.












          Regev Cohen is a new contributor. Be nice, and check out our Code of Conduct.
















          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.





          Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


          Please pay close attention to the following guidance:


          • 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.


          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%2f3053058%2fweak-jacobian-of-proximal-operator%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

          Human spaceflight

          Can not write log (Is /dev/pts mounted?) - openpty in Ubuntu-on-Windows?

          File:DeusFollowingSea.jpg