/*
?*?cgf.c:	conj.?grad.?routine?for?finding?optimal?s?-?fast!
?*/
#include?
#include?
#include?“mex.h“

extern?float?frprmn（）;


#define?sgn（x）?（x>0???1?:?（x<0???-1?:?0））

extern??__declspec（?dllexport?）
void?cgf（double?*Sout?double?*nits?double?*nf?double?*ng
	???????double?*Sin?double?*X?int?npats?
	???????double?tol?int?maxiter?int?numflag）;

/*?Input?&?Output?Arguments?*/

#define	A_IN		prhs[0]		/*?basis?matrix?*/
#define	X_IN		prhs[1]		/*?data?vectors?*/
#define	S_IN		prhs[2]		/*?initial?guess?for?S?*/
#define	LAMBDA_IN	prhs[3]		/*?precision?*/
#define?BETA_IN		prhs[4]		/*?prior?steepness?*/
#define?SIGMA_IN????????prhs[5]?????????/*?scaling?parameter?for?prior?*/
#define	TOL_IN		prhs[6]		/*?tolerance?*/
#define?MAXITER_IN??????prhs[7]		/*?maximum?iterations?for?dfrpmin?*/
#define	OUTFLAG_IN	prhs[8]		/*?output?flag?*/
#define	NUMFLAG_IN	prhs[9]		/*?pattern?number?output?flag?*/

#define	S_OUT???????????plhs[0]		/*?basis?coeffs?for?each?data?vector?*/
#define?NITS_OUT????????plhs[1]?????????/*?total?iterations?done?by?cg?*/
#define?NF_OUT??????????plhs[2]?????????/*?total?P（s|xA）?calcs?*/
#define?NG_OUT??????????plhs[3]?????????/*?total?d/ds?P（s|xA）?calcs?*/

/*?Define?indexing?macros?for?matricies?*/

/*?L?=?dimension?of?input?vectors
?*?M?=?number?of?basis?functions
?*/

#define?A_（ij）		A[（i）?+?（j）*L]		/*?A?is?L?x?M?*/
#define?X_（in）		X[（i）?+?（n）*L]		/*?X?is?L?x?npats?*/

#define?Sout_（in）	Sout[（i）?+?（n）*M]	/*?S?is?M?x?npats?*/
#define?Sin_（in）	Sin[（i）?+?（n）*M]	/*?S?is?M?x?npats?*/

#define?AtA_（ij）	AtA[（i）?+?（j）*M]	/*?AtA?is?M?x?M?*/

/*?Globals?for?using?with?frprmin?*/

static?double?*A;		/*?basis?matrix?*/
static?int?L;			/*?data?dimension?*/
static?int?M;			/*?number?of?basis?vectors?*/
static?double?lambda;		/*?1/noise_var?*/
static?double?beta;		/*?prior?steepness?*/
static?double?sigma;		/*?prior?scaling?*/
static?double?k1k2k3;		/*?precomputed?constants?for?f1dim?*/

static?double?*x;		/*?current?data?vector?being?fitted?*/
static?double?*s0;		/*?init?coefficient?vector?（1:M）?*/
static?double?*d;		/*?search?dir.?coefficient?vector?（1:M）?*/
static?int?outflag;		/*?print?search?progress?*/

static?double?*AtA;		/*?Only?compute?A‘*A?once?（1:M1:M）?*/
static?double?*Atx;		/*?A*x?（1:M）?*/

static?int?fcount?gcount;

static?void?init_global_arrays（）?
{
??int?		ijk;
??double	*Ai?*Aj?sum;

??x??????=?（double?*）malloc（L*sizeof（double））;
??s0?????=?（double?*）malloc（M*sizeof（double））;
??d??????=?（double?*）malloc（M*sizeof（double））;
??AtA????=?（double?*）malloc（M*M*sizeof（double））;
??Atx????=?（double?*）malloc（M*sizeof（double））;

??/*?Calc??A‘*A?*/
??for?（i?=?0;?i?????Ai=A+i*L;
????for?（j?=?0;?j???????Aj=A+j*L;
??????sum=0.0;
??????for?（k?=?0;?k?	sum?+=?Ai[k]*Aj[k];
??????}
??????AtA_（ij）?=?sum;
????}
??}
}

static?void?free_global_arrays（）?{

??free（（double?*）x）;
??free（（double?*）s0）;
??free（（double?*）d）;
??free（（double?*）AtA）;
??free（（double?*）Atx）;
}



float?init_f1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2002-03-11?13:41??Code\Matlab\Bruno\sparsenet\
?????文件??????524472??2002-02-28?11:52??Code\Matlab\Bruno\sparsenet\A16.mat
?????文件????????7457??2002-02-27?12:22??Code\Matlab\Bruno\sparsenet\cgf.c
?????文件???????40960??2002-02-27?13:50??Code\Matlab\Bruno\sparsenet\cgf.dll
?????文件???????35250??2002-02-27?13:50??Code\Matlab\Bruno\sparsenet\cgf.map
?????文件????????1419??2002-02-27?14:35??Code\Matlab\Bruno\sparsenet\cgf_fitS.m
?????文件?????????442??1998-04-14?21:43??Code\Matlab\Bruno\sparsenet\display_network.m
?????文件????20971712??2000-07-12?12:05??Code\Matlab\Bruno\sparsenet\IMAGES.mat
?????文件?????????975??2000-07-12?12:07??Code\Matlab\Bruno\sparsenet\make-your-own-images
?????文件??????????97??2002-02-27?13:42??Code\Matlab\Bruno\sparsenet\makefile
?????目錄???????????0??2002-03-11?13:26??Code\Matlab\Bruno\sparsenet\nrf\
?????文件????????1494??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\brent.c
?????文件????????2010??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\brent.obj
?????文件????????1204??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\frprmn.c
?????文件????????1603??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\frprmn.obj
?????文件?????????447??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\linmin.c
?????文件?????????664??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\linmin.obj
?????文件?????????386??2002-02-27?13:44??Code\Matlab\Bruno\sparsenet\nrf\makefile
?????文件?????????176??2002-02-26?14:23??Code\Matlab\Bruno\sparsenet\nrf\Makefile.original
?????文件?????????386??2002-02-27?13:46??Code\Matlab\Bruno\sparsenet\nrf\Makefile.vc
?????文件?????????419??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\makefile.vcG6
?????文件?????????509??2002-03-11?13:26??Code\Matlab\Bruno\sparsenet\nrf\makefile.vcGB
?????文件????????1309??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\mnbrak.c
?????文件????????1814??2002-03-11?13:25??Code\Matlab\Bruno\sparsenet\nrf\mnbrak.obj
?????文件???????45056??2002-02-27?13:44??Code\Matlab\Bruno\sparsenet\nrf\nrfopt.dll
?????文件????????1073??2002-02-27?13:44??Code\Matlab\Bruno\sparsenet\nrf\nrfopt.exp
?????文件?????????609??2002-02-27?10:41??Code\Matlab\Bruno\sparsenet\nrf\nrfopt.h
?????文件????????2860??2002-02-27?13:44??Code\Matlab\Bruno\sparsenet\nrf\nrfopt.lib
?????文件???????34854??2002-02-27?13:44??Code\Matlab\Bruno\sparsenet\nrf\nrfopt.map
?????文件???????45056??2002-03-11?13:26??Code\Matlab\Bruno\sparsenet\nrf\nrfoptGB.dll
?????文件????????1081??2002-03-11?13:26??Code\Matlab\Bruno\sparsenet\nrf\nrfoptGB.exp
............此處省略17個文件信息