Return to Answer

NB: I'm assuming from the way you worded the question that your '2D manifold' is supposed to be a surface in $xyz$-space and that the 'projection' is the projection to the $xy$-plane. You may have had a more general situation in mind, such as a surface in $\mathbb{R}^n$ for $n>2$ or a more general projection than the obvious one. In that case, the answer would change, but there is still a way of finding the answer. Let me know if you want to know about these more general situations.

You are asking a special case of an inverse problem in the calculus of variations, and there is a standard method for solving this problem. Note that, it is not true that every 2-parameter family $\mathscr{F}$ of curves in the plane (or an open domain in the plane) that has the property that there is a unique curve in the family passing through any given point in any given direction is the set of geodesics of some Riemannian metric, so it could easily be that most of the curve families $\mathscr{F}$ you might consider aren't the projections of geodesics of a metric of the form that you describe, i.e., a metric of the form $$ g = \sqrt{dx^2+dy^2 + df^2} $$ where $f = f(x,y)$ is a function on the given domain in the plane.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution.

Here, more explicitly, is what this amounts to: If your curve family is described by a second order equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) \\cr &\ \ \ \ - \bigl(2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\cr &\ \ \ \ +\bigl(2f_x(x,y)f_{xy}(x,y){-}f_y(x,y)f_{yy}(x,y)\bigr)(y')^2\\cr &\ \ \ \ +f_x(x,y)f_{yy}(x,y)(y')^3. \end{align} $$ Note that, one thing you'd know right off the bat is that if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient. To see more conditions, suppose that the curve family is the set of solutions of $$ y''=G(x,y,y')=a_0(x,y)+a_1(x,y)\ y'+a_2(x,y)\ (y')^2+a_3(x,y)\ (y')^3.\tag{1} $$ Substituting this into the above equation and equating like powers of $y'$ will yield $4$ second order equations for $f$ in terms of the functions $a_i$, namely $$ \begin{align} - f_y\ f_{xx} &= a_0\ (1{+}{f_x}^2{+}{f_y}^2)\\cr -2f_y\ f_{xy}{+}f_x\ f_{xx} &= a_1\ (1{+}{f_x}^2{+}{f_y}^2)\\cr 2f_x\ f_{xy}{-}f_y\ f_{yy} &= a_2\ (1{+}{f_x}^2{+}{f_y}^2)\\cr f_x\ f_{yy} &= a_3\ (1{+}{f_x}^2{+}{f_y}^2). \end{align}\tag{2} $$ For most choices of $a_i$ these equations are incompatible.

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions). In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this system either.

One useful observation for solving system (2) is that it implies the first order equation $$ a_0\ {f_x}^3 + a_1\ {f_x}^2f_y + a_2\ f_x{f_y}^2 + a_3\ {f_y}^3=0,\tag{3} $$ which shows that, when the functions $a_i$ are not all zero, the gradient of $f$ must point in one of three possible directions. Differentiating (3) and plugging this back into the system (2) and doing a bit more work, one can derive the necessary and sufficient conditions on the $a_i$ that there exist a solution $f$.

NB: I'm assuming from the way you worded the question that your '2D manifold' is supposed to be a surface in $xyz$-space and that the 'projection' is the projection to the $xy$-plane. You may have had a more general situation in mind, such as a surface in $\mathbb{R}^n$ for $n>2$ or a more general projection than the obvious one. In that case, the answer would change, but there is still a way of finding the answer. Let me know if you want to know about these more general situations.

You are asking a special case of an inverse problem in the calculus of variations, and there is a standard method for solving this problem. Note that, it is not true that every 2-parameter family $\mathscr{F}$ of curves in the plane (or an open domain in the plane) that has the property that there is a unique curve in the family passing through any given point in any given direction is the set of geodesics of some Riemannian metric, so it could easily be that most of the curve families $\mathscr{F}$ you might consider aren't the projections of geodesics of a metric of the form that you describe, i.e., a metric of the form $$ g = \sqrt{dx^2+dy^2 + df^2} $$ where $f = f(x,y)$ is a function on the given domain in the plane.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution.

