Solving quadratic Diophantine Equation $x_1^2+2x_2^2+3x_3^2=7y^2$. UPDATE: find all primitive solutions...

Question

After failing at stereographic projection, I opted to write things in terms of vectors. Definitely got hazy, but I was able to solve it. Currently, I haven't been able to come up with a parameterization that covers all solutions listed in Will Jagy's answer.

A good reference got me started: Equal Sums of Squares by C. J. Bradley Source: The Mathematical Gazette, Vol. 82, No. 493 (Mar., 1998), pp. 80-85 https://www.jstor.org/stable/3620159

Bradley Approach:
To solve $$x_1^2+2x_2^2+3x_3^2=7y_1^2 \tag{1}$$

First I went for: $$x_1^2+2x_2^2+3x_3^2=y_1^2+2y_2^2+3y_3^2 \tag{2}$$

this will correspond to $$ap+2bq+3cr=0$$ which becomes $$pf+2qg+3rh=0$$ with $(f,g,h)=1$, while $a=mf$, $b=mg$, $c=mh$ and $p, q, r$ being the cofactors of $$ \left( \begin{array}{ccc} f & 2g & 3h \\ t & u & v \end{array} \right)$$

So $$\begin{cases} p=2vg-3uh \\ q= -vf+3th\\ r=uf-2tg \end{cases}$$ Continuing $$\begin{cases} a=x_1-y_1 \\ p=x_1+y_1 \\ b=x_2-y_2 \\ q= x_2+y_2\\ c= x_3-y_3\\ r=x_3+y_3 \end{cases} \implies \begin{cases} x_1=\frac 1 2 (p+a)=\frac 1 2 (2vg-3uh+mf) \\ x_2=\frac 1 2 (q+b)=\frac 1 2 (-vf+3th+mg) \\ x_3 =\frac 1 2 (r+c) =\frac 1 2 (uf-2tg+mh)\\ y_1=\frac 1 2 (p-a)=\frac 1 2 (2vg-3uh-mf) =-\frac 1 2 (mf-2vg+3uh)\\ y_2= \frac 1 2 (q-b)=\frac 1 2 (-vf+3th-mg)=-\frac 1 2 (vf+mg-3th)\\ y_3=\frac 1 2 (r-c)=\frac 1 2 (uf-2tg-mh)=-\frac 1 2 (-uf+2tg+mh) \end{cases}$$

Eliminate the $1/2$ and multiply the $y$ values by $-1$. So this is $A$: $$x=Af \iff \left[ \begin{array}{c} x_1 \\ x_2 \\ x_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} m & 2v & -3u\\ -v & m & 3t\\ u & -2t & m\\ \end{array} \right]\left[ \begin{array}{c} f \\ g \\ h \\ \end{array} \right]$$

and this is $B$, $$y=Bf \iff \left[ \begin{array}{c} y_1 \\ y_2 \\ y_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} m & -2v & 3u\\ v & m & -3t\\ -u & 2t & m\\ \end{array} \right]\left[ \begin{array}{c} f \\ g \\ h \\ \end{array} \right]$$

$B$ is the negated triangle matrices that make up $A$

To test so far, let $m=1, f=-1, g=2, h=-2, t=3, u=-3, v=5$

$$\left[ \begin{array}{c} x_1 \\ x_2 \\ x_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} 1 & 10 & 9\\ -5 & 1 & 9\\ -3 & -6 & 1\\ \end{array} \right]\left[ \begin{array}{c} -1 \\ 2 \\ -2 \\ \end{array} \right]=\left[\begin{array}{c} 1 \\ -11\\ -11\end{array}\right]$$ $$\left[ \begin{array}{c} y_1 \\ y_2 \\ y_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} 1 & -10 & -9\\ 5 & 1 & -9\\ 3 & 6 & 1\\ \end{array} \right]\left[ \begin{array}{c} -1 \\ 2\\ -2 \\ \end{array} \right]=\left[\begin{array}{c} -3\\ 15\\ 7\end{array}\right]$$ And this IS a solution to $x_1^2+2x_2^2+3x_3^2=y_1^2+2y_2^2+3y_3^2$.

