Hessian-vector products












0












$begingroup$


Can someone explain why this is true?



$$g(x + Delta x) = g(x) + H(x) Delta x$$



where g is the gradient of function f(x) with respect to x, and H is the hessian of f(x) with respect to x.



I would really appreciate a detailed derivation because I don't understand what it means to take the gradient/derivative of $Delta x$ which represents the change in x.



Source towards the top here: https://justindomke.wordpress.com/2009/01/17/hessian-vector-products/










share|cite|improve this question











$endgroup$












  • $begingroup$
    Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
    $endgroup$
    – Emil
    Nov 28 '18 at 7:14












  • $begingroup$
    Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
    $endgroup$
    – Emil
    Nov 28 '18 at 7:22












  • $begingroup$
    $Delta x$ is a vector I believe
    $endgroup$
    – user3180
    Nov 28 '18 at 7:58












  • $begingroup$
    I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
    $endgroup$
    – user3180
    Nov 28 '18 at 8:01










  • $begingroup$
    user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
    $endgroup$
    – Emil
    Nov 28 '18 at 13:26


















0












$begingroup$


Can someone explain why this is true?



$$g(x + Delta x) = g(x) + H(x) Delta x$$



where g is the gradient of function f(x) with respect to x, and H is the hessian of f(x) with respect to x.



I would really appreciate a detailed derivation because I don't understand what it means to take the gradient/derivative of $Delta x$ which represents the change in x.



Source towards the top here: https://justindomke.wordpress.com/2009/01/17/hessian-vector-products/










share|cite|improve this question











$endgroup$












  • $begingroup$
    Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
    $endgroup$
    – Emil
    Nov 28 '18 at 7:14












  • $begingroup$
    Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
    $endgroup$
    – Emil
    Nov 28 '18 at 7:22












  • $begingroup$
    $Delta x$ is a vector I believe
    $endgroup$
    – user3180
    Nov 28 '18 at 7:58












  • $begingroup$
    I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
    $endgroup$
    – user3180
    Nov 28 '18 at 8:01










  • $begingroup$
    user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
    $endgroup$
    – Emil
    Nov 28 '18 at 13:26
















0












0








0





$begingroup$


Can someone explain why this is true?



$$g(x + Delta x) = g(x) + H(x) Delta x$$



where g is the gradient of function f(x) with respect to x, and H is the hessian of f(x) with respect to x.



I would really appreciate a detailed derivation because I don't understand what it means to take the gradient/derivative of $Delta x$ which represents the change in x.



Source towards the top here: https://justindomke.wordpress.com/2009/01/17/hessian-vector-products/










share|cite|improve this question











$endgroup$




Can someone explain why this is true?



$$g(x + Delta x) = g(x) + H(x) Delta x$$



where g is the gradient of function f(x) with respect to x, and H is the hessian of f(x) with respect to x.



I would really appreciate a detailed derivation because I don't understand what it means to take the gradient/derivative of $Delta x$ which represents the change in x.



Source towards the top here: https://justindomke.wordpress.com/2009/01/17/hessian-vector-products/







optimization machine-learning






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Nov 28 '18 at 20:48







user3180

















asked Nov 28 '18 at 6:41









user3180user3180

1185




1185












  • $begingroup$
    Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
    $endgroup$
    – Emil
    Nov 28 '18 at 7:14












  • $begingroup$
    Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
    $endgroup$
    – Emil
    Nov 28 '18 at 7:22












  • $begingroup$
    $Delta x$ is a vector I believe
    $endgroup$
    – user3180
    Nov 28 '18 at 7:58












  • $begingroup$
    I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
    $endgroup$
    – user3180
    Nov 28 '18 at 8:01










  • $begingroup$
    user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
    $endgroup$
    – Emil
    Nov 28 '18 at 13:26




















  • $begingroup$
    Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
    $endgroup$
    – Emil
    Nov 28 '18 at 7:14












  • $begingroup$
    Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
    $endgroup$
    – Emil
    Nov 28 '18 at 7:22












  • $begingroup$
    $Delta x$ is a vector I believe
    $endgroup$
    – user3180
    Nov 28 '18 at 7:58












  • $begingroup$
    I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
    $endgroup$
    – user3180
    Nov 28 '18 at 8:01










  • $begingroup$
    user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
    $endgroup$
    – Emil
    Nov 28 '18 at 13:26


















$begingroup$
Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
$endgroup$
– Emil
Nov 28 '18 at 7:14






$begingroup$
Is it like the third term in the taylor expansion? en.m.wikipedia.org/wiki/…
$endgroup$
– Emil
Nov 28 '18 at 7:14














$begingroup$
Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
$endgroup$
– Emil
Nov 28 '18 at 7:22






$begingroup$
Is the lhs supposed to be the gradient taken at a point or is it a scalar? The hessian is a second order tensor, so it should eat two dx to become a scalar, not one, I believe.
$endgroup$
– Emil
Nov 28 '18 at 7:22














