947 static void ecc_point_double_jacobian(u64 *x1, u64 *y1, u64 *z1,  in ecc_point_double_jacobian()  argument
954 	if (vli_is_zero(z1, ndigits))  in ecc_point_double_jacobian()
963 	/* t2 = y1*z1 = z3 */  in ecc_point_double_jacobian()
964 	vli_mod_mult_fast(y1, y1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
965 	/* t3 = z1^2 */  in ecc_point_double_jacobian()
966 	vli_mod_square_fast(z1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
968 	/* t1 = x1 + z1^2 */  in ecc_point_double_jacobian()
969 	vli_mod_add(x1, x1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
970 	/* t3 = 2*z1^2 */  in ecc_point_double_jacobian()
971 	vli_mod_add(z1, z1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
972 	/* t3 = x1 - z1^2 */  in ecc_point_double_jacobian()
973 	vli_mod_sub(z1, x1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
974 	/* t1 = x1^2 - z1^4 */  in ecc_point_double_jacobian()
975 	vli_mod_mult_fast(x1, x1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
977 	/* t3 = 2*(x1^2 - z1^4) */  in ecc_point_double_jacobian()
978 	vli_mod_add(z1, x1, x1, curve_prime, ndigits);  in ecc_point_double_jacobian()
979 	/* t1 = 3*(x1^2 - z1^4) */  in ecc_point_double_jacobian()
980 	vli_mod_add(x1, x1, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
989 	/* t1 = 3/2*(x1^2 - z1^4) = B */  in ecc_point_double_jacobian()
992 	vli_mod_square_fast(z1, x1, curve_prime, ndigits);  in ecc_point_double_jacobian()
994 	vli_mod_sub(z1, z1, t5, curve_prime, ndigits);  in ecc_point_double_jacobian()
996 	vli_mod_sub(z1, z1, t5, curve_prime, ndigits);  in ecc_point_double_jacobian()
998 	vli_mod_sub(t5, t5, z1, curve_prime, ndigits);  in ecc_point_double_jacobian()
1004 	vli_set(x1, z1, ndigits);  in ecc_point_double_jacobian()
1005 	vli_set(z1, y1, ndigits);  in ecc_point_double_jacobian()