found that: $$x=Af \iff \left[ \begin{array}{c} x_1 \\ x_2 \\ x_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} m & 2v & -3u\\ -v & m & 3t\\ u & -2t & m\\ \end{array} \right]\left[ \begin{array}{c} f \\ g \\ h \\ \end{array} \right]$$

$$y=Bf \iff \left[ \begin{array}{c} y_1 \\ y_2 \\ y_3 \\ \end{array} \right]=\left[ \begin{array}{ccc} m & -2v & 3u\\ v & m & -3t\\ -u & 2t & m\\ \end{array} \right]\left[ \begin{array}{c} f \\ g \\ h \\ \end{array} \right]$$

I wrote $(2)$ in vector notation as: $x^t(Tx)=y^t(Ty)$ such that $T=\left[\begin{array}{ccc} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3\end{array}\right]$ So $$x^tT(Af)=y^tT(Bf)$$

Then, changing $f \to f_r$ by imposing a condition, $c=(2,0,1)$.

UPDATE: Right now the question stands as how can one find a parameterization containing all primitive solutions? The below was the answer I found with the approach I used. You will notice $(x_1, x_2)$ are permanently even. I was able to correct this, but my solution doesn't contain primitive forms of all solutions in Will Jagy's program.

$$\begin{cases} x_1=2(m^2+6t^2-3u^2-2v^2)+(-6mu+12tv) \\ x_2=2(-2mv+6tu)+(6mt+6uv) \\ x_3=2(2mu+4tv)+(m^2-6t^2-3u^2+2v^2) \\ y_1=m^2+6t^2+3u^2+2v^2 \end{cases}$$

Derivation: $x_1^2+2x_2^2+3x_3^2=y_1^2+2y_2^2+3y_3^2 \to (Af)^t(TAf)=(Bf)^t(TBf)$. If you set $Bf=c$, then $f=B^{-1}c$, and we have $$(AB^{-1}c)^t(TAB^{-1}c)=c^t(Tc)=2^2+2(0)^2+3(1)^2=7$$ Thus $$x=AB^{-1}c$$

Answer 1

Nevermind I found it $$\begin{cases} x_1=2(m^2+6t^2-3u^2-2v^2)+(-6mu+12tv) \\ x_2=2(-2mv+6tu)+(6mt+6uv) \\ x_3=2(2mu+4tv)+(m^2-6t^2-3u^2+2v^2) \\ y_1=m^2+6t^2+3u^2+2v^2 \end{cases}$$

Derivation: $x_1^2+2x_2^2+3x_3^2=y_1^2+2y_2^2+3y_3^2 \to (Af)^t(TAf)=(Bf)^t(TBf)$. If you set $Bf=c$, then $f=B^{-1}c$, and we have $$(AB^{-1}c)^t(TAB^{-1}c)=c^t(Tc)=2^2+2(0)^2+3(1)^2=7$$ I didn't interpret it right that here $x=AB^{-1}c$, I also didn't impose $c$ correctly in redefining $f$. When you work that out it is correct. It's just that calculating $AB^{-1}c$ was a doozy for me. I may have just made a small error the first time. Totally fine with that!

Answer 2

I'm afraid the permanent even variables are not doing the entire job:

     -7      1      2      3
