lambda.c File Reference

#include <string.h>
#include <math.h>
#include <stdio.h>
#include <cblas.h>
#include <clapack.h>
#include "amb_kf.h"

Include dependency graph for lambda.c:

Go to the source code of this file.

Defines
#define	LOOPMAX   10000 /* maximum count of search loop */
#define	ROUND(x)   (floor((x)+0.5))
#define	SGN(x)   ((x)<=0.0?-1.0:1.0)
#define	SWAP(x, y)   do {double tmp_; tmp_=x; x=y; y=tmp_;} while (0)
Functions
static void	eye2 (u32 n, double *M)
void	gauss (int n, double *L, double *Z, int i, int j)
int	lambda_reduction (int n, const double *Q, double *Z)
int	lambda_solution (int n, int m, const double *a, const double *Q, double *F, double *s)
int	LD (int n, const double *Q, double *L, double *D)
void	matmul (const char *tr, integer n, integer k, integer m, double alpha, const double *A, const double *B, double beta, double *C)
void	perm (int n, double *L, double *D, int j, double del, double *Z)
void	reduction (int n, double *L, double *D, double *Z)
static int	search (int n, int m, const double *L, const double *D, const double *zs, double *zn, double *s)
int	solve (const char *tr, const double *A, const double *Y, integer n, integer m, double *X)
static void	triu2 (u32 n, double *M)
static s8	udu2 (u32 n, double *M, double *U, double *D)

---

Define Documentation

#define LOOPMAX   10000 /* maximum count of search loop */

Definition at line 26 of file lambda.c.

#define ROUND

(

x

)

(floor((x)+0.5))

Definition at line 29 of file lambda.c.

#define SGN

(

x

)

((x)<=0.0?-1.0:1.0)

Definition at line 28 of file lambda.c.

#define SWAP	(		x,
			y
	)		do {double tmp_; tmp_=x; x=y; y=tmp_;} while (0)

Definition at line 30 of file lambda.c.

---

Function Documentation

static void eye2	(	u32	n,
		double *	M
	)		[static]

Definition at line 42 of file lambda.c.

void gauss	(	int	n,
		double *	L,
		double *	Z,
		int	i,
		int	j
	)

Definition at line 104 of file lambda.c.

int lambda_reduction	(	int	n,
		const double *	Q,
		double *	Z
	)

Definition at line 223 of file lambda.c.

int lambda_solution	(	int	n,
		int	m,
		const double *	a,
		const double *	Q,
		double *	F,
		double *	s
	)

Definition at line 306 of file lambda.c.

int LD	(	int	n,
		const double *	Q,
		double *	L,
		double *	D
	)

Definition at line 79 of file lambda.c.

void matmul	(	const char *	tr,
		integer	n,
		integer	k,
		integer	m,
		double	alpha,
		const double *	A,
		const double *	B,
		double	beta,
		double *	C
	)

Definition at line 259 of file lambda.c.

void perm	(	int	n,
		double *	L,
		double *	D,
		int	j,
		double	del,
		double *	Z
	)

Definition at line 114 of file lambda.c.

void reduction	(	int	n,
		double *	L,
		double *	D,
		double *	Z
	)

Definition at line 132 of file lambda.c.

static int search	(	int	n,
		int	m,
		const double *	L,
		const double *	D,
		const double *	zs,
		double *	zn,
		double *	s
	)		[static]

Definition at line 149 of file lambda.c.

int solve	(	const char *	tr,
		const double *	A,
		const double *	Y,
		integer	n,
		integer	m,
		double *	X
	)

Definition at line 279 of file lambda.c.

static void triu2	(	u32	n,
		double *	M
	)		[static]

Definition at line 33 of file lambda.c.

static s8 udu2	(	u32	n,
		double *	M,
		double *	U,
		double *	D
	)		[static]

Definition at line 50 of file lambda.c.