Index: uspace/lib/math/arch/amd64/src/cos.S
===================================================================
--- uspace/lib/math/arch/amd64/src/cos.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/amd64/src/cos.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -35,31 +35,31 @@
 	pushq %rbp
 	movq %rsp, %rbp
-	
+
 	# compute cosine (using red zone)
-	
+
 	movsd %xmm0, -8(%rbp)
 	fldl -8(%rbp)
-	
+
 	fcos
-	
+
 	# detect if source operand is out of range
-	
+
 	fnstsw %ax
 	andw $X87_STATUS_WORD_C2_MASK, %ax
 	jnz fix_range
-	
+
 		fstpl -8(%rbp)
 		movsd -8(%rbp), %xmm0
-		
+
 		leave
 		retq
-	
+
 	# argument reduction
-	
+
 	fix_range:
 		fldpi
 		fadd %st(0)
 		fxch %st(1)
-	
+
 	reduce:
 		fprem1
@@ -67,11 +67,11 @@
 		andw $X87_STATUS_WORD_C2_MASK, %ax
 		jnz reduce
-	
+
 	fstp %st(1)
 	fcos
-	
+
 	fstpl -8(%rbp)
 	movsd -8(%rbp), %xmm0
-	
+
 	leave
 	retq
Index: uspace/lib/math/arch/amd64/src/sin.S
===================================================================
--- uspace/lib/math/arch/amd64/src/sin.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/amd64/src/sin.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -35,31 +35,31 @@
 	pushq %rbp
 	movq %rsp, %rbp
-	
+
 	# compute sine (using red zone)
-	
+
 	movsd %xmm0, -8(%rbp)
 	fldl -8(%rbp)
-	
+
 	fsin
-	
+
 	# detect if source operand is out of range
-	
+
 	fnstsw %ax
 	andw $X87_STATUS_WORD_C2_MASK, %ax
 	jnz fix_range
-	
+
 		fstpl -8(%rbp)
 		movsd -8(%rbp), %xmm0
-		
+
 		leave
 		retq
-	
+
 	# argument reduction
-	
+
 	fix_range:
 		fldpi
 		fadd %st(0)
 		fxch %st(1)
-	
+
 	reduce:
 		fprem1
@@ -67,11 +67,11 @@
 		andw $X87_STATUS_WORD_C2_MASK, %ax
 		jnz reduce
-	
+
 	fstp %st(1)
 	fsin
-	
+
 	fstpl -8(%rbp)
 	movsd -8(%rbp), %xmm0
-	
+
 	leave
 	retq
Index: uspace/lib/math/arch/amd64/src/trunc.S
===================================================================
--- uspace/lib/math/arch/amd64/src/trunc.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/amd64/src/trunc.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -35,19 +35,19 @@
 	pushq %rbp
 	movq %rsp, %rbp
-	
+
 	# store x87 control word in the red zone
-	
+
 	fnstcw -8(%rbp)
 	movw -8(%rbp), %ax
-	
+
 	# set rounding control to truncate
 	# (no masking necessary for this flag)
-	
+
 	orw $X87_CONTROL_WORD_RC_TRUNCATE, %ax
 	movw %ax, -16(%rbp)
 	fldcw -16(%rbp)
-	
+
 	# truncate
-	
+
 	movsd %xmm0, -16(%rbp)
 	fldl -16(%rbp)
@@ -55,9 +55,9 @@
 	fstpl -16(%rbp)
 	movsd -16(%rbp), %xmm0
-	
+
 	# restore original control word
-	
+
 	fldcw -8(%rbp)
-	
+
 	leave
 	retq
Index: uspace/lib/math/arch/ia32/src/cos.S
===================================================================
--- uspace/lib/math/arch/ia32/src/cos.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/ia32/src/cos.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -34,23 +34,23 @@
 FUNCTION_BEGIN(cos_f64)
 	# compute cosine (no stack frame)
-	
+
 	fldl 4(%esp)
 	fcos
