return (sqrt(dot(vec, beg, end, vec)));
}
-/* inline */
void orthog1(int n, double *vec /* vector to be orthogonalized against 1 */
)
{
#define RANGE 500
-/* inline */
void init_vec_orth1(int n, double *vec)
{
/* randomly generate a vector orthogonal to 1 (i.e., with mean 0) */
orthog1(n, vec);
}
-/* inline */
void
right_mult_with_vector(vtx_data * matrix, int n, double *vector,
double *result)
/* orthog1(n,vector); */
}
-/* inline */
void
right_mult_with_vector_f(float **matrix, int n, double *vector,
double *result)
/* orthog1(n,vector); */
}
-/* inline */
void
vectors_subtraction(int n, double *vector1, double *vector2,
double *result)
}
}
-/* inline */
void
vectors_addition(int n, double *vector1, double *vector2, double *result)
{
}
}
-/* inline */
void
vectors_scalar_mult(int n, double *vector, double alpha, double *result)
{
}
}
-/* inline */
void copy_vector(int n, double *source, double *dest)
{
int i;
dest[i] = source[i];
}
-/* inline */
double vectors_inner_product(int n, double *vector1, double *vector2)
{
int i;
return result;
}
-/* inline */
double max_abs(int n, double *vector)
{
double max_val = -1e50;
return max_val;
}
-/* inline */
void
right_mult_with_vector_transpose(double **matrix,
int dim1, int dim2,
}
}
-/* inline */
void
right_mult_with_vector_d(double **matrix,
int dim1, int dim2,
** version **
*****************************/
-/* inline */
void orthog1f(int n, float *vec)
{
int i;
}
}
-/* inline */
void right_mult_with_vector_ff
(float *packed_matrix, int n, float *vector, float *result) {
/* packed matrix is the upper-triangular part of a symmetric matrix arranged in a vector row-wise */
}
}
-/* inline */
void
vectors_substractionf(int n, float *vector1, float *vector2, float *result)
{
}
}
-/* inline */
void
vectors_additionf(int n, float *vector1, float *vector2, float *result)
{
}
}
-/* inline */
void
vectors_mult_additionf(int n, float *vector1, float alpha, float *vector2)
{
}
}
-/* inline */
void vectors_scalar_multf(int n, float *vector, float alpha, float *result)
{
int i;
}
}
-/* inline */
void copy_vectorf(int n, float *source, float *dest)
{
int i;
dest[i] = source[i];
}
-/* inline */
double vectors_inner_productf(int n, float *vector1, float *vector2)
{
int i;
return result;
}
-/* inline */
void set_vector_val(int n, double val, double *result)
{
int i;
result[i] = val;
}
-/* inline */
void set_vector_valf(int n, float val, float* result)
{
int i;
result[i] = val;
}
-/* inline */
double max_absf(int n, float *vector)
{
int i;
return max_val;
}
-/* inline */
void square_vec(int n, float *vec)
{
int i;
}
}
-/* inline */
void invert_vec(int n, float *vec)
{
int i;
}
}
-/* inline */
void sqrt_vec(int n, float *vec)
{
int i;
}
}
-/* inline */
void sqrt_vecf(int n, float *source, float *target)
{
int i;
}
}
-/* inline */
void invert_sqrt_vec(int n, float *vec)
{
int i;