vsip_dtoepsol_p - Solve a Toeplitz System of Equations

Core

This manual is preliminary and incomplete.
While our Core implementation implements all functions given in the standard we are still working on completing this documentation.

Please refer to the VSIPL standard for a complete function reference of the Core profile until we have completed work on this documentation.

[next] [prev] [prev-tail] [tail] [up]

6.2.1 vsip_dtoepsol_p - Solve a Toeplitz System of Equations

int vsip_toepsol_f(const vsip_vview_f* t, const vsip_vview_f* r, const vsip_vview_f* b, const vsip_vview_f* x); 
int vsip_ctoepsol_f(const vsip_cvview_f* t, const vsip_cvview_f* r, const vsip_cvview_f* b, const vsip_cvview_f* x);

Description

This function solves a real Toeplitz system of linear equations $Tx =b$ , where T is a symmetric Toeplitz matrix defined by its first column t and first row r, b is the right-hand side vector, and x is the solution vector. The Toeplitz matrix has constant diagonals, with the first column t and first row r defining the matrix structure.

Parameters

const vsip_dvview_p* t: Pointer to the vector view containing the first column of the Toeplitz matrix.
const vsip_dvview_p* r: Pointer to the vector view containing the first row of the Toeplitz matrix.
const vsip_dvview_p* b: Pointer to the vector view containing the right-hand side vector.
const vsip_dvview_p* x: Pointer to the vector view where the solution will be stored.

Return Value

Returns 0 on success.
Returns a non-zero value on error (e.g., if the Toeplitz matrix is singular).

Example

vsip_vview_f *toeplitz_col; 
vsip_vview_f *toeplitz_row; 
vsip_vview_f *rhs_vector; 
vsip_vview_f *solution_vector; 
int result; 
 
// Assuming all vector views have been properly initialized 
result = vsip_toepsol_f(toeplitz_col, toeplitz_row, rhs_vector, solution_vector); 
 
if (result != 0) { 
    // Handle error (e.g., singular matrix) 
}

[next] [prev] [prev-tail] [front] [up]