mul_add
mul_add(rp[0], ap[0], w, c1);
mul_add(rp[1], ap[1], w, c1);
mul_add(rp[2], ap[2], w, c1);
mul_add(rp[3], ap[3], w, c1);
mul_add(rp[0], ap[0], w, c1);
mul_add(rp[1], ap[1], w, c1);
mul_add(rp[2], ap[2], w, c1);
mul_add(rp[0], ap[0], bl, bh, c);
mul_add(rp[1], ap[1], bl, bh, c);
mul_add(rp[2], ap[2], bl, bh, c);
mul_add(rp[3], ap[3], bl, bh, c);
mul_add(rp[0], ap[0], bl, bh, c);
mul_add(rp[0], ap[0], w, c1);
mul_add(rp[1], ap[1], w, c1);
mul_add(rp[2], ap[2], w, c1);
mul_add(rp[3], ap[3], w, c1);
mul_add(rp[0], ap[0], w, c1);
mul_add(tp[j], ap[j], ml, mh, c0);
mul_add(tp[j], ap[j], ml, c0);
mul_add(c1, np[0], ml, mh, c0);
mul_add(c1, ml, np[0], c0);
mul_add(c1, np[j], ml, mh, c0);
mul_add(c1, ml, np[j], c0);