Here, more explicitly, is what this amounts to: If your curve family is described by a second order equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) \\ &\ \ \ \ - \bigl(2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\ &\ \ \ \ +\bigl(2f_x(x,y)f_{xy}(x,y){-}f_y(x,y)f_{yy}(x,y)\bigr)(y')^2\\ &\ \ \ \ +f_x(x,y)f_{yy}(x,y)(y')^3. \end{align} $$ Note that, one thing you'd know right off the bat is that if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient. To see more conditions, suppose that the curve family is the set of solutions of $$ y''=G(x,y,y')=a_0(x,y)+a_1(x,y)\ y'+a_2(x,y)\ (y')^2+a_3(x,y)\ (y')^3.\tag{1} $$ Substituting this into the above equation and equating like powers of $y'$ will yield $4$ second order equations for $f$ in terms of the functions $a_i$, namely $$ \begin{align} - f_y\ f_{xx} &= a_0\ (1{+}{f_x}^2{+}{f_y}^2)\\ -2f_y\ f_{xy}{+}f_x\ f_{xx} &= a_1\ (1{+}{f_x}^2{+}{f_y}^2)\\ 2f_x\ f_{xy}{-}f_y\ f_{yy} &= a_2\ (1{+}{f_x}^2{+}{f_y}^2)\\ f_x\ f_{yy} &= a_3\ (1{+}{f_x}^2{+}{f_y}^2). \end{align}\tag{2} $$ For most choices of $a_i$ these equations are incompatible.

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions). In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this system either.

One useful observation for solving system (2) is that it implies the first order equation $$ a_0\ {f_x}^3 + a_1\ {f_x}^2f_y + a_2\ f_x{f_y}^2 + a_3\ {f_y}^3=0,\tag{3} $$ which shows that, when the functions $a_i$ are not all zero, the gradient of $f$ must point in one of three possible directions. Differentiating (3) and plugging this back into the system (2) and doing a bit more work, one can derive the necessary and sufficient conditions on the $a_i$ that there exist a solution $f$.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution. Note that, one thing you'd know right off the bat is that any $F$ that you get from the above process must be at most cubic in $y'$, so if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient.

Here, more explicitly, is what this amounts to: If your curve family is described by an equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) \\cr &\ \ \ \ - (2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\cr &\ \ \ \ +\bigl(2f_x(x,y)f_{xy}(x,y){-}f_y(x,y)f_{yy}(x,y)\bigr)(y')^2\\cr &\ \ \ \ +f_x(x,y)f_{yy}(x,y)(y')^3 \end{align} $$

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions).

  In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this system either.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution.

Here, more explicitly, is what this amounts to: If your curve family is described by a second order equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) \\cr &\ \ \ \ - \bigl(2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\cr &\ \ \ \ +\bigl(2f_x(x,y)f_{xy}(x,y){-}f_y(x,y)f_{yy}(x,y)\bigr)(y')^2\\cr &\ \ \ \ +f_x(x,y)f_{yy}(x,y)(y')^3. \end{align} $$ Note that, one thing you'd know right off the bat is that if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient. To see more conditions, suppose that the curve family is the set of solutions of $$ y''=G(x,y,y')=a_0(x,y)+a_1(x,y)\ y'+a_2(x,y)\ (y')^2+a_3(x,y)\ (y')^3.\tag{1} $$ Substituting this into the above equation and equating like powers of $y'$ will yield $4$ second order equations for $f$ in terms of the functions $a_i$, namely $$ \begin{align} - f_y\ f_{xx} &= a_0\ (1{+}{f_x}^2{+}{f_y}^2)\\cr -2f_y\ f_{xy}{+}f_x\ f_{xx} &= a_1\ (1{+}{f_x}^2{+}{f_y}^2)\\cr 2f_x\ f_{xy}{-}f_y\ f_{yy} &= a_2\ (1{+}{f_x}^2{+}{f_y}^2)\\cr f_x\ f_{yy} &= a_3\ (1{+}{f_x}^2{+}{f_y}^2). \end{align}\tag{2} $$ For most choices of $a_i$ these equations are incompatible.

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions). In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this system either.

One useful observation for solving system (2) is that it implies the first order equation $$ a_0\ {f_x}^3 + a_1\ {f_x}^2f_y + a_2\ f_x{f_y}^2 + a_3\ {f_y}^3=0,\tag{3} $$ which shows that, when the functions $a_i$ are not all zero, the gradient of $f$ must point in one of three possible directions. Differentiating (3) and plugging this back into the system (2) and doing a bit more work, one can derive the necessary and sufficient conditions on the $a_i$ that there exist a solution $f$.

NB: I'm assuming from the way you worded the question that your '2D manifold' is supposed to be a surface in $xyz$-space and that the 'projection' is the projection to the $xy$-plane. You may have had a more general situation in mind, such as a surface in $\mathbb{R}^n$ for $n>2$ or a more general projection than the obvious one. In that case, the answer would change, but there is still a way of finding the answer. Let me know if you want to know about these more general situations.

