void
shape_conv(real h[], real shen[])
{
    int j;
    real h0 = h[0], h1 = h[1], h2 = h[2], h3 = h[3], h4 = h[4], tmp;
    for (j=0; j<NCWD; j++) {
        real energy=0;
        tmp = h0*cb_shape[j][0];
        energy += tmp*tmp;
        tmp = h0*cb_shape[j][1] + h1*cb_shape[j][0];
        energy += tmp*tmp;
        tmp = h0*cb_shape[j][2] + h1*cb_shape[j][1] + h2*cb_shape[j][0];
        energy += tmp*tmp;
        tmp = h0*cb_shape[j][3] + h1*cb_shape[j][2] + h2*cb_shape[j][1] +
              h3*cb_shape[j][0];
        energy += tmp*tmp;
        tmp = h0*cb_shape[j][4] + h1*cb_shape[j][3] + h2*cb_shape[j][2] +
              h3*cb_shape[j][1] + h4*cb_shape[j][0];
        energy += tmp*tmp;
        shen[j] = energy;
    }
}