crecord.cpp

LCL分群法 改善cure kmean.......等

#include "stdafx.h"
#include <xmmintrin.h>
#include <math.h>
#include "CRecord.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

IMPLEMENT_SERIAL(CRecord,CObject,VERSIONABLE_SCHEMA | 2)

CRecord::CRecord()
:mpVector(NULL),mlDim(0)
{
}

CRecord::CRecord(const char *label,float *vector,long dim)
:mpVector(NULL),mlDim(0)
{
	moLabel = label;
	init(vector,dim);
}

CRecord::~CRecord()
{
	if(mpVector) { _aligned_free(mpVector); }
}

void CRecord::init(float *vector,long dim)
{
	alloc(dim);
	if(vector){
		memcpy(mpVector,vector,sizeof(float)*dim);
	} else {
		memset(mpVector,0,sizeof(float)*dim);
	}
}

void CRecord::alloc(long dim)
{
	long mod = dim % 4;
	long count = dim;
	if(mod != 0){
		count += (4 - mod);
	}
	if(mpVector != NULL) {
		mpVector = (float *)_aligned_realloc(mpVector,sizeof(float)*count,16);
	} else {
		mpVector = (float *)_aligned_malloc(sizeof(float)*count,16);
	}
	if(mod != 0){
		memset(mpVector+mlDim,0,sizeof(float)*(4 - mod));
	}
	mlDim = dim;
}

float CRecord::euclidean_distance(CRecord *record,BOOL sqrt)
{
	float dis,distance=0.0f;
	for(long i=0;i<mlDim;i++){
		dis = mpVector[i] - record->mpVector[i];
		distance += pow(dis,2);
	}
	if(sqrt){
		return ::sqrt(distance);
	}
	return distance;
/*
	__declspec (align(16)) float value[4] = {0.0f,0.0f,0.0f,0.0f};
	__m128 distance = _mm_load_ps(value);
	__m128 dis;
	for(long i=0;i<mlDim;i+=4){
		dis = _mm_sub_ps(*((__m128 *)(mpVector+i)),*((__m128 *)(record->mpVector+i)));
		dis = _mm_mul_ps(dis,dis);
		distance = _mm_add_ps(distance,dis);
	}
	_mm_store_ps(value,distance);
	if(sqrt){
		return ::sqrt(value[0] + value[1] + value[2] + value[3]);
	}
	return value[0] + value[1] + value[2] + value[3];
*/
}