Sage code to check independence of rational points on elliptic curve












0












$begingroup$


Suppose I have three rational points $(x_1,y_1),(x_2,y_2),(x_3,y_3)$ on certain elliptic curve. Then they are linearly independent if and only if the determinant of matrix $(<P_i,P_j>)_{i,j}$ is non zero where $<_,_ >$ is Neron-Tate height pairing.



How to write a code in Sage or Pari to compute this determinant. Note that I don't want to write an equation of my elliptic curve in the code to compute rational points first. Rather I know my rational points and want to directly put rational points and compute the determinant.










share|cite|improve this question









$endgroup$








  • 2




    $begingroup$
    If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
    $endgroup$
    – Hw Chu
    Jan 20 at 20:39










  • $begingroup$
    @Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
    $endgroup$
    – ersh
    Jan 20 at 21:28










  • $begingroup$
    Try asking this question in the coding branch of Stack Exchange.
    $endgroup$
    – BadAtGeometry
    Jan 21 at 19:18
















0












$begingroup$


Suppose I have three rational points $(x_1,y_1),(x_2,y_2),(x_3,y_3)$ on certain elliptic curve. Then they are linearly independent if and only if the determinant of matrix $(<P_i,P_j>)_{i,j}$ is non zero where $<_,_ >$ is Neron-Tate height pairing.



How to write a code in Sage or Pari to compute this determinant. Note that I don't want to write an equation of my elliptic curve in the code to compute rational points first. Rather I know my rational points and want to directly put rational points and compute the determinant.










share|cite|improve this question









$endgroup$








  • 2




    $begingroup$
    If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
    $endgroup$
    – Hw Chu
    Jan 20 at 20:39










  • $begingroup$
    @Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
    $endgroup$
    – ersh
    Jan 20 at 21:28










  • $begingroup$
    Try asking this question in the coding branch of Stack Exchange.
    $endgroup$
    – BadAtGeometry
    Jan 21 at 19:18














0












0








0





$begingroup$


Suppose I have three rational points $(x_1,y_1),(x_2,y_2),(x_3,y_3)$ on certain elliptic curve. Then they are linearly independent if and only if the determinant of matrix $(<P_i,P_j>)_{i,j}$ is non zero where $<_,_ >$ is Neron-Tate height pairing.



How to write a code in Sage or Pari to compute this determinant. Note that I don't want to write an equation of my elliptic curve in the code to compute rational points first. Rather I know my rational points and want to directly put rational points and compute the determinant.










share|cite|improve this question









$endgroup$




Suppose I have three rational points $(x_1,y_1),(x_2,y_2),(x_3,y_3)$ on certain elliptic curve. Then they are linearly independent if and only if the determinant of matrix $(<P_i,P_j>)_{i,j}$ is non zero where $<_,_ >$ is Neron-Tate height pairing.



How to write a code in Sage or Pari to compute this determinant. Note that I don't want to write an equation of my elliptic curve in the code to compute rational points first. Rather I know my rational points and want to directly put rational points and compute the determinant.







algebraic-geometry elliptic-curves sagemath






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Jan 20 at 19:51









ershersh

423113




423113








  • 2




    $begingroup$
    If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
    $endgroup$
    – Hw Chu
    Jan 20 at 20:39










  • $begingroup$
    @Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
    $endgroup$
    – ersh
    Jan 20 at 21:28










  • $begingroup$
    Try asking this question in the coding branch of Stack Exchange.
    $endgroup$
    – BadAtGeometry
    Jan 21 at 19:18














  • 2




    $begingroup$
    If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
    $endgroup$
    – Hw Chu
    Jan 20 at 20:39










  • $begingroup$
    @Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
    $endgroup$
    – ersh
    Jan 20 at 21:28










  • $begingroup$
    Try asking this question in the coding branch of Stack Exchange.
    $endgroup$
    – BadAtGeometry
    Jan 21 at 19:18








2




2




$begingroup$
If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
$endgroup$
– Hw Chu
Jan 20 at 20:39




$begingroup$
If you do not provide the equation of the curve, how do you compare its Neron-Tate height? Given three points on the plane, you can always write down an elliptic curve which passes through all three which make the points dependent.
$endgroup$
– Hw Chu
Jan 20 at 20:39












$begingroup$
@Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
$endgroup$
– ersh
Jan 20 at 21:28




$begingroup$
@Hw Chu Yes. Mistake! So in general, how do we write a code in Sage to find that determinant. Can you please illustrate me though some example.
$endgroup$
– ersh
Jan 20 at 21:28












$begingroup$
Try asking this question in the coding branch of Stack Exchange.
$endgroup$
– BadAtGeometry
Jan 21 at 19:18