-----------------------------
      1      2      0      1
      2      1      0      3
      2      5      0      1
      3      1      5      2
      3      2      4      3
      3      7      1      2
      5      1      9      2
      5      7      3      6
      5     10      6      1
      6      1      2      9
      6      5     10      3
      6      7      8      5
      6      7     10      1
      6     11      8      1
      6     13      2      5
      7      5      3     10
      7     10      0      9
      7     13      9      2
      7     14      6      5
      9      2     10     11
      9      5     11     10
      9     10     14      5
      9     11     13      6
      9     13      7     10
      9     14      2     11
      9     14      8      9
      9     19      7      6
      9     22      2      5
     10      5     18      3
     10      7     12     11
     10      7     18      1
     10     11      6     13
     10     13     12      9
     10     25      6      1
     11     10      6     15
     11     14     12     11
     11     14     18      1
     11     17     15      6
     11     22     12      5
     11     23      3     10
     13      1     21     10
     13      2     24      3
     13      7      9     18
     13     10      0     19
     13     17     21      2
     13     25     15      6
     13     34      0      3
     14      1     12     19
     14      7     18     15
     14     17      0     19
     14     19     18     11
     14     25      6     15
     14     29     12      9
     14     35      6      5
     14     37      0      1
     15      2     28      1
     15      5      7     22
     15      5     17     18
     15     10     14     19
     15     10     22     13
     15     11      1     22
     15     14     16     17
     15     14     26      3
     15     19     11     18
     15     22      2     19
     15     26      4     17
     15     26     14     13
     15     29     19      2
     15     34     14      3
     15     35     11      6
     15     35     13      2
     15     37      7      6
     15     38      8      1
     17      2     24     17
     17     11     15     22
     17     13     21     18
     17     14     24     15
     17     14     30      3
     17     29     21     10
     17     38      6     13
     17     43      9      2
     18      1     14     25
     18      5     28     15
     18      7      4     27
     18      7     26     17
     18     13     16     23
     18     17     14     23
     18     19      4     25
     18     25     22     15
     18     25     28      5
     18     29     26      5
     18     31     20     13
     18     35     20      9
     18     35     22      5
     18     37      4     17
     18     37     14     13
     18     41     16      5
     18     43     14      3
     18     47      4      3
     19      2      0     29
     19      2     18     25
     19      7     15     26
     19      7     33     10
     19     22     30      9
     19     25     15     22
     19     31     27      6
     19     34     12     19
     19     47      3     10
     19     49      3      6
     19     50      0      3
     21      1     23     26
     21      2     10     31
     21      7     13     30
     21      7     37     10
     21     10     20     27
     21     10     34     15
     21     14      2     31
     21     14     38      1
     21     17      7     30
     21     23     35      6
     21     25     35      2
     21     26     32     11
     21     31      7     26
     21     34     28     11
     21     38     22     15
     21     38     28      5
     21     46     22      1
     21     47     17     10
     21     49     17      6
     21     50     16      5
     21     55      5      2
     22      1     30     23
     22      7      6     33
     22      7     24     27
     22     11     36     15
     22     17     36     13
     22     19     24     25
     22     35     12     25
     22     35     30     11
     22     37     24     17
     22     47     24      3
     22     53      6     13
     22     55     12      5

==============================

Answer 3

Well, it is possible to solve $p^2 + 3 q^2 + 3 r^2 = 7 s^2$ using quaternions. The first recipe, using 1+3+3 = 7, comes out

---

