/** Compute y <- alpha * x + y (elementwise addition) */
public static DoubleMatrix axpy(double da, DoubleMatrix dx, DoubleMatrix dy) {
// NativeBlas.daxpy(dx.length, da, dx.data, 0, 1, dy.data, 0, 1);
JavaBlas.raxpy(dx.length, da, dx.data, 0, 1, dy.data, 0, 1);
return dy;
}
/** Compute y <- alpha * x + y (elementwise addition) */
public static FloatMatrix axpy(float da, FloatMatrix dx, FloatMatrix dy) {
// NativeBlas.saxpy(dx.length, da, dx.data, 0, 1, dy.data, 0, 1);
JavaBlas.raxpy(dx.length, da, dx.data, 0, 1, dy.data, 0, 1);
return dy;
}