Svenska ElektronikForumet

Kod: Markera allt

// Repeat, substract, power and then sum
// MATLAB:
// Q = sum(pow(P-repmat(Q, 1, size(P, 2)), 2), 1);
static void repmat_substract_power_sum(float Q[], const float P[]){
	// First row of Repeat -> Substract ->Ppower^2 -> Sum
	float q = *Q; // First column and first row value
	for(uint16_t j = 0; j < Pcolumns; j++){
		*Q = *P - q; // Substract
		*Q = *Q * *Q; // Power^2
		Q++;
		P++;
	}

	// The rest of the rows of Repeat -> Substract -> Power^2 -> Sum
	for(uint16_t i = 1; i < Prows; i++){
		q = *Q; // First column value at row i
		for(uint16_t j = 0; j < Pcolumns; j++){
			*Q = *P - q; // Substract
			Q[j] += *Q * *Q; //  Power^2 -> Sum to the top row at column index j
			Q++;
			P++;
		}
	}
}

Kod: Markera allt

Q[j] += *Q * *Q; //  Power^2 -> Sum to the top row at column index j

Kod: Markera allt

Q[-Pcolumns*i] += *Q * *Q; //  Power^2 -> Sum to the top row at column index j

Kod: Markera allt

Q[-Pcolumns*i] += *Q * *Q;

Kod: Markera allt

float 213*8 
float 65536*8
uint16_t 1*213
float 213*8 
uint16_t 1*213

Kod: Markera allt

float 166*13 
float 77760*13
uint16_t 1*166
float 166*13 
uint16_t 1*166