$begingroup$
$Delta x$ is a vector I believe
$endgroup$
– user3180
Nov 28 '18 at 7:58






$begingroup$
$Delta x$ is a vector I believe
$endgroup$
– user3180
Nov 28 '18 at 7:58














$begingroup$
I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
$endgroup$
– user3180
Nov 28 '18 at 8:01




$begingroup$
I don't know why you're talking about scalars, the Hessian H(x) is a vector and so is $Delta x$. Is it because the source link is $g(x) = frac{d}{dx} f(x)$ is using scalar notation? He intends this to be the gradient of a vector
$endgroup$
– user3180
Nov 28 '18 at 8:01












$begingroup$
user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
$endgroup$
– Emil
Nov 28 '18 at 13:26






$begingroup$
user3180: to know if there is an error and if so if it is on lhs or rhs (hessian is a second order tensor not a vector, according to some definitions, or maybe even third order tensor)
$endgroup$
– Emil
Nov 28 '18 at 13:26












1 Answer
1






active

oldest

votes


















0












$begingroup$

For a scalar function f,



$f(x+dx) approx f(x)+ frac{d f(x)}{dx}* dx$



Vector case [$x$ is (n,1) dimensional vector]:



$g(x + Delta x) approx g(x) + g(Delta x)$



$g(x) = nabla f(x) = begin{bmatrix}
g_1{(x)} \
g_2{(x)} \
vdots \
g_n{(x)}
end{bmatrix}$
$mid g_i{(x)} = frac{partial f(x)}{dx_i}$



$H(x) = nabla g(x) = begin{bmatrix}
nabla g_1{(x)}^T \
nabla g_2{(x)}^T \
vdots \
nabla g_n{(x)}^T
end{bmatrix}$



$g(x) +g(Delta x) = g(x) + Delta g(x)approx g(x)+ begin{bmatrix}
nabla g_1{(x)}^T Delta x \
nabla g_2{(x)}^T Delta x\
vdots \
nabla g_n{(x)}^T Delta x
end{bmatrix}$



Translation:
$g(Delta x)$ represents g with a small change in x as input, which is equivalent to $Delta g(x)$ [the infinitesimal change in g(x)].



If you look at the last column vector, you will see that each row $nabla g_i(x)^T Delta x$ represents the (rate of change of $g_i(x)$/rate of change of x) dotted with (infinitesimal change in x) := infinitesimal change in $g_i(x)$. In aggregate, the column vector represents the infinitesimal change in g(x) := $Delta g(x)$.