? p = w^2 + x^2 - 3 * y^2 - 3 * z^2 - 0 *w*x - 6*w*y +6*w*z +6*x*y + 6*x*z +0*y*z 
%1 = x^2 + (6*y + 6*z)*x + (w^2 + (-6*y + 6*z)*w + (-3*y^2 - 3*z^2))
? q = w^2 - x^2 + 3 * y^2 - 3 * z^2 +2 *w*x - 0*w*y -2*w*z +2*x*y + 0*x*z +6*y*z 
%2 = -x^2 + (2*w + 2*y)*x + (w^2 - 2*z*w + (3*y^2 + 6*z*y - 3*z^2))
? r = w^2 - x^2 - 3 * y^2 + 3 * z^2 -2 *w*x +2*w*y -0*w*z +0*x*y + 2*x*z +6*y*z 
%3 = -x^2 + (-2*w + 2*z)*x + (w^2 + 2*y*w + (-3*y^2 + 6*z*y + 3*z^2))
? 
? s = w^2 + x^2 + 3*y^2 + 3 * z^2
%4 = x^2 + (w^2 + (3*y^2 + 3*z^2))
? 7* s^2
%5 = 7*x^4 + (14*w^2 + (42*y^2 + 42*z^2))*x^2 + (7*w^4 + (42*y^2 + 42*z^2)*w^2 + (63*y^4 + 126*z^2*y^2 + 63*z^4))
? 
? 
? p^2 + 3 * q^2 + 3 * r^2
%6 = 7*x^4 + (14*w^2 + (42*y^2 + 42*z^2))*x^2 + (7*w^4 + (42*y^2 + 42*z^2)*w^2 + (63*y^4 + 126*z^2*y^2 + 63*z^4))
? p^2 + 3 * q^2 + 3 * r^2 - 7 * s^2
%7 = 0
? 
? 

---

The second comes from writing 4+0+3 = 7.

? 
? p
%19 = 2*x^2 + 6*y*x + (2*w^2 + 6*z*w + (-6*y^2 - 6*z^2))
? 
? q
%20 = -x^2 + 4*y*x + (w^2 - 4*z*w + (3*y^2 - 3*z^2))
? 
? r
%21 = (-2*w + 4*z)*x + (4*y*w + 6*z*y)
? 
? s
%22 = x^2 + (w^2 + (3*y^2 + 3*z^2))
? 
? p^2 + 3 * q^2 + 3 * r^2 
%23 = 7*x^4 + (14*w^2 + (42*y^2 + 42*z^2))*x^2 + (7*w^4 + (42*y^2 + 42*z^2)*w^2 + (63*y^4 + 126*z^2*y^2 + 63*z^4))
? 
? 
? 7 * s^2
%24 = 7*x^4 + (14*w^2 + (42*y^2 + 42*z^2))*x^2 + (7*w^4 + (42*y^2 + 42*z^2)*w^2 + (63*y^4 + 126*z^2*y^2 + 63*z^4))
? 
? 
? p^2 + 3 * q^2 + 3 * r^2 - 7 * s^2
%25 = 0
? 

---

Answer 4

Above equation shown below:

$x_1^2+2x_2^2+3x_3^2=7y_1^2 \tag{1}$

"OP" & Will Jagy gave four variable solution to equation $(1)$

There is a one variable solution & is given below:

$x_1=(7k^2+8k-35)$

$x_2=(k+1)(k+13)$

$x_3=2(k^2-7k+1)$

$y_1=3(k^2+5)$

For, $k=0$ we get:

$(x_1,x_2,x_3,y_1)=(35,13,2,15)$

Answer 5

[Too long for a comment…]

Did you solve to this point? \begin{align} x_1 &= \frac{4g^2v^2-12ghuv-12ghmv-8g^2mv+9h^2u^2+18h^2mu+12ghmu-3h^2m^2-2g^2m^2}{ 2\left(2gv-3hu-3hm-2gm\right)}, \\[0.75em] x_2 &= \frac{4g^2v^2-12ghuv+4g^2mv+9h^2u^2-6ghmu+3h^2m^2-6ghm^2-2g^2m^2}{2\left(2gv-3hu-3hm-2gm\right)}, \\[0.75em] x_3 &= \frac{4g^2v^2-12ghuv+4ghmv+9h^2u^2-6h^2mu-3h^2m^2-4ghm^2+2g^2m^2}{2\left(2gv-3hu-3hm-2gm\right)}, \\[0.75em] y_1 &= \frac{4g^2v^2-12ghuv+9h^2u^2+3h^2m^2+2g^2m^2}{2\left(2gv-3hu-3hm-2gm\right)} \end{align}

It’s a rational [not-necessarily-integral] parameterization, but evidently allows different parities for the $x_i$. If you could absorb the denominator into the numerators, you’d have what you want!