Progress output in MAGMA

[evanlezar]

Hi,

We use MAGMA predominantly for the LU decomposition.

One nice to have feature would be support for displaying some form of progress for the computational processes -- something like the mkl_progress or acml_progress available in mkl and acml. This does not have to be particularly fine-grained.

Is it feasible to consider something like this?

I don't mind doing some of the legwork -- at least for the LU decomposition implementation -- but would appreciate some input on the best way to go about this.

Regards
Evan

[mgates3]

You can accomplish this by adding a call at the end of the loop and nb0 if-block in src/*getrf.cpp. E.g.,

Code: Select all

        for( i = 0; i < s; i++ )
        {
            [ ... ]
            your_progress( (i+1)*nb );
        }
        magma_int_t nb0 = min(m - s*nb, n - s*nb);
        if ( nb0 > 0 ) {
            [ ... ]
            your_progress( n );
        }


Aborting the computation is harder, as you would have to free various resources. If you use non-resident or multi-GPU codes, those routines would have to be modified as well.

I don't know if we would be interested to add this as a general feature in MAGMA, though. We will consider it.

-mark

[evanlezar]

Thanks mark.

I will add the progress code as you suggest. For now I am not concerned with terminating the calculation.

I will add the code with the relevant defines and submit a patch to see if it is to you liking.

Regards