vecview.c File Reference

Detailed Description

Multivector view implementation.

#include <assert.h>
#include <math.h>
#include <oski/config.h>
#include <oski/common.h>
#include <oski/vecview.h>

Defines

#define MIN(a, b)   ((a) < (b) ? (a) : (b))
 Min of two side-effect-free values.

Functions

static void SetRowColIncMultiVecView (oski_vecstruct_t *vec)
 Set the redundantly maintained row/column stride parameters.
static oski_vecview_t CreateMultiVecView_NoError (oski_value_t *X, oski_index_t length, oski_index_t num_vecs, oski_storage_t orient, oski_index_t stride)
 Create a multivector view object, but without error checking.
void oski_PrintDebugVecView (const oski_vecview_t x, const char *val_fmt)
 Dump a vector view object to the debug log (FOR DEBUGGING ONLY).
oski_vecview_t oski_CreateVecView (oski_value_t *x, oski_index_t length, oski_index_t inc)
 Creates a valid view on a single dense column vector $x$.
oski_vecview_t oski_CreateMultiVecView (oski_value_t *X, oski_index_t length, oski_index_t num_vecs, oski_storage_t orient, oski_index_t stride)
 Creates a multivector view on $k$ dense column vectors $X = \left(\begin{array}{c} x_1 \cdots x_k \end{array}\right)$, stored as a submatrix in the user's data.
int oski_DestroyVecView (oski_vecview_t x)
 Destroy an existing vector view.
oski_vecview_t oski_CopyVecView (const oski_vecview_t x)
 Creates a copy of the given (multi)vector view.
double oski_CalcVecViewInfNorm (const oski_vecview_t x, oski_index_t j)
 Compute the infinity norm of a specified column vector of a vector view.
double oski_CalcVecViewInfNormDiff (const oski_vecview_t x0, oski_index_t j0, const oski_vecview_t x1, oski_index_t j1)
 Compute the infinity norm of the difference between two vectors.
int oski_InitSubVecView (const oski_vecview_t x, oski_index_t i, oski_index_t j, oski_index_t R, oski_index_t C, oski_vecstruct_t *y)
 Given an $m\times n$ multivector object $x$, returns a shallow copy submatrix $y = x(i_0:(i_0+R-1), j_0:(j_0+C-1)$.

Variables

static oski_vecstruct_t OBJ_SYMBOLIC_VEC
 Statically declared symbolic single vector object.
static oski_vecstruct_t OBJ_SYMBOLIC_MULTIVEC
 Statically declared symbolic multivector object.
const oski_vecview_t SYMBOLIC_VEC = &OBJ_SYMBOLIC_VEC
 Symbolic single-vector view.
const oski_vecview_t SYMBOLIC_MULTIVEC = &OBJ_SYMBOLIC_MULTIVEC
 Symbolic multiple-vector view.

Function Documentation

static oski_vecview_t CreateMultiVecView_NoError (  oski_value_t *  X,
oski_index_t  length,
oski_index_t  num_vecs,
oski_storage_t  orient,
oski_index_t  stride
)  [static]
 

Create a multivector view object, but without error checking.

This routine is equivalent to oski_CreateMultiVecView(), but without the parameter checking.

See also:
oski_CreateMultiVecView

double oski_CalcVecViewInfNorm (  const oski_vecview_t  x,
oski_index_t  j
) 
 

Compute the infinity norm of a specified column vector of a vector view.

Parameters:
[in] x A valid vector view.
[in] j Column for which to compute the infinity norm (in 1-based indices).
Returns:
Let $x_j$ denote column $j$ of $x$. Then this routine returns $||x_j||_\infty$, or -1 on error.

double oski_CalcVecViewInfNormDiff (  const oski_vecview_t  x0,
oski_index_t  j0,
const oski_vecview_t  x1,
oski_index_t  j1
) 
 

Compute the infinity norm of the difference between two vectors.

Parameters:
[in] x0 A valid vector view.
[in] j0 Column of x0 to use (1-based).
[in] x1 A valid vector view.
[in] j1 Column of x1 to use (1-based).
Returns:
$||x0_{j0} - x1_{j1}||_\infty$, where $x0_{j0}$ is column $j0$ of $x0$, and $x1_{j1}$ is column $j1$ of $x1$. Returns -1 on error.

int oski_InitSubVecView (  const oski_vecview_t  x,
oski_index_t  i,
oski_index_t  j,
oski_index_t  R,
oski_index_t  C,
oski_vecstruct_t *  y
) 
 

Given an $m\times n$ multivector object $x$, returns a shallow copy submatrix $y = x(i_0:(i_0+R-1), j_0:(j_0+C-1)$.

Parameters:
[in] x The multivector object, $x$.
[in] i Starting row (1-based index).
[in] j Starting column (1-based index).
[in] R Number of rows in the submatrix, or 0 to include all rows.
[in] C Number of columns in the submatrix, or 0 to include all columns.
[in,out] y Pointer to a structure in which to store the shallow copy of the submatrix.
Returns:
0 on success, or an error code on error.

void oski_PrintDebugVecView (  const oski_vecview_t  x,
const char *  val_fmt
) 
 

Dump a vector view object to the debug log (FOR DEBUGGING ONLY).

Parameters:
[in] x Vector view object.
[in] val_fmt Format string to use when displaying values.
This format string MUST contain the following fields, in this order: "%f", "%c". That is, the implementation will print a single character to the right of the value. If val_mft is NULL, then the routines uses the default format of "%5.2f%c".

If oski_value_t is complex, then the format string must contain "%f", "%f", and "%c", to display the real and imaginary parts.

The extra fields are used to print special characters to distinguish between values of the logical matrix, and values lying in the "padding" when the stride is not equal to the matrix dimension.

static void SetRowColIncMultiVecView (  oski_vecstruct_t *  vec  )  [static]
 

Set the redundantly maintained row/column stride parameters.

Parameters:
[in,out] vec Partially initialized vector view (see precondition).
Precondition:
vec->orient and vec->stride have been initialized.
Postcondition:
Sets vec->row_stride and vec->col_stride.

Variable Documentation

oski_vecstruct_t OBJ_SYMBOLIC_MULTIVEC [static]
 

Initial value:

MAKE_VECSTRUCT_INIT (0, 1, LAYOUT_COLMAJ, 0, 1, 0, NULL)
Statically declared symbolic multivector object.

oski_vecstruct_t OBJ_SYMBOLIC_VEC [static]
 

Initial value:

MAKE_VECSTRUCT_INIT (0, 1, LAYOUT_COLMAJ, 0, 1, 0, NULL)
Statically declared symbolic single vector object.

Generated on Wed Sep 19 16:41:23 2007 for BeBOP Optimized Sparse Kernel Interface Library by  doxygen 1.4.6