cov_scale.c

最大似然估计算法

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

void cov_scale(double *C_unit, int nrow, double *params, double *C) {

	int j, k;

	double pl2, wh2;

	pl2 = params[1] * params[1];
	wh2 = params[0] * params[0];

	for (j = 0; j < nrow; j++) {
		for (k = 0; k < nrow; k++) {
			C[j + k * nrow] = C_unit[j + k * nrow] * pl2;
			if (j == k) C[j + k * nrow] += wh2;
		}
	}
}