A Survey of Gaussian Convolution Algorithms
Functions
invert_matrix.c File Reference

Invert matrix through QR decomposition. More...

Detailed Description

Invert matrix through QR decomposition.

Author
Pascal Getreuer getre.nosp@m.uer@.nosp@m.cmla..nosp@m.ens-.nosp@m.cacha.nosp@m.n.fr

Copyright (c) 2012-2013, Pascal Getreuer All rights reserved.

This program is free software: you can redistribute it and/or modify it under, at your option, the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version, or the terms of the simplified BSD license.

You should have received a copy of these licenses along with this program. If not, see http://www.gnu.org/licenses/ and http://www.opensource.org/licenses/bsd-license.html.

Definition in file invert_matrix.c.

#include "invert_matrix.h"
#include <assert.h>
#include <math.h>
#include <stdlib.h>

Go to the source code of this file.

Functions

int invert_matrix (double *inv_A, double *A, int N)
 Invert matrix through QR decomposition. More...
 

Function Documentation

int invert_matrix ( double *  inv_A,
double *  A,
int  N 
)

Invert matrix through QR decomposition.

Parameters
inv_Apointer to memory for holding the result
Apointer to column-major matrix data
Nthe number of dimensions
Returns
1 on success, 0 on failure

The input data is overwritten during the computation. inv_A should be allocated before calling this function with space for at least N^2 doubles. Matrices are represented in column-major format, meaning A(i,j) = A[i + N*j], 0 <= i, j < N.

Definition at line 37 of file invert_matrix.c.