You are asking a special case of an inverse problem in the calculus of variations, and there is a standard method for solving this problem. Note that, it is not true that every 2-parameter family $\mathscr{F}$ of curves in the plane (or an open domain in the plane) that has the property that there is a unique curve in the family passing through any given point in any given direction is the set of geodesics of some Riemannian metric, so it could easily be that most of the curve families $\mathscr{F}$ you might consider aren't the projections of geodesics of a metric of the form that you describe, i.e., a metric of the form $$ g = \sqrt{dx^2+dy^2 + df^2} $$ where $f = f(x,y)$ is a function on the given domain in the plane.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution. Note that, one thing you'd know right off the bat is that any $F$ that you get from the above process must be at most cubic in $y'$, so if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient.

Here, more explicitly, is what this amounts to: If your curve family is described by an equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) - (2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\cr &\ +\bigl(2f_x(x,y)f_{xy}(x,y)-f_y(x,y)f_{yy}(x,y)\bigr)(y')^2+f_x(x,y)f_{yy}(x,y)(y')^3 \end{align} $$

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions).

In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this equation either.

NB: I'm assuming from the way you worded the question that your '2D manifold' is supposed to be a surface in $xyz$-space and that the 'projection' is the projection to the $xy$-plane. You may have had a more general situation in mind, such as a surface in $\mathbb{R}^n$ for $n>2$ or a more general projection than the obvious one. In that case, the answer would change, but there is still a way of finding the answer. Let me know if you want to know about these more general situations.

You are asking a special case of an inverse problem in the calculus of variations, and there is a standard method for solving this problem. Note that, it is not true that every 2-parameter family $\mathscr{F}$ of curves in the plane (or an open domain in the plane) that has the property that there is a unique curve in the family passing through any given point in any given direction is the set of geodesics of some Riemannian metric, so it could easily be that most of the curve families $\mathscr{F}$ you might consider aren't the projections of geodesics of a metric of the form that you describe, i.e., a metric of the form $$ g = \sqrt{dx^2+dy^2 + df^2} $$ where $f = f(x,y)$ is a function on the given domain in the plane.

The reconstruction method involves computing the Euler-Lagrange equation for the above $g$ as a Lagrangian for curves in the plane and writing it in the form $$ y'' = F(x,y,y') $$ for a certain function $F$. Then you set this $F$ equal to the equation that defines your curve family. For example, for the family $y = c_1 + c_2 e^x$, you'd have $y'' = y'= F(x,y,y')$. This will give you an overdetermined equation for $f$ and then you apply the usual theory of overdetermined systems to determine whether or not there is a solution. Note that, one thing you'd know right off the bat is that any $F$ that you get from the above process must be at most cubic in $y'$, so if your defining equation for the curve family isn't of the form $y'' = G(x,y,y')$ where $G$ is at most cubic in $y'$, then there will be no solution. This is only a first necessary condition, though, it is not sufficient.

Here, more explicitly, is what this amounts to: If your curve family is described by an equation of the form $y'' = G(x,y,y')$, then any function $f(x,y)$ on a domain such that the above metric has the curve family as geodesics must satisfy the nonlinear equation $$ \begin{align} (1{+}f_x(x,y)^2{+}f_y(x,y)^2)G(x,y,y') &= - f_y(x,y)f_{xx}(x,y) \\cr &\ \ \ \ - (2f_y(x,y)f_{xy}(x,y){-}f_x(x,y)f_{xx}(x,y)\bigr)y'\\cr &\ \ \ \ +\bigl(2f_x(x,y)f_{xy}(x,y){-}f_y(x,y)f_{yy}(x,y)\bigr)(y')^2\\cr &\ \ \ \ +f_x(x,y)f_{yy}(x,y)(y')^3 \end{align} $$

Thus, for example, if $G\equiv0$, $f$ must satisfy the four second order equations $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0. $$ Of course, this implies that, at any point where $(f_x,f_y)$ is nonzero, one must have $f_{xx}=f_{xy}=f_{yy}=0$, so the only solutions to these equations are to have $f$ be a linear function of $x$ and $y$ (and all linear functions are solutions).

In the case in which $G\equiv 2$ (the monic quadratic polynomials), the equations become $$ f_y\ f_{xx} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_x\ f_{xx}-2f_y\ f_{xy} = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is not difficult to see that there are no solutions to this overdetermined system. In the case in which $G\equiv y'$ (your other case), the equations become $$ f_y\ f_{xx} = f_x\ f_{xx}-2f_y\ f_{xy} + 2(1{+}{f_x}^2{+}{f_y}^2) = f_y\ f_{yy}-2f_x\ f_{xy} = f_x\ f_{yy} = 0, $$ and it is easy to see that there are no solutions to this system either.