See also: http://www.friesian.com/calculus.htm






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%2f3016823%2fhessian-vector-products%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$

    For a scalar function f,



    $f(x+dx) approx f(x)+ frac{d f(x)}{dx}* dx$



    Vector case [$x$ is (n,1) dimensional vector]:



    $g(x + Delta x) approx g(x) + g(Delta x)$



    $g(x) = nabla f(x) = begin{bmatrix}
    g_1{(x)} \
    g_2{(x)} \
    vdots \
    g_n{(x)}
    end{bmatrix}$
    $mid g_i{(x)} = frac{partial f(x)}{dx_i}$



    $H(x) = nabla g(x) = begin{bmatrix}
    nabla g_1{(x)}^T \
    nabla g_2{(x)}^T \
    vdots \
    nabla g_n{(x)}^T
    end{bmatrix}$



    $g(x) +g(Delta x) = g(x) + Delta g(x)approx g(x)+ begin{bmatrix}
    nabla g_1{(x)}^T Delta x \
    nabla g_2{(x)}^T Delta x\
    vdots \
    nabla g_n{(x)}^T Delta x
    end{bmatrix}$



    Translation:
    $g(Delta x)$ represents g with a small change in x as input, which is equivalent to $Delta g(x)$ [the infinitesimal change in g(x)].



    If you look at the last column vector, you will see that each row $nabla g_i(x)^T Delta x$ represents the (rate of change of $g_i(x)$/rate of change of x) dotted with (infinitesimal change in x) := infinitesimal change in $g_i(x)$. In aggregate, the column vector represents the infinitesimal change in g(x) := $Delta g(x)$.



    See also: http://www.friesian.com/calculus.htm






    share|cite|improve this answer











    $endgroup$


















      0












      $begingroup$

      For a scalar function f,



      $f(x+dx) approx f(x)+ frac{d f(x)}{dx}* dx$



      Vector case [$x$ is (n,1) dimensional vector]:



      $g(x + Delta x) approx g(x) + g(Delta x)$



      $g(x) = nabla f(x) = begin{bmatrix}
      g_1{(x)} \
      g_2{(x)} \
      vdots \
      g_n{(x)}
      end{bmatrix}$
      $mid g_i{(x)} = frac{partial f(x)}{dx_i}$



      $H(x) = nabla g(x) = begin{bmatrix}
      nabla g_1{(x)}^T \
      nabla g_2{(x)}^T \
      vdots \
      nabla g_n{(x)}^T
      end{bmatrix}$



      $g(x) +g(Delta x) = g(x) + Delta g(x)approx g(x)+ begin{bmatrix}
      nabla g_1{(x)}^T Delta x \
      nabla g_2{(x)}^T Delta x\
      vdots \
      nabla g_n{(x)}^T Delta x
      end{bmatrix}$



      Translation:
      $g(Delta x)$ represents g with a small change in x as input, which is equivalent to $Delta g(x)$ [the infinitesimal change in g(x)].



      If you look at the last column vector, you will see that each row $nabla g_i(x)^T Delta x$ represents the (rate of change of $g_i(x)$/rate of change of x) dotted with (infinitesimal change in x) := infinitesimal change in $g_i(x)$. In aggregate, the column vector represents the infinitesimal change in g(x) := $Delta g(x)$.



      See also: http://www.friesian.com/calculus.htm






      share|cite|improve this answer











      $endgroup$
















        0












        0








        0





        $begingroup$

        For a scalar function f,



        $f(x+dx) approx f(x)+ frac{d f(x)}{dx}* dx$



        Vector case [$x$ is (n,1) dimensional vector]:



        $g(x + Delta x) approx g(x) + g(Delta x)$



        $g(x) = nabla f(x) = begin{bmatrix}
        g_1{(x)} \
        g_2{(x)} \
        vdots \
        g_n{(x)}
        end{bmatrix}$
        $mid g_i{(x)} = frac{partial f(x)}{dx_i}$



        $H(x) = nabla g(x) = begin{bmatrix}
        nabla g_1{(x)}^T \
        nabla g_2{(x)}^T \
        vdots \
        nabla g_n{(x)}^T
        end{bmatrix}$



        $g(x) +g(Delta x) = g(x) + Delta g(x)approx g(x)+ begin{bmatrix}
        nabla g_1{(x)}^T Delta x \
        nabla g_2{(x)}^T Delta x\
        vdots \
        nabla g_n{(x)}^T Delta x
        end{bmatrix}$



        Translation:
        $g(Delta x)$ represents g with a small change in x as input, which is equivalent to $Delta g(x)$ [the infinitesimal change in g(x)].



        If you look at the last column vector, you will see that each row $nabla g_i(x)^T Delta x$ represents the (rate of change of $g_i(x)$/rate of change of x) dotted with (infinitesimal change in x) := infinitesimal change in $g_i(x)$. In aggregate, the column vector represents the infinitesimal change in g(x) := $Delta g(x)$.



        See also: http://www.friesian.com/calculus.htm






        share|cite|improve this answer











        $endgroup$



        For a scalar function f,



        $f(x+dx) approx f(x)+ frac{d f(x)}{dx}* dx$



        Vector case [$x$ is (n,1) dimensional vector]:



        $g(x + Delta x) approx g(x) + g(Delta x)$



        $g(x) = nabla f(x) = begin{bmatrix}
        g_1{(x)} \
        g_2{(x)} \
        vdots \
        g_n{(x)}
        end{bmatrix}$
        $mid g_i{(x)} = frac{partial f(x)}{dx_i}$



        $H(x) = nabla g(x) = begin{bmatrix}
        nabla g_1{(x)}^T \
        nabla g_2{(x)}^T \
        vdots \
        nabla g_n{(x)}^T
        end{bmatrix}$



        $g(x) +g(Delta x) = g(x) + Delta g(x)approx g(x)+ begin{bmatrix}
        nabla g_1{(x)}^T Delta x \
        nabla g_2{(x)}^T Delta x\
        vdots \
        nabla g_n{(x)}^T Delta x
        end{bmatrix}$



        Translation:
        $g(Delta x)$ represents g with a small change in x as input, which is equivalent to $Delta g(x)$ [the infinitesimal change in g(x)].



        If you look at the last column vector, you will see that each row $nabla g_i(x)^T Delta x$ represents the (rate of change of $g_i(x)$/rate of change of x) dotted with (infinitesimal change in x) := infinitesimal change in $g_i(x)$. In aggregate, the column vector represents the infinitesimal change in g(x) := $Delta g(x)$.



        See also: http://www.friesian.com/calculus.htm







        share|cite|improve this answer














        share|cite|improve this answer



        share|cite|improve this answer








        edited Nov 28 '18 at 21:42

























        answered Nov 28 '18 at 9:48









        user3180user3180

        1185




        1185






























            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%2f3016823%2fhessian-vector-products%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

            Biblatex bibliography style without URLs when DOI exists (in Overleaf with Zotero bibliography)

            How to change which sound is reproduced for terminal bell?

            Can I use Tabulator js library in my java Spring + Thymeleaf project?