Method handleMAD_AND_REDUCE

src/library/blas/gens/kprintf.cpp:773–833

}

void kprintf::handleMAD_AND_REDUCE(char **_src, char **_dst)
{
    int numCharsWritten = 0;
    char id1[256], id2[256], id3[256];
    char * ptr;
    char *src = *_src;
    char *dst = *_dst;
    int vwidth=1;

    ptr = mystrtok( src, "(,)");
    ptr = mystrtok( NULL, "(,)"); // Get first ID
    strcpy( id1, ptr);
    ptr = mystrtok( NULL, "(,)"); // Get second ID
    strcpy( id2, ptr);
    ptr = mystrtok( NULL, "(,)"); // Get third ID
    strcpy( id3, ptr);
    *_src = ptr + strlen(ptr) + 1;
    //std::cout << id1 << "  " << id2 << "  " << id3 << std::endl;
    if ( (strcmp(id1, id2) == 0) || (strcmp(id1, id3)==0) || (strcmp(id2,id3) == 0) )
    {
        std::cout << "%MAD_AND_REDUCE( C, A, B) : C , A and B have to be UNIQUE" << std::endl;
        throw -1;
    }

    switch(s_or_v)
    {
        case SCALAR:
            //
            // e.g. float += float4*float4
            //      We will use only the first vector component
            //
            if (vectorWidth == 1)
            {
                numCharsWritten = sprintf(dst, "%s = mad(%s,%s,%s);\n\t", id1, id2, id3, id1);
                dst += numCharsWritten;

            } else {
                numCharsWritten = sprintf(dst, "%s = mad(%s.%s,%s.%s,%s);\n\t", id1, id2, vecIndices[0], id3,
                                                                               vecIndices[0], id1);
                dst += numCharsWritten;
            }
            break;

        case VECTOR:
            numCharsWritten = sprintf(dst, "%s.S0 = mad((%s).S0,(%s).S0,%s.S0);\n\t", id1, id2, id3, id1);
            dst += numCharsWritten;
            numCharsWritten = sprintf(dst, "%s.S0 = mad((%s).S1,-(%s.S1),%s.S0);\n\t", id1, id2, id3, id1);
            dst += numCharsWritten;

            numCharsWritten = sprintf(dst, "%s.S1 = mad((%s).S0,(%s).S1,%s.S1);\n\t", id1, id2, id3, id1);
            dst += numCharsWritten;
            numCharsWritten = sprintf(dst, "%s.S1 = mad((%s).S1,(%s.S0),%s.S1);\n\t", id1, id2, id3, id1);
            dst += numCharsWritten;
            break;

        default:
            std::cout << "handleMAD_AND_REDUCE: s_or_v is neither scalar nor a vector" << std::endl;
            throw -1;