Head

GCC Code Coverage Report

Directory:	./		Exec	Total	Coverage
File:	usr.bin/calendar/pesach.c	Lines:	0	13	0.0 %
Date:	2016-12-06	Branches:	0	20	0.0 %


/*	$OpenBSD: pesach.c,v 1.4 2015/03/15 00:41:28 millert Exp $	*/

/*
 * Copyright (c) 2004 Michael Shalayeff
 * All rights reserved.
 *
 * Permission to use, copy, modify, and distribute this software for any
 * purpose with or without fee is hereby granted, provided that the above
 * copyright notice and this permission notice appear in all copies.
 *
 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
 * WHATSOEVER RESULTING FROM LOSS OF MIND, USE, DATA OR PROFITS, WHETHER IN
 * AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT
 * OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
 */

#include <stdio.h>

#include "calendar.h"

/* Calculate the Julian date of Pesach using the Gauss formula */

#define	T	(33. + 14. / 24.)
#define	L	((1.  + 485. / 1080.) / 24. / 19.)
#define	K	((29. + (12. + 793. / 1080.) / 24. ) / 19.)

int
pesach(int R)
{
	int a, b, y, cumdays;
	double d;

	y = R + 3760;

	a = (12 * y + 17) % 19;
	b = y % 4;
	d = (T - 10 * K + L + 14) + K * a +  b / 4. - L * y;
	cumdays = d;

	/* the postponement */
	switch ((int)(cumdays + 3 * y + 5 * b + 5) % 7) {
	case 1:
		if (a > 6 && d - cumdays >= (15. + 204. / 1080.) / 24.)
			cumdays += 2;
		break;

	case 0:
		if (a <= 11 || d - cumdays < (21. + 589. / 1080.) / 24.)
			break;
		/* FALLTHROUGH */
	case 2:
	case 4:
	case 6:
		cumdays++;
		break;
	}

	if (R > 1582)
		cumdays += R / 100 - R /400 - 2;

	return (31 + 28 + cumdays + (isleap(R)? 1 : 0));
}


Generated by: GCOVR (Version 3.3)

Line	Branch	Exec	Source
1			/* $OpenBSD: pesach.c,v 1.4 2015/03/15 00:41:28 millert Exp $ */
2
3			/*
4			* Copyright (c) 2004 Michael Shalayeff
5			* All rights reserved.
6			*
7			* Permission to use, copy, modify, and distribute this software for any
8			* purpose with or without fee is hereby granted, provided that the above
9			* copyright notice and this permission notice appear in all copies.
10			*
11			* THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
12			* WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
13			* MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
14			* ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
15			* WHATSOEVER RESULTING FROM LOSS OF MIND, USE, DATA OR PROFITS, WHETHER IN
16			* AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT
17			* OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
18			*/
19
20			#include <stdio.h>
21
22			#include "calendar.h"
23
24			/* Calculate the Julian date of Pesach using the Gauss formula */
25
26			#define T (33. + 14. / 24.)
27			#define L ((1. + 485. / 1080.) / 24. / 19.)
28			#define K ((29. + (12. + 793. / 1080.) / 24. ) / 19.)
29
30			int
31			pesach(int R)
32			{
33			int a, b, y, cumdays;
34			double d;
35
36			y = R + 3760;
37
38			a = (12 * y + 17) % 19;
39			b = y % 4;
40			d = (T - 10 * K + L + 14) + K * a + b / 4. - L * y;
41			cumdays = d;
42
43			/* the postponement */
44			switch ((int)(cumdays + 3 * y + 5 * b + 5) % 7) {
45			case 1:
46			if (a > 6 && d - cumdays >= (15. + 204. / 1080.) / 24.)
47			cumdays += 2;
48			break;
49
50			case 0:
51			if (a <= 11 \|\| d - cumdays < (21. + 589. / 1080.) / 24.)
52			break;
53			/* FALLTHROUGH */
54			case 2:
55			case 4:
56			case 6:
57			cumdays++;
58			break;
59			}
60
61			if (R > 1582)
62			cumdays += R / 100 - R /400 - 2;
63
64			return (31 + 28 + cumdays + (isleap(R)? 1 : 0));
65			}