Function executeImageStep

src/library/blas/generic/solution_seq.c:252–421

}

static cl_int
executeImageStep(
    SolutionStep *step,
    cl_uint numEventsInWaitList,
    const cl_event *eventWaitList,
    cl_event *event)
{
    SolutionStep outerStep, innerStep, execStep;
    cl_int err = CL_SUCCESS;
    int currImg = 0;
    size_t imgWidth, imgHeight;
    size_t ha, hb;
    size_t maxPanels[MATRIX_ROLES_NUMBER], maxBlocks[MATRIX_ROLES_NUMBER];
    size_t off;
    SubproblemDim wholeDim;
    MatrixRole mrole;
    CLBlasKargs *kargs = &step->args;
    cl_mem *imgs = kargs->scimage;
    MemoryPattern *mempat = &clblasSolvers[step->funcID].memPatterns[step->patternID];
    ProblemIterator innerIter, outerIter;
    int oend = 0, iend;
    SDimComponent comp[2];
    bool backward;
    ListHead doneSteps;
    CLBlasKernelType ktype;

    kargsToProbDims(&wholeDim, step->funcID, kargs, false);
    memset(maxPanels, 0, sizeof(maxPanels));
    memset(maxBlocks, 0, sizeof(maxPanels));

    memcpy(&outerStep, step, sizeof(SolutionStep));
    memcpy(&execStep, step, sizeof(SolutionStep));
    listInitHead(&doneSteps);

    /*
     * Cover the whole problem with dimension which matrix blocks are
     * fitted to images at.
     */

    for (mrole = MATRIX_A; mrole < MATRIX_C; mrole++) {
        if (!isMatrixInImage(mempat, mrole)) {
            continue;
        }

        clGetImageInfo(imgs[currImg], CL_IMAGE_WIDTH, sizeof(imgWidth),
                       &imgWidth, NULL);
        clGetImageInfo(imgs[currImg], CL_IMAGE_HEIGHT, sizeof(imgHeight),
                       &imgHeight, NULL);

        if (step->funcID == CLBLAS_TRSM) {
            maxPanels[mrole] = 0;
            maxBlocks[mrole] = 0;
        } else {
            maxPanels[mrole] = imgHeight / matrBlockHeight(step->subdims, mrole,
                                                           clblasLeft);
        }
        currImg++;
    }