$begingroup$
Try asking this question in the coding branch of Stack Exchange.
$endgroup$
– BadAtGeometry
Jan 21 at 19:18










3 Answers
3






active

oldest

votes


















1












$begingroup$

In SageMath this can be accomplished using the command height_pairing_matrix. For example:



sage: E = EllipticCurve([1, 0, 0, 0, 1])
sage: E
Elliptic Curve defined by y^2 + x*y = x^3 + 1 over Rational Field
sage: pts = E.gens()
sage: pts
[(-1 : 1 : 1), (0 : -1 : 1)]
sage: M = E.height_pairing_matrix(points = pts)
sage: M
[0.541484699372469 0.161800998136433]
[0.161800998136433 0.463307138146013]
sage: M.det()
0.224694163418167


You can find more information, including many examples, in the documentation.






share|cite|improve this answer









$endgroup$













  • $begingroup$
    3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
    $endgroup$
    – ersh
    Jan 21 at 18:09






  • 1




    $begingroup$
    If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
    $endgroup$
    – André 3000
    Jan 21 at 18:12



















2












$begingroup$

Something similar was already asked at Independence of points on Elliptic curve the answers there are the same as above but maybe a little more extensive.



Edit: I just noticed the question asker is the same as the previous question. So maybe this question really is trying to ask something different. The height of a point will depend on the elliptic curve itself so you will need to put the equation of the curve somewhere.






share|cite|improve this answer











$endgroup$













  • $begingroup$
    Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
    $endgroup$
    – ersh
    Jan 21 at 19:46








  • 1




    $begingroup$
    If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
    $endgroup$
    – Alex J Best
    Jan 21 at 20:11












  • $begingroup$
    Thanks for more clarification.
    $endgroup$
    – ersh
    Jan 21 at 20:25



















1












$begingroup$

I am not a sage expert. But I know you can call pari within sage and here you can compute the Gram/height matrix from the points and determine its determinant. Suppose you e is your elliptic curve (in pari this created by ellinit) and suppose P1,P2,P3 are your points then you just use the following commands:



v=[P1,P2,P3];
m = ellheightmatrix(e,v);
print(matdet(m));


