ascon-x3-c64.c

Go to the documentation of this file.

/*
 * Copyright (C) 2022 Southern Storm Software, Pty Ltd.
 *
 * Permission is hereby granted, free of charge, to any person obtaining a
 * copy of this software and associated documentation files (the "Software"),
 * to deal in the Software without restriction, including without limitation
 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
 * and/or sell copies of the Software, and to permit persons to whom the
 * Software is furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included
 * in all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
 * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
 * DEALINGS IN THE SOFTWARE.
 */
 
#include "ascon-masked-state.h"
#include "ascon-masked-backend.h"
#include "core/ascon-util.h"
 
#if defined(ASCON_MASKED_X3_BACKEND_C64) && ASCON_MASKED_MAX_SHARES >= 3
 
#define and_not_xor(x, y, z) \
    do { \
        x##_a ^= (~(y##_a) & z##_a); \
        x##_a ^= ((y##_a) & ascon_mask64_unrotate_share1_0(z##_b)); \
        x##_a ^= (y##_a & ascon_mask64_unrotate_share2_0(z##_c)); \
        \
        x##_b ^= (y##_b & ascon_mask64_rotate_share1_0(z##_a)); \
        x##_b ^= ((~y##_b) & z##_b); \
        x##_b ^= (y##_b & ascon_mask64_unrotate_share2_1(z##_c)); \
        \
        x##_c ^= (y##_c & ascon_mask64_rotate_share2_0(~z##_a)); \
        x##_c ^= (y##_c & ascon_mask64_rotate_share2_1(z##_b)); \
        x##_c ^= (y##_c | z##_c); \
    } while (0)
 
/* Generate a pre-inverted round constant so that we can
 * avoid NOT'ing x2 in the S-box during the rounds */
#define ROUND_CONSTANT(round)   \
        (~(uint64_t)(((0x0F - (round)) << 4) | (round)))
 
