mphell_tuto_conversion_jacobi_quartic_to_weierstrass.c

Go to the documentation of this file.

/* 
  MPHELL-5.0 
  Author(s): The MPHELL team
 
 (C) Copyright 2015-2021 - Institut Fourier / Univ. Grenoble Alpes (France)
 
 This file is part of the MPHELL Library.
 MPHELL is free software: you can redistribute it and/or modify
 it under the terms of the GNU Lesser General Public License as published by
 the Free Software Foundation, version 3 of the License.
 
 MPHELL is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU Lesser General Public License for more details.
 
 You should have received a copy of the GNU Lesser General Public License
 along with MPHELL.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#include <stdio.h>
#include "mphell/mphell.h"
 
int main()
{
    /* Initialise MPHELL with 256 bits of security strength for the entropy, RANDOM_AES256 as DRBG and DEVURANDOM as entropy source */
    
    init_mphell(256, RANDOM_AES256, DEVURANDOM);
    
    /* Allocate a field of size 6*block_SIZE = 6*64 = 384 on 64 bits architecture */
    
    field k;
    field_alloc(k, FP, bits_to_nblock(384), NULL);
    
    /* Allocate temorary memory */
    
    number p;
    field_elt a, b;
    ec_point G, res, res_jq;
    number n, h;
    
    number_init(&p, bits_to_nblock(384));
 
    /* Set the number p from a string in base 16 */
    
    number_set_str(p, "fffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffeffffffff0000000000000000ffffffff", 16);
#if MPHELL_USE_AMNS == 1
    amns AMNS;
#if MPHELL_USE_AMNS_32 == 0
    amns_alloc_init_str(&AMNS, "[0, 7, [-2, 0, 0, 0, 0, 0, 0, 1], 59, 8010427042960414838884725536003025378504357914088932113335334539926383793809451548018236227209625700052722361381087, [-8576931003528147, 12963241998501462, -11685779024318585, 12198217901242760, 5469056521909520, -16952339206744690, -7801587214757611], [936600828956320511, 17864700281010200544, 7457034104920941375, 4906978974859566336, 1521310033682434027, 11977106235113612084, 6613199331606154624], [-48661600919137196, -16406260014740521, -9043881017757559, -15840327204754192, -4503903233538141, -10931797952018020, 3933578076642458], [-10700965224081570, -31747207631521220, -21134995866173053, -3985617456248275, -932744058138743, -25427794319760237, -2489368063126092]]", p);
#else
    amns_alloc_init_str(&AMNS, "[0, 18, [-2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1], 26, 28816340648335589686749390875939188706868735452517709312097740458242732054770156839198296353641579295263032233716200, [-51353, 1463368, 1158979, 544183, -726450, 976213, 144497, -886193, -264198, 688818, 1070182, -170086, 864472, -144241, -360825, -510521, 1017098, -59072], [3993315613, 3730700630, 867896355, 1312786307, 471958917, 4205419201, 2925107995, 1257333470, 1009072108, 3882493844, 334530385, 923966745, 2062755457, 1054721524, 1062123069, 2491481070, 813468425, 1351925765], [1499105, 1805368, 3513276, 2380622, 4551927, 3187369, 3810913, 2129326, 3493684, 1909307, 1981401, 2013148, 1633911, 2996996, 2975679, 3661385, 2455278, 1909836], [5854047, 6627253, 2246346, 3861327, 2251173, 2334257, 1163302, 3551220, 4498564, 4312125, 4221299, 3184071, 2952719, 307544, 1633107, 2468475, 2132615, 2940733]]", p);
#endif
    field_set_amns(k, AMNS);
#endif
 
    /* Create the field of characteristic p */
    
    field_create(k, "", STACK_1, 1, p);
 
    ec_point_alloc(res, k);
    ec_point_init(res, k);
    ec_point_alloc(res_jq, k);
    ec_point_init(res_jq, k);
    field_elt_alloc(&a, k);
    field_elt_init(a, k);
    field_elt_alloc(&b, k);
    field_elt_init(b, k);
    ec_point_alloc(G, k);
    ec_point_init(G, k);
    number_init(&n, bits_to_nblock(384));
    number_init(&h, bits_to_nblock(384));
 
    
    
    field_elt_set_str(a, "a3aeb56058fa59051926a1559fe934bf64545ee7dee231a8d65e7211702ac9bf9801cfdae54ad7a050e20e612d4c080", 16, false, k, STACK_1);
    field_elt_set_str(b, "1", 16, false, k, STACK_1);
    ec_point_set_str(G, "1541e4e6fec9431b561ad9d815306961cbda834477c2884b93bb5fb3a4c865d1638d6b1e1eb5d91afea443da57ec3c95", "e55245714090febf2f55c3dee76267c5aa27edb22a95ffa0328adecc2690ffef38c5cb63423401492f8e75bd4753f0b6","1","a7c8a5c9364fb8f2142801d10e8321c885c9fe7b88a7c76eabc474f635b090ac8001e9efcfc308ff06960a29869ebf8", false, 16, k, STACK_1);
    number_set_str(h, "15ad3c423d6df1572016199efaeee94ca5242a87a4", 16);
    number_set_str(n, "bcf569c46c00fa8b5f4c7d41b55e14295468dac4097918b25e24499", 16);
    
    /* Allocate curves */
    
    ec_curve E, E_jq;
    ec_alloc (E, k);
    ec_init(E, k);
    ec_alloc (E_jq, k);
    ec_init(E_jq, k);
    
    /* Create curve */
    
    ec_create (E_jq, "RandP384JQ", k, a, b, G, h, n, JACOBI_QUARTIC, PROJECTIVE, STACK_1);
    
    printf("The Jacobi Quartic departing curve \nE_jq: \n"); ec_curve_print(E_jq, 16, STACK_1); printf("\n");
 
    /* Creation of a random point on the Jacobi Quartic elliptic curve */
    
    jacobi_quartic_point_random(res_jq, E_jq, STACK_1);
    printf("random point, res_jq = "); ec_point_print(res_jq, 16, true, k, STACK_1); printf("\n");
    printf("res_jq belongs to E_jq : %d\n\n", jacobi_quartic_belongs(res_jq, E_jq, STACK_1));
 
    /* Conversion to the Weierstrass Elliptic curve */
    
    jacobi_quartic_to_weierstrass(E, E_jq, STACK_1);
 
    printf("The Weierstrass curve isomorphic to E_jq \n"); 
    printf("E: \n"); ec_curve_print(E, 16, STACK_1); printf("\n");
 
    /* Conversion of a point on a Weierstrass Elliptic curve*/
    
    jacobi_quartic_point_to_weierstrass_point (res, res_jq, E_jq, STACK_1);
 
    printf("Weierstrass random point, res = "); ec_point_print(res, 16, true, k, STACK_1); printf("\n");
    printf("res belongs to E : %d\n\n", weierstrass_belongs(res, E, STACK_1));
    
    /* Free allocated memory */
    
    field_elt_free(&a, k);
    field_elt_free(&b, k);
    ec_point_free(G, k);
    number_free(&n);
    number_free(&h);
    ec_point_free(res, k);
    ec_point_free(res_jq, k);
    number_free(&p);
#if MPHELL_USE_AMNS == 1
    amns_free(&AMNS);
#endif
    field_free(k);
    ec_free(E_jq); 
    ec_free(E);
    
    free_mphell();
 
    return 0;
}