If the last line is $0$ then your points are dependent. Either something similar is done in sage or just call pari commands from sage.






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%2f3081066%2fsage-code-to-check-independence-of-rational-points-on-elliptic-curve%23new-answer', 'question_page');
    }
    );

    Post as a guest















    Required, but never shown

























    3 Answers
    3






    active

    oldest

    votes








    3 Answers
    3






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    1












    $begingroup$

    In SageMath this can be accomplished using the command height_pairing_matrix. For example:



    sage: E = EllipticCurve([1, 0, 0, 0, 1])
    sage: E
    Elliptic Curve defined by y^2 + x*y = x^3 + 1 over Rational Field
    sage: pts = E.gens()
    sage: pts
    [(-1 : 1 : 1), (0 : -1 : 1)]
    sage: M = E.height_pairing_matrix(points = pts)
    sage: M
    [0.541484699372469 0.161800998136433]
    [0.161800998136433 0.463307138146013]
    sage: M.det()
    0.224694163418167


    You can find more information, including many examples, in the documentation.






    share|cite|improve this answer









    $endgroup$













    • $begingroup$
      3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
      $endgroup$
      – ersh
      Jan 21 at 18:09






    • 1




      $begingroup$
      If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
      $endgroup$
      – André 3000
      Jan 21 at 18:12
















    1












    $begingroup$

    In SageMath this can be accomplished using the command height_pairing_matrix. For example:



    sage: E = EllipticCurve([1, 0, 0, 0, 1])
    sage: E
    Elliptic Curve defined by y^2 + x*y = x^3 + 1 over Rational Field
    sage: pts = E.gens()
    sage: pts
    [(-1 : 1 : 1), (0 : -1 : 1)]
    sage: M = E.height_pairing_matrix(points = pts)
    sage: M
    [0.541484699372469 0.161800998136433]
    [0.161800998136433 0.463307138146013]
    sage: M.det()
    0.224694163418167


    You can find more information, including many examples, in the documentation.






    share|cite|improve this answer









    $endgroup$













    • $begingroup$
      3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
      $endgroup$
      – ersh
      Jan 21 at 18:09






    • 1




      $begingroup$
      If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
      $endgroup$
      – André 3000
      Jan 21 at 18:12














    1












    1








    1





    $begingroup$

    In SageMath this can be accomplished using the command height_pairing_matrix. For example:



    sage: E = EllipticCurve([1, 0, 0, 0, 1])
    sage: E
    Elliptic Curve defined by y^2 + x*y = x^3 + 1 over Rational Field
    sage: pts = E.gens()
    sage: pts
    [(-1 : 1 : 1), (0 : -1 : 1)]
    sage: M = E.height_pairing_matrix(points = pts)
    sage: M
    [0.541484699372469 0.161800998136433]
    [0.161800998136433 0.463307138146013]
    sage: M.det()
    0.224694163418167


    You can find more information, including many examples, in the documentation.






    share|cite|improve this answer









    $endgroup$



    In SageMath this can be accomplished using the command height_pairing_matrix. For example:



    sage: E = EllipticCurve([1, 0, 0, 0, 1])
    sage: E
    Elliptic Curve defined by y^2 + x*y = x^3 + 1 over Rational Field
    sage: pts = E.gens()
    sage: pts
    [(-1 : 1 : 1), (0 : -1 : 1)]
    sage: M = E.height_pairing_matrix(points = pts)
    sage: M
    [0.541484699372469 0.161800998136433]
    [0.161800998136433 0.463307138146013]
    sage: M.det()
    0.224694163418167


    You can find more information, including many examples, in the documentation.







    share|cite|improve this answer












    share|cite|improve this answer



    share|cite|improve this answer










    answered Jan 21 at 18:04









    André 3000André 3000

    12.7k22243




    12.7k22243












    • $begingroup$
      3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
      $endgroup$
      – ersh
      Jan 21 at 18:09






    • 1




      $begingroup$
      If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
      $endgroup$
      – André 3000
      Jan 21 at 18:12


















    • $begingroup$
      3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
      $endgroup$
      – ersh
      Jan 21 at 18:09






    • 1




      $begingroup$
      If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
      $endgroup$
      – André 3000
      Jan 21 at 18:12
















    $begingroup$
    3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
    $endgroup$
    – ersh
    Jan 21 at 18:09




    $begingroup$
    3000, So using above code, sage will first generate these rational points, right? What If I already know my rational points, I don't want Sage to generate the points because it would take my time. How should I write a code in that case?
    $endgroup$
    – ersh
    Jan 21 at 18:09




    1




    1




    $begingroup$
    If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
    $endgroup$
    – André 3000
    Jan 21 at 18:12




    $begingroup$
    If you already have a list of the points, you can input it as points = in the height_pairing_matrix function. I tried to show this in my example; I could've just put E.height_pairing_matrix() and it would've computed the generators automatically.
    $endgroup$
    – André 3000
    Jan 21 at 18:12











    2












    $begingroup$

    Something similar was already asked at Independence of points on Elliptic curve the answers there are the same as above but maybe a little more extensive.



    Edit: I just noticed the question asker is the same as the previous question. So maybe this question really is trying to ask something different. The height of a point will depend on the elliptic curve itself so you will need to put the equation of the curve somewhere.






    share|cite|improve this answer











    $endgroup$













    • $begingroup$
      Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
      $endgroup$
      – ersh
      Jan 21 at 19:46








    • 1




      $begingroup$
      If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
      $endgroup$
      – Alex J Best
      Jan 21 at 20:11












    • $begingroup$
      Thanks for more clarification.
      $endgroup$
      – ersh
      Jan 21 at 20:25
















    2












    $begingroup$

    Something similar was already asked at Independence of points on Elliptic curve the answers there are the same as above but maybe a little more extensive.



    Edit: I just noticed the question asker is the same as the previous question. So maybe this question really is trying to ask something different. The height of a point will depend on the elliptic curve itself so you will need to put the equation of the curve somewhere.






    share|cite|improve this answer











    $endgroup$













    • $begingroup$
      Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
      $endgroup$
      – ersh
      Jan 21 at 19:46








    • 1




      $begingroup$
      If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
      $endgroup$
      – Alex J Best
      Jan 21 at 20:11












    • $begingroup$
      Thanks for more clarification.
      $endgroup$
      – ersh
      Jan 21 at 20:25














    2












    2








    2





    $begingroup$

    Something similar was already asked at Independence of points on Elliptic curve the answers there are the same as above but maybe a little more extensive.



    Edit: I just noticed the question asker is the same as the previous question. So maybe this question really is trying to ask something different. The height of a point will depend on the elliptic curve itself so you will need to put the equation of the curve somewhere.






    share|cite|improve this answer











    $endgroup$



    Something similar was already asked at Independence of points on Elliptic curve the answers there are the same as above but maybe a little more extensive.



    Edit: I just noticed the question asker is the same as the previous question. So maybe this question really is trying to ask something different. The height of a point will depend on the elliptic curve itself so you will need to put the equation of the curve somewhere.







    share|cite|improve this answer














    share|cite|improve this answer



    share|cite|improve this answer








    edited Jan 21 at 19:11

























    answered Jan 21 at 19:05









    Alex J BestAlex J Best

    2,24111226




    2,24111226












    • $begingroup$
      Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
      $endgroup$
      – ersh
      Jan 21 at 19:46








    • 1




      $begingroup$
      If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
      $endgroup$
      – Alex J Best
      Jan 21 at 20:11












    • $begingroup$
      Thanks for more clarification.
      $endgroup$
      – ersh
      Jan 21 at 20:25


















    • $begingroup$
      Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
      $endgroup$
      – ersh
      Jan 21 at 19:46








    • 1




      $begingroup$
      If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
      $endgroup$
      – Alex J Best
      Jan 21 at 20:11












    • $begingroup$
      Thanks for more clarification.
      $endgroup$
      – ersh
      Jan 21 at 20:25
















    $begingroup$
    Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
    $endgroup$
    – ersh
    Jan 21 at 19:46






    $begingroup$
    Yes, I was actually expecting you to comment as I had worked with your answer you provided me last time. You are right my question is different from the last one. I wanted to feed Sage only may rational points even including the equation of elliptic curve. But I don't want sage to compute the rational points. I already have a rational points in my hand and I want to compute determinant for my specific rational points.
    $endgroup$
    – ersh
    Jan 21 at 19:46






    1




    1




    $begingroup$
    If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
    $endgroup$
    – Alex J Best
    Jan 21 at 20:11






    $begingroup$
    If you give sage the equation of and elliptical curve it will not compute any rational points by default. You can write P = E((1,0)) to define a point with coordinates (1,0) on E if you know it. Then a list of such can be given to height_pairing_matrix . you will have go define the elliptic curve though as the matrix depends on the equation of the curve.
    $endgroup$
    – Alex J Best
    Jan 21 at 20:11














    $begingroup$
    Thanks for more clarification.
    $endgroup$
    – ersh
    Jan 21 at 20:25




    $begingroup$
    Thanks for more clarification.
    $endgroup$
    – ersh
    Jan 21 at 20:25











    1












    $begingroup$

    I am not a sage expert. But I know you can call pari within sage and here you can compute the Gram/height matrix from the points and determine its determinant. Suppose you e is your elliptic curve (in pari this created by ellinit) and suppose P1,P2,P3 are your points then you just use the following commands:



    v=[P1,P2,P3];
    m = ellheightmatrix(e,v);
    print(matdet(m));


    If the last line is $0$ then your points are dependent. Either something similar is done in sage or just call pari commands from sage.






    share|cite|improve this answer









    $endgroup$


















      1












      $begingroup$

      I am not a sage expert. But I know you can call pari within sage and here you can compute the Gram/height matrix from the points and determine its determinant. Suppose you e is your elliptic curve (in pari this created by ellinit) and suppose P1,P2,P3 are your points then you just use the following commands:



      v=[P1,P2,P3];
      m = ellheightmatrix(e,v);
      print(matdet(m));


      If the last line is $0$ then your points are dependent. Either something similar is done in sage or just call pari commands from sage.






      share|cite|improve this answer









      $endgroup$
















        1












        1








        1





        $begingroup$

        I am not a sage expert. But I know you can call pari within sage and here you can compute the Gram/height matrix from the points and determine its determinant. Suppose you e is your elliptic curve (in pari this created by ellinit) and suppose P1,P2,P3 are your points then you just use the following commands:



        v=[P1,P2,P3];
        m = ellheightmatrix(e,v);
        print(matdet(m));


        If the last line is $0$ then your points are dependent. Either something similar is done in sage or just call pari commands from sage.






        share|cite|improve this answer









        $endgroup$



        I am not a sage expert. But I know you can call pari within sage and here you can compute the Gram/height matrix from the points and determine its determinant. Suppose you e is your elliptic curve (in pari this created by ellinit) and suppose P1,P2,P3 are your points then you just use the following commands:



        v=[P1,P2,P3];
        m = ellheightmatrix(e,v);
        print(matdet(m));


        If the last line is $0$ then your points are dependent. Either something similar is done in sage or just call pari commands from sage.







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered Jan 20 at 22:10









        JoseJose

        29919




        29919






























            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%2f3081066%2fsage-code-to-check-independence-of-rational-points-on-elliptic-curve%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

            Can a sorcerer learn a 5th-level spell early by creating spell slots using the Font of Magic feature?

            Does disintegrating a polymorphed enemy still kill it after the 2018 errata?

            A Topological Invariant for $pi_3(U(n))$