void ascon_x3_permute
    (ascon_masked_state_t *state, uint8_t first_round, uint64_t preserve[2])
{
    static const uint64_t RC[12] = {
        ROUND_CONSTANT(0),
        ROUND_CONSTANT(1),
        ROUND_CONSTANT(2),
        ROUND_CONSTANT(3),
        ROUND_CONSTANT(4),
        ROUND_CONSTANT(5),
        ROUND_CONSTANT(6),
        ROUND_CONSTANT(7),
        ROUND_CONSTANT(8),
        ROUND_CONSTANT(9),
        ROUND_CONSTANT(10),
        ROUND_CONSTANT(11)
    };
    uint64_t x0_a, x1_a, x2_a, x3_a, x4_a;
    uint64_t x0_b, x1_b, x2_b, x3_b, x4_b;
    uint64_t x0_c, x1_c, x2_c, x3_c, x4_c;
    uint64_t t0_a, t0_b, t0_c, t1_a, t1_b, t1_c;
 
    /* Start with the randomness that the caller provided */
    t0_a = preserve[0];
    t0_b = preserve[1];
 
    /* Load the state into local variables */
#if defined(ASCON_MASKED_WORD_BACKEND_DIRECT_XOR)
    x0_a = be_load_word64(&(state->M[0].B[0]));
    x0_b = be_load_word64(&(state->M[0].B[8]));
    x0_c = be_load_word64(&(state->M[0].B[16]));
    x1_a = be_load_word64(&(state->M[1].B[0]));
    x1_b = be_load_word64(&(state->M[1].B[8]));
    x1_c = be_load_word64(&(state->M[1].B[16]));
    x2_a = be_load_word64(&(state->M[2].B[0]));
    x2_b = be_load_word64(&(state->M[2].B[8]));
    x2_c = be_load_word64(&(state->M[2].B[16]));
    x3_a = be_load_word64(&(state->M[3].B[0]));
    x3_b = be_load_word64(&(state->M[3].B[8]));
    x3_c = be_load_word64(&(state->M[3].B[16]));
    x4_a = be_load_word64(&(state->M[4].B[0]));
    x4_b = be_load_word64(&(state->M[4].B[8]));
    x4_c = be_load_word64(&(state->M[4].B[16]));
#else
    x0_a = state->M[0].S[0];
    x0_b = state->M[0].S[1];
    x0_c = state->M[0].S[2];
    x1_a = state->M[1].S[0];
    x1_b = state->M[1].S[1];
    x1_c = state->M[1].S[2];
    x2_a = state->M[2].S[0];
    x2_b = state->M[2].S[1];
    x2_c = state->M[2].S[2];
    x3_a = state->M[3].S[0];
    x3_b = state->M[3].S[1];
    x3_c = state->M[3].S[2];
    x4_a = state->M[4].S[0];
    x4_b = state->M[4].S[1];
    x4_c = state->M[4].S[2];
#endif
 
    /* The round constants invert x2 as part of the rounds so that we
     * don't need an explicit "x2 = ~x2" step in the S-box.  Pre-invert
     * x2 before the first round to compensate. */
    x2_a = ~x2_a;
 
    /* Perform all encryption rounds */
    while (first_round < 12) {
        /* Add the inverted round constant to x2 */
        x2_a ^= RC[first_round++];
 
        /* Start of the substitution layer, first share */
        x0_a ^= x4_a;
        x4_a ^= x3_a;
        x2_a ^= x1_a;
        t1_a  = x0_a;
 
        /* Start of the substitution layer, second share */
        x0_b ^= x4_b;
        x4_b ^= x3_b;
        x2_b ^= x1_b;
        t1_b  = x0_b;
 
        /* Start of the substitution layer, third share */
        x0_c ^= x4_c;
        x4_c ^= x3_c;
        x2_c ^= x1_c;
        t1_c  = x0_c;
 
        /* Middle part of the substitution layer, Chi5 */
        t0_c = ascon_mask64_rotate_share2_0(t0_a) ^ /* t0 = random shares */
               ascon_mask64_rotate_share2_1(t0_b);
        and_not_xor(t0, x0, x1);                    /* t0 ^= (~x0) & x1; */
        and_not_xor(x0, x1, x2);                    /* x0 ^= (~x1) & x2; */
        and_not_xor(x1, x2, x3);                    /* x1 ^= (~x2) & x3; */
        and_not_xor(x2, x3, x4);                    /* x2 ^= (~x3) & x4; */
        and_not_xor(x3, x4, t1);                    /* x3 ^= (~x4) & t1; */
        x4_a ^= t0_a;                               /* x4 ^= t0; */
        x4_b ^= t0_b;
        x4_c ^= t0_c;
 
        /* End of the substitution layer */
        x1_a ^= x0_a;
        x0_a ^= x4_a;
        x3_a ^= x2_a;
        x1_b ^= x0_b;
        x0_b ^= x4_b;
        x3_b ^= x2_b;
        x1_c ^= x0_c;
        x0_c ^= x4_c;
        x3_c ^= x2_c;
 
        /* NOT'ing x2 is done as part of the next round constant */
        /* x2_a = ~x2_a; */
 
        /* Linear diffusion layer, third share */
        x0_c ^= rightRotate19_64(x0_c) ^ rightRotate28_64(x0_c);
        x1_c ^= rightRotate61_64(x1_c) ^ rightRotate39_64(x1_c);
        x2_c ^= rightRotate1_64(x2_c)  ^ rightRotate6_64(x2_c);
        x3_c ^= rightRotate10_64(x3_c) ^ rightRotate17_64(x3_c);
        x4_c ^= rightRotate7_64(x4_c)  ^ rightRotate41_64(x4_c);
 
        /* Linear diffusion layer, second share */
        x0_b ^= rightRotate19_64(x0_b) ^ rightRotate28_64(x0_b);
        x1_b ^= rightRotate61_64(x1_b) ^ rightRotate39_64(x1_b);
        x2_b ^= rightRotate1_64(x2_b)  ^ rightRotate6_64(x2_b);
        x3_b ^= rightRotate10_64(x3_b) ^ rightRotate17_64(x3_b);
        x4_b ^= rightRotate7_64(x4_b)  ^ rightRotate41_64(x4_b);
 
        /* Linear diffusion layer, first share */
        x0_a ^= rightRotate19_64(x0_a) ^ rightRotate28_64(x0_a);
        x1_a ^= rightRotate61_64(x1_a) ^ rightRotate39_64(x1_a);
        x2_a ^= rightRotate1_64(x2_a)  ^ rightRotate6_64(x2_a);
        x3_a ^= rightRotate10_64(x3_a) ^ rightRotate17_64(x3_a);
        x4_a ^= rightRotate7_64(x4_a)  ^ rightRotate41_64(x4_a);
 
        /* Rotate the randomness in t0 before the next round */
        t0_a = rightRotate13_64(t0_a);
        t0_b = rightRotate29_64(t0_b);
    }
 
    /* Return the final randomness to the caller to preserve it */
    preserve[0] = t0_a;
    preserve[1] = t0_b;
 
    /* Store the local variables back to the state with a final invert of x2 */
#if defined(ASCON_MASKED_WORD_BACKEND_DIRECT_XOR)
    be_store_word64(&(state->M[0].B[0]),  x0_a);
    be_store_word64(&(state->M[0].B[8]),  x0_b);
    be_store_word64(&(state->M[0].B[16]), x0_c);
    be_store_word64(&(state->M[1].B[0]),  x1_a);
    be_store_word64(&(state->M[1].B[8]),  x1_b);
    be_store_word64(&(state->M[1].B[16]), x1_c);
    be_store_word64(&(state->M[2].B[0]),  ~x2_a);
    be_store_word64(&(state->M[2].B[8]),  x2_b);
    be_store_word64(&(state->M[2].B[16]), x2_c);
    be_store_word64(&(state->M[3].B[0]),  x3_a);
    be_store_word64(&(state->M[3].B[8]),  x3_b);
    be_store_word64(&(state->M[3].B[16]), x3_c);
    be_store_word64(&(state->M[4].B[0]),  x4_a);
    be_store_word64(&(state->M[4].B[8]),  x4_b);
    be_store_word64(&(state->M[4].B[16]), x4_c);
#else
    state->M[0].S[0] = x0_a;
    state->M[0].S[1] = x0_b;
    state->M[0].S[2] = x0_c;
    state->M[1].S[0] = x1_a;
    state->M[1].S[1] = x1_b;
    state->M[1].S[2] = x1_c;
    state->M[2].S[0] = ~x2_a;
    state->M[2].S[1] = x2_b;
    state->M[2].S[2] = x2_c;
    state->M[3].S[0] = x3_a;
    state->M[3].S[1] = x3_b;
    state->M[3].S[2] = x3_c;
    state->M[4].S[0] = x4_a;
    state->M[4].S[1] = x4_b;
    state->M[4].S[2] = x4_c;
#endif
}
 
#endif /* ASCON_MASKED_X3_BACKEND_C64 */