-	
+
 	# detect if source operand is out of range
-	
+
 	fnstsw %ax
 	andw $X87_STATUS_WORD_C2_MASK, %ax
 	jnz fix_range
-	
+
 		ret
-	
+
 	# argument reduction
-	
+
 	fix_range:
 		fldpi
 		fadd %st(0)
 		fxch %st(1)
-	
+
 	reduce:
 		fprem1
@@ -58,8 +58,8 @@
 		andw $X87_STATUS_WORD_C2_MASK, %ax
 		jnz reduce
-	
+
 	fstp %st(1)
 	fcos
-	
+
 	ret
 FUNCTION_END(cos_f64)
Index: uspace/lib/math/arch/ia32/src/sin.S
===================================================================
--- uspace/lib/math/arch/ia32/src/sin.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/ia32/src/sin.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -34,23 +34,23 @@
 FUNCTION_BEGIN(sin_f64)
 	# compute sine (no stack frame)
-	
+
 	fldl 4(%esp)
 	fsin
-	
+
 	# detect if source operand is out of range
-	
+
 	fnstsw %ax
 	andw $X87_STATUS_WORD_C2_MASK, %ax
 	jnz fix_range
-	
+
 		ret
-	
+
 	# argument reduction
-	
+
 	fix_range:
 		fldpi
 		fadd %st(0)
 		fxch %st(1)
-	
+
 	reduce:
 		fprem1
@@ -58,8 +58,8 @@
 		andw $X87_STATUS_WORD_C2_MASK, %ax
 		jnz reduce
-	
+
 	fstp %st(1)
 	fsin
-	
+
 	ret
 FUNCTION_END(sin_f64)
Index: uspace/lib/math/arch/ia32/src/trunc.S
===================================================================
--- uspace/lib/math/arch/ia32/src/trunc.S	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/arch/ia32/src/trunc.S	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -36,26 +36,26 @@
 	movl %esp, %ebp
 	subl $8, %esp
-	
+
 	# store x87 control word
-	
+
 	fnstcw -4(%ebp)
 	movw -4(%ebp), %ax
-	
+
 	# set rounding control to truncate
 	# (no masking necessary for this flag)
-	
+
 	orw $X87_CONTROL_WORD_RC_TRUNCATE, %ax
 	movw %ax, -8(%ebp)
 	fldcw -8(%ebp)
-	
+
 	# truncate
-	
+
 	fldl 8(%ebp)
 	frndint
-	
+
 	# restore original control word
-	
+
 	fldcw -4(%ebp)
-	
+
 	leave
 	ret
Index: uspace/lib/math/generic/asin.c
===================================================================
--- uspace/lib/math/generic/asin.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/asin.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -49,10 +49,10 @@
 {
 	float32_t aval;
-	
+
 	if (arg < -1.0 || arg > 1.0) {
 		errno = EDOM;
 		return FLOAT32_NAN;
 	}
-	
+
 	aval = 2.0 * atan_f32(arg / (1.0 + sqrt_f32(1.0 - arg*arg)));
 	if (arg > 0.0)
@@ -74,10 +74,10 @@
 {
 	float64_t aval;
-	
+
 	if (arg < -1.0 || arg > 1.0) {
 		errno = EDOM;
 		return FLOAT64_NAN;
 	}
-	
+
 	aval = 2.0 * atan_f64(arg / (1.0 + sqrt_f64(1.0 - arg*arg)));
 	if (arg > 0.0)
Index: uspace/lib/math/generic/atan.c
===================================================================
--- uspace/lib/math/generic/atan.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/atan.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -54,10 +54,10 @@
 	float32_t sum = 0;
 	float32_t a = arg / (1.0 + arg * arg);
-	
+
 	/*
 	 * atan(z) = sum(n=0, +inf) [ (2^2n) * (n!)^2 / (2n + 1)! *
 	 *    z^(2n+1) / (1 + z^2)^(n+1) ]
 	 */
-	
+
 	for (unsigned int n = 0; n < SERIES_DEGREE_32; n++) {
 		if (n > 0) {
@@ -68,5 +68,5 @@
 		a = a * 4.0 * arg * arg / (1.0 + arg * arg);
 	}
-	
+
 	return sum;
 }
@@ -86,10 +86,10 @@
 	float64_t sum = 0;
 	float64_t a = arg / (1.0 + arg * arg);
-	
+
 	/*
 	 * atan(z) = sum(n=0, +inf) [ (2^2n) * (n!)^2 / (2n + 1)! *
 	 *    z^(2n+1) / (1 + z^2)^(n+1) ]
 	 */
-	
+
 	for (unsigned int n = 0; n < SERIES_DEGREE_64; n++) {
 		if (n > 0) {
@@ -100,5 +100,5 @@
 		a = a * 4.0 * arg * arg / (1.0 + arg * arg);
 	}
-	
+
 	return sum;
 }
Index: uspace/lib/math/generic/ceil.c
===================================================================
--- uspace/lib/math/generic/ceil.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/ceil.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -48,8 +48,8 @@
 	float32_u v;
 	float32_u r;
-	
+
 	v.val = val;
 	t.val = trunc_f32(val);
-	
+
 	if (v.data.parts.sign == 1 || val == t.val) {
 		r = t;
@@ -57,5 +57,5 @@
 		r.val = t.val + 1.0;
 	}
-	
+
 	return r.val;
 }
@@ -72,8 +72,8 @@
 	float64_u v;
 	float64_u r;
-	
+
 	v.val = val;
 	t.val = trunc_f64(val);
-	
+
 	if (v.data.parts.sign == 1 || val == t.val) {
 		r = t;
@@ -81,5 +81,5 @@
 		r.val = t.val + 1.0;
 	}
-	
+
 	return r.val;
 }
Index: uspace/lib/math/generic/exp.c
===================================================================
--- uspace/lib/math/generic/exp.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/exp.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -64,10 +64,10 @@
 	float32_t ret = 1;
 	float32_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_32; i++) {
 		nom *= arg;
 		ret += nom / factorials[i];
 	}
-	
+
 	return ret;
 }
@@ -89,10 +89,10 @@
 	float64_t ret = 1;
 	float64_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_64; i++) {
 		nom *= arg;
 		ret += nom / factorials[i];
 	}
-	
+
 	return ret;
 }
Index: uspace/lib/math/generic/floor.c
===================================================================
--- uspace/lib/math/generic/floor.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/floor.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -48,8 +48,8 @@
 	float32_t t;
 	float32_u v;
-	
+
 	v.val = val;
 	t = trunc_f32(val);
-	
+
 	if (v.data.parts.sign == 0 || val == t)
 		return t;
@@ -69,8 +69,8 @@
 	float64_t t;
 	float64_u v;
-	
+
 	v.val = val;
 	t = trunc_f64(val);
-	
+
 	if (v.data.parts.sign == 0 || val == t)
 		return t;
Index: uspace/lib/math/generic/fmod.c
===================================================================
--- uspace/lib/math/generic/fmod.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/fmod.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -55,7 +55,7 @@
 {
 	// FIXME: replace with exact arithmetics
-	
+
 	float32_t quotient = trunc_f32(dividend / divisor);
-	
+
 	return (dividend - quotient * divisor);
 }
@@ -80,7 +80,7 @@
 {
 	// FIXME: replace with exact arithmetics
-	
+
 	float64_t quotient = trunc_f64(dividend / divisor);
-	
+
 	return (dividend - quotient * divisor);
 }
Index: uspace/lib/math/generic/log.c
===================================================================
--- uspace/lib/math/generic/log.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/log.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -55,8 +55,8 @@
 	float32_t ret = 0;
 	float32_t num = 1;
-	
+
 	for (unsigned int i = 1; i <= TAYLOR_DEGREE_32; i++) {
 		num *= arg;
-		
+
 		if ((i % 2) == 0)
 			ret += num / i;
@@ -64,5 +64,5 @@
 			ret -= num / i;
 	}
-	
+
 	return ret;
 }
@@ -83,8 +83,8 @@
 	float64_t ret = 0;
 	float64_t num = 1;
-	
+
 	for (unsigned int i = 1; i <= TAYLOR_DEGREE_64; i++) {
 		num *= arg;
-		
+
 		if ((i % 2) == 0)
 			ret += num / i;
@@ -92,5 +92,5 @@
 			ret -= num / i;
 	}
-	
+
 	return ret;
 }
@@ -107,5 +107,5 @@
 	float32_u m;
 	int e;
-	
+
 	m.val = arg;
 	/*
@@ -118,5 +118,5 @@
 	e = m.data.parts.exp - (FLOAT32_BIAS - 1);
 	m.data.parts.exp = FLOAT32_BIAS - 1;
-	
+
 	/*
 	 * arg = m * 2^e ; log(arg) = log(m) + log(2^e) =
@@ -137,7 +137,7 @@
 	float64_u m;
 	int e;
-	
+
 	m.val = arg;
-	
+
 	/*
 	 * Factor arg into m * 2^e where m has exponent -1,
@@ -149,5 +149,5 @@
 	e = m.data.parts.exp - (FLOAT64_BIAS - 1);
 	m.data.parts.exp = FLOAT64_BIAS - 1;
-	
+
 	/*
 	 * arg = m * 2^e ; log(arg) = log(m) + log(2^e) =
Index: uspace/lib/math/generic/trig.c
===================================================================
--- uspace/lib/math/generic/trig.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/trig.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -64,8 +64,8 @@
 	float32_t ret = 0;
 	float32_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_32; i++) {
 		nom *= arg;
-		
+
 		if ((i % 4) == 0)
 			ret += nom / factorials[i];
@@ -73,5 +73,5 @@
 			ret -= nom / factorials[i];
 	}
-	
+
 	return ret;
 }
@@ -93,8 +93,8 @@
 	float64_t ret = 0;
 	float64_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_64; i++) {
 		nom *= arg;
-		
+
 		if ((i % 4) == 0)
 			ret += nom / factorials[i];
@@ -102,5 +102,5 @@
 			ret -= nom / factorials[i];
 	}
-	
+
 	return ret;
 }
@@ -122,8 +122,8 @@
 	float32_t ret = 1;
 	float32_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_32; i++) {
 		nom *= arg;
-		
+
 		if ((i % 4) == 1)
 			ret -= nom / factorials[i];
@@ -131,5 +131,5 @@
 			ret += nom / factorials[i];
 	}
-	
+
 	return ret;
 }
@@ -151,8 +151,8 @@
 	float64_t ret = 1;
 	float64_t nom = 1;
-	
+
 	for (unsigned int i = 0; i < TAYLOR_DEGREE_64; i++) {
 		nom *= arg;
-		
+
 		if ((i % 4) == 1)
 			ret -= nom / factorials[i];
@@ -160,5 +160,5 @@
 			ret += nom / factorials[i];
 	}
-	
+
 	return ret;
 }
@@ -179,5 +179,5 @@
 {
 	unsigned int period = arg / (M_PI / 4);
-	
+
 	switch (period) {
 	case 0:
@@ -212,5 +212,5 @@
 {
 	unsigned int period = arg / (M_PI / 4);
-	
+
 	switch (period) {
 	case 0:
@@ -245,5 +245,5 @@
 {
 	unsigned int period = arg / (M_PI / 4);
-	
+
 	switch (period) {
 	case 0:
@@ -278,5 +278,5 @@
 {
 	unsigned int period = arg / (M_PI / 4);
-	
+
 	switch (period) {
 	case 0:
@@ -308,8 +308,8 @@
 {
 	float32_t base_arg = fmod_f32(arg, 2 * M_PI);
-	
+
 	if (base_arg < 0)
 		return -base_sin_32(-base_arg);
-	
+
 	return base_sin_32(base_arg);
 }
@@ -327,8 +327,8 @@
 {
 	float64_t base_arg = fmod_f64(arg, 2 * M_PI);
-	
+
 	if (base_arg < 0)
 		return -base_sin_64(-base_arg);
-	
+
 	return base_sin_64(base_arg);
 }
@@ -346,8 +346,8 @@
 {
 	float32_t base_arg = fmod_f32(arg, 2 * M_PI);
-	
+
 	if (base_arg < 0)
 		return base_cos_32(-base_arg);
-	
+
 	return base_cos_32(base_arg);
 }
@@ -365,8 +365,8 @@
 {
 	float64_t base_arg = fmod_f64(arg, 2 * M_PI);
-	
+
 	if (base_arg < 0)
 		return base_cos_64(-base_arg);
-	
+
 	return base_cos_64(base_arg);
 }
Index: uspace/lib/math/generic/trunc.c
===================================================================
--- uspace/lib/math/generic/trunc.c	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/generic/trunc.c	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -56,8 +56,8 @@
 	float32_u v;
 	int32_t exp;
-	
+
 	v.val = val;
 	exp = v.data.parts.exp - FLOAT32_BIAS;
-	
+
 	if (exp < 0) {
 		/* -1 < val < 1 => result is +0 or -0 */
@@ -69,5 +69,5 @@
 			// FIXME TODO
 		}
-		
+
 		/* All bits in val are relevant for the result */
 	} else {
@@ -75,5 +75,5 @@
 		v.data.parts.fraction &= ~(UINT32_C(0x007fffff) >> exp);
 	}
-	
+
 	return v.val;
 }
@@ -98,8 +98,8 @@
 	float64_u v;
 	int32_t exp;
-	
+
 	v.val = val;
 	exp = v.data.parts.exp - FLOAT64_BIAS;
-	
+
 	if (exp < 0) {
 		/* -1 < val < 1 => result is +0 or -0 */
@@ -111,5 +111,5 @@
 			// FIXME TODO
 		}
-		
+
 		/* All bits in val are relevant for the result */
 	} else {
@@ -117,5 +117,5 @@
 		v.data.parts.fraction &= ~(UINT64_C(0x000fffffffffffff) >> exp);
 	}
-	
+
 	return v.val;
 }
Index: uspace/lib/math/include/mathtypes.h
===================================================================
--- uspace/lib/math/include/mathtypes.h	(revision 9d8b12da8c497c2a0001860cd8545439328c6071)
+++ uspace/lib/math/include/mathtypes.h	(revision 615e83d68dc1cfdb51a3e2557d4776a37010a7ca)
@@ -91,5 +91,5 @@
 typedef union {
 	uint32_t bin;
-	
+
 	struct {
 		uint32_t sign : 1;
@@ -101,5 +101,5 @@
 typedef union {
 	uint64_t bin;
-	
+
 	struct {
 		uint64_t sign : 1;
@@ -114,5 +114,5 @@
 		uint32_t lo;
 	} bin __attribute__((packed));
-	
+
 	struct {
 		uint64_t padding : 16;
@@ -128,5 +128,5 @@
 		uint64_t lo;
 	} bin __attribute__((packed));
-	
+
 	struct {
 		uint64_t sign : 1;
@@ -141,5 +141,5 @@
 typedef union {
 	uint32_t bin;
-	
+
 	struct {
 		uint32_t fraction : 23;
@@ -151,5 +151,5 @@
 typedef union {
 	uint64_t bin;
-	
+
 	struct {
 		uint64_t fraction : 52;
@@ -164,5 +164,5 @@
 		uint64_t hi;
 	} bin __attribute__((packed));
-	
+
 	struct {
 		uint64_t fraction : 64;
@@ -178,5 +178,5 @@
 		uint64_t hi;
 	} bin __attribute__((packed));
-	
+
 	struct {
 		uint64_t frac_lo : 64;