Commit 438029ae authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Redefine magma_izamax() returned indices as 0-based

parent a47e3dbc

src/libnptm/magma_calls.cpp

+4 −4
Original line number Diff line number Diff line
@@ -408,9 +408,9 @@ magma_int_t magma_newton( 
    exit(1);

  }

  double max_residue, target_residue;

  magma_int_t maxindex = magma_izamax(mm, d_a, 1, queue);

  magma_int_t maxindex = magma_izamax(mm, d_a, 1, queue) - 1;

  magmaDoubleComplex magmamax = magma_mone;

  magma_zgetvector(1, d_a + maxindex - 1, 1, &magmamax, 1, queue);

  magma_zgetvector(1, d_a + maxindex, 1, &magmamax, 1, queue);

  curmax = MAGMA_Z_ABS(magmamax); //cabs(magmamax.x + I * magmamax.y);

  target_residue = curmax * rs.accuracy_goal;

  sprintf(buffer, "INFO: largest matrix value has modulus %.5le; target residue is %.5le.\n", curmax, target_residue);
@@ -425,8 +425,8 @@ magma_int_t magma_newton( 
    );

    // Transform -A*X into (I - A*X)

    magma_zaxpy(m, magma_one, d_id_diag, 1, d_ax, m + 1, queue);

    maxindex = magma_izamax(mm, d_ax, 1, queue);

    magma_zgetvector(1, d_ax + maxindex - 1, 1, &magmamax, 1, queue);

    maxindex = magma_izamax(mm, d_ax, 1, queue) - 1;

    magma_zgetvector(1, d_ax + maxindex, 1, &magmamax, 1, queue);

    curmax = cabs(magmamax.x + I * magmamax.y);

    sprintf(buffer, "DEBUG: iteration %d has residue %.5le; target residue is %.5le.\n", ri, curmax, target_residue);

    message = buffer;