#ifndef __LSSYSTEM_H__
#define __LSSYSTEM_H__

#include <models.h>

void computeLeastSquaresMatrix(double* A,phi_function* phi, int basis_size, r3_t* X,int n, double w[], double wpos[],void* l_data);

void computeLeastSquaresRHSVector(double* b,phi_function* phi, int basis_size,r3_t* X,double* F,int n , double w[], double wpos[],void* l_data);

void FitOneLight(double* A,
	double* b ,
	double* c,
	phi_function* phi,
	int basis_size,
	validate_function* validateResults,
	r3_t* X,
	double* F,
	int n,
	double* w,
	double* wpos,
	void* l_data
	);
	
void FitGaugeToModel(r3_t* X, double* F, int n,
		     lighting_model_t* lm,
		     r3_t view_dir,
		     void* l_data
		     );
#endif