static void typeError(const char *expected, const RtlFieldInfo *field)
{
VStringBuffer msg("mysql: type mismatch - %s expected", expected);
if (field)
msg.appendf(" for field %s", field->name->queryStr());
rtlFail(0, msg.str());
}
extern void failx(const char *message, ...)
{
va_list args;
va_start(args,message);
StringBuffer msg;
msg.appendf("%s: ", g_moduleName).valist_appendf(message,args);
va_end(args);
rtlFail(0, msg.str());
}
ECLBLAS_CALL void dpotf2(bool & __isAllResult, size32_t & __lenResult,
void * & __result, uint8_t tri, uint32_t r,
bool isAllA, size32_t lenA, const void * A,
bool clear) {
unsigned int cells = r*r;
__isAllResult = false;
__lenResult = cells * sizeof(double);
double *new_a = (double*) rtlMalloc(__lenResult);
memcpy(new_a, A, __lenResult);
double ajj;
// x and y refer to the embedded vectors for the multiply, not an axis
unsigned int diag, a_pos, x_pos, y_pos;
unsigned int col_step = r; // between columns
unsigned int row_step = 1; // between rows
unsigned int x_step = (tri==UPPER_TRIANGLE) ? row_step : col_step;
unsigned int y_step = (tri==UPPER_TRIANGLE) ? col_step : row_step;
for (unsigned int j=0; j<r; j++) {
diag = (j * r) + j; // diagonal
x_pos = j * ((tri==UPPER_TRIANGLE) ? col_step : row_step);
a_pos = (j+1) * ((tri==UPPER_TRIANGLE) ? col_step : row_step);
y_pos = diag + y_step;
// ddot.value <- x'*y
ajj = new_a[diag] - cblas_ddot(j, (new_a+x_pos), x_step, (new_a+x_pos), x_step);
//if ajj is 0, negative or NaN, then error
if (ajj <= 0.0) {
rtlFree(new_a);
rtlFail(0, "Not a positive definite matrix");
}
ajj = sqrt(ajj);
new_a[diag] = ajj;
if ( j < r-1) {
// y <- alpha*op(A)*x + beta*y
cblas_dgemv(CblasColMajor,
(tri==UPPER_TRIANGLE) ? CblasTrans : CblasNoTrans,
(tri==UPPER_TRIANGLE) ? j : r-1-j, // M
(tri==UPPER_TRIANGLE) ? r-1-j : j, // N
-1.0, // alpha
(new_a+a_pos), r, //A
(new_a+x_pos), x_step, //X
1.0, (new_a+y_pos), y_step); // beta and Y
// x <- alpha * x
cblas_dscal(r-1-j, 1.0/ajj, (new_a+y_pos), y_step);
}
// clear lower or upper part if clear flag set
for(unsigned int k=1; clear && k<r-j; k++) new_a[(k*x_step)+diag] = 0.0;
}
__result = (void*) new_a;
}
extern void fail(const char *message)
{
StringBuffer msg;
msg.appendf("%s: ", g_moduleName).append(message);
rtlFail(0, msg.str());
}
static void fail(const char *message)
{
StringBuffer msg;
msg.append("mysql: ").append(message);
rtlFail(0, msg.str());
}
