Line data    Source code

       1             : /*
       2             :  * Copyright (C) 2013 Andrea Mazzoleni
       3             :  *
       4             :  * This program is free software: you can redistribute it and/or modify
       5             :  * it under the terms of the GNU General Public License as published by
       6             :  * the Free Software Foundation, either version 2 of the License, or
       7             :  * (at your option) any later version.
       8             :  *
       9             :  * This program is distributed in the hope that it will be useful,
      10             :  * but WITHOUT ANY WARRANTY; without even the implied warranty of
      11             :  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
      12             :  * GNU General Public License for more details.
      13             :  */
      14             : 
      15             : #include "internal.h"
      16             : #include "gf.h"
      17             : 
      18             : /*
      19             :  * GENz (triple parity with powers of 2^-1) 32bit C implementation
      20             :  */
      21          21 : void raid_genz_int32(int nd, size_t size, void **vv)
      22             : {
      23          21 :         uint8_t **v = (uint8_t**)vv;
      24             :         uint8_t *p;
      25             :         uint8_t *q;
      26             :         uint8_t *r;
      27             :         int d, l;
      28             :         size_t i;
      29             : 
      30             :         uint32_t d0, r0, q0, p0;
      31             :         uint32_t d1, r1, q1, p1;
      32             : 
      33          21 :         l = nd - 1;
      34          21 :         p = v[nd];
      35          21 :         q = v[nd + 1];
      36          21 :         r = v[nd + 2];
      37             : 
      38      655413 :         for (i = 0; i < size; i += 8) {
      39      655392 :                 r0 = q0 = p0 = v_32(v[l][i]);
      40      655392 :                 r1 = q1 = p1 = v_32(v[l][i + 4]);
      41     5243904 :                 for (d = l - 1; d >= 0; --d) {
      42     4588512 :                         d0 = v_32(v[d][i]);
      43     4588512 :                         d1 = v_32(v[d][i + 4]);
      44             : 
      45     4588512 :                         p0 ^= d0;
      46     4588512 :                         p1 ^= d1;
      47             : 
      48     4588512 :                         q0 = x2_32(q0);
      49     4588512 :                         q1 = x2_32(q1);
      50             : 
      51     4588512 :                         q0 ^= d0;
      52     4588512 :                         q1 ^= d1;
      53             : 
      54     4588512 :                         r0 = d2_32(r0);
      55     4588512 :                         r1 = d2_32(r1);
      56             : 
      57     4588512 :                         r0 ^= d0;
      58     4588512 :                         r1 ^= d1;
      59             :                 }
      60      655392 :                 v_32(p[i]) = p0;
      61      655392 :                 v_32(p[i + 4]) = p1;
      62      655392 :                 v_32(q[i]) = q0;
      63      655392 :                 v_32(q[i + 4]) = q1;
      64      655392 :                 v_32(r[i]) = r0;
      65      655392 :                 v_32(r[i + 4]) = r1;
      66             :         }
      67          21 : }
      68             : 
      69             : /*
      70             :  * GENz (triple parity with powers of 2^-1) 64bit C implementation
      71             :  */
      72          38 : void raid_genz_int64(int nd, size_t size, void **vv)
      73             : {
      74          38 :         uint8_t **v = (uint8_t**)vv;
      75             :         uint8_t *p;
      76             :         uint8_t *q;
      77             :         uint8_t *r;
      78             :         int d, l;
      79             :         size_t i;
      80             : 
      81             :         uint64_t d0, r0, q0, p0;
      82             :         uint64_t d1, r1, q1, p1;
      83             : 
      84          38 :         l = nd - 1;
      85          38 :         p = v[nd];
      86          38 :         q = v[nd + 1];
      87          38 :         r = v[nd + 2];
      88             : 
      89      606262 :         for (i = 0; i < size; i += 16) {
      90      606224 :                 r0 = q0 = p0 = v_64(v[l][i]);
      91      606224 :                 r1 = q1 = p1 = v_64(v[l][i + 8]);
      92     4850176 :                 for (d = l - 1; d >= 0; --d) {
      93     4243952 :                         d0 = v_64(v[d][i]);
      94     4243952 :                         d1 = v_64(v[d][i + 8]);
      95             : 
      96     4243952 :                         p0 ^= d0;
      97     4243952 :                         p1 ^= d1;
      98             : 
      99     4243952 :                         q0 = x2_64(q0);
     100     4243952 :                         q1 = x2_64(q1);
     101             : 
     102     4243952 :                         q0 ^= d0;
     103     4243952 :                         q1 ^= d1;
     104             : 
     105     4243952 :                         r0 = d2_64(r0);
     106     4243952 :                         r1 = d2_64(r1);
     107             : 
     108     4243952 :                         r0 ^= d0;
     109     4243952 :                         r1 ^= d1;
     110             :                 }
     111      606224 :                 v_64(p[i]) = p0;
     112      606224 :                 v_64(p[i + 8]) = p1;
     113      606224 :                 v_64(q[i]) = q0;
     114      606224 :                 v_64(q[i + 8]) = q1;
     115      606224 :                 v_64(r[i]) = r0;
     116      606224 :                 v_64(r[i + 8]) = r1;
     117             :         }
     118          38 : }
     119             :