Example #1
0
static void randomizeMatrix(Matrix3d &matrix) {
	ssize_t size[3] = { (ssize_t)matrix.size(0), (ssize_t)matrix.size(1), (ssize_t)matrix.size(2) };
	const ssize_t rim = matrix.rim();
	for(ssize_t ix=-rim;ix<size[0]+rim;ix++) {
		for(ssize_t iy=-rim;iy<size[1]+rim;iy++) {
			for(ssize_t iz=-rim;iz<size[2]+rim;iz++) {
				matrix(ix,iy,iz) = matrix_func(ix,iy,iz); //(ix*ix) + ix*iy + iz;
			}
		}
	}
}
Example #2
0
static void diagonalizeMatrix(Matrix3d &matrix) {
	ssize_t size[3] = { (ssize_t)matrix.size(0), (ssize_t)matrix.size(1), (ssize_t)matrix.size(2) };
		const ssize_t rim = matrix.rim();
		for(ssize_t ix=-rim;ix<size[0]+rim;ix++) {
			for(ssize_t iy=-rim;iy<size[1]+rim;iy++) {
				for(ssize_t iz=-rim;iz<size[2]+rim;iz++) {
					bool isDiagonal = ix==iy;// && iy==iz;
					matrix(ix,iy,iz) = isDiagonal?1.0:0.0;
				}
			}
	}
}