Java code:

	public static double roundHalfUp(double value) {
		// parts with CAREFULLY can compile to wrong code in broken compilers
		
		// next integer after largest positive fractional double, even
		double allIntsEven = 4503599627370496.0; // 2^52
		// next integer after largest positive fractional double, but odd
		double allIntsOdd = 4503599627370497.0; // 2^52 + 1
		
		// CAREFULLY, turn signed zeroes and other values to infinities of the
		// same sign, by using 1/0 => INF and 1/(-0) => -INF
		double negative = 1.0 / (0.0 * value) < 0.0 ? -1.0 : 1.0;
		
		// if needed flip the sign to work with positive values only
		double absValue = value * negative;
		
		// avoid round-off with large values which are already integers
		double use = absValue < allIntsEven ? 1.0 : 0.0;
		
		// CAREFULLY, induce round-off in an intermediate value of the same
		// evenness as input, ties broken towards even values
		// CAREFULLY, flip evenness of intermediate value and repeat
		double even = use * -allIntsEven + (use * allIntsEven + absValue);
		double odd = use * -allIntsOdd + (use * allIntsOdd + absValue);
		
		// choose larger of tie-breaks of unknown evenness
		double up = even > odd ? even : odd;
		
		// CAREFULLY, add back sign or use special values as is, by using
		// INF-INF <=> -INF-(-INF) <=> NaN-NaN => NaN != 0.0
		double signed = value - value == 0.0 ? up * negative : value;
		return signed;
	}
	public static void main(String[] args) {
		double[][] tests = {
			{ -0.49999999999999994, -0.0 },
			{ -0.5, -1.0 },
			{ -0.5000000000000001, -1.0 },
			{ 0.0, 0.0 },
			{ 0.49999999999999994, 0.0 },
			{ 0.5, 1.0 },
			{ 0.5000000000000001, 1.0 },
			{ 1.390671161567e-309, 0.0 },
			{ 2.2517998136852485e+15, 2.251799813685249e+15 },
			{ 4.503599627370497e+15, 4.503599627370497e+15 },
			{ Double.NEGATIVE_INFINITY, Double.NEGATIVE_INFINITY },
			{ Double.POSITIVE_INFINITY, Double.POSITIVE_INFINITY },
			{ Double.NaN, Double.NaN },
			{ -0.0, -0.0 },
		};
		for(double[] test : tests) {
			double got = roundHalfUp(test[0]);
			boolean ok = Double.doubleToLongBits(got) == Double.doubleToLongBits(test[1]);
			System.out.printf("%22s ~%s %22s %s= %22s\n", test[0], ok ? '=' : '?', got, ok ? '=' : '!', test[1]);
		}
	}