<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en_US" xml:lang="en_US">
<head>
<title>libkholidays: lunarphase.cpp Source File (libkholidays)</title>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta http-equiv="Content-Style-Type" content="text/css" />
<meta http-equiv="pics-label" content='(pics-1.1 "http://www.icra.org/ratingsv02.html" comment "ICRAonline DE v2.0" l gen true for "http://www.kde.org" r (nz 1 vz 1 lz 1 oz 1 cb 1) "http://www.rsac.org/ratingsv01.html" l gen true for "http://www.kde.org" r (n 0 s 0 v 0 l 0))' />
<meta name="trademark" content="KDE e.V." />
<meta name="description" content="K Desktop Environment Homepage, KDE.org" />
<meta name="MSSmartTagsPreventParsing" content="true" />
<meta name="robots" content="all" />
<link rel="shortcut icon" href="../../favicon.ico" />
<link rel="stylesheet" media="screen" type="text/css" title="APIDOX" href="doxygen.css" />
<style type="text/css">
<!--
hr { display: none; }
#content h2 { margin-left: 0px; }
table.mdTable { background-color: #f8f8f8; border: .2em solid #d7d7d7; }
td.mdRow { padding: 8px 20px; }
td.md { font-weight: bold; }
td.mdname1 { font-weight: bold; color: #602020; }
td.mdname { font-weight: bold; color: #602020; }
-->
</style>
</head>
<body>
<div id="nav_header_top" align="right">
<a href="#content" class="doNotDisplay" accesskey="2">Skip to main content ::</a>
<a href="../.."><img id="nav_header_logo" alt="Home" align="left" src="../../kde_gear_64.png" border="0" /></a>
<span class="doNotDisplay">::</span>
<img id="nav_header_logo_right" alt="" align="right" src="../../pimlogo.png" border="0" />
<div id="nav_header_title" align="left">KDE PIM API Reference</div>
</div>
<div id="nav_header_bottom" align="right">
<span class="doNotDisplay">:: <a href="#navigation" accesskey="5">Skip to Link Menu</a><br/></span>
<div id="nav_header_bottom_left" style="text-align: left;">
/ <a href="../../">API Reference</a>
/ <a href=".">libkholidays</a>
</div>
</div>
<table id="main" border="0" cellpadding="0" cellspacing="0" width="100%">
<tr>
<td valign="top" class="menuheader" height="0"></td>
<td id="contentcolumn" valign="top" rowspan="2" >
<div id="content" style="padding-top: 0px;"><div style="width:100%; margin: 0px; padding: 0px;">
<h2><a name="content"></a>libkholidays</h2>
<!-- Generated by Doxygen 1.5.5 -->
<h1>lunarphase.cpp</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*</span>
<a name="l00002"></a>00002 <span class="comment"> This file is part of libkholidays.</span>
<a name="l00003"></a>00003 <span class="comment"> Copyright (c) 2004 Allen Winter <winter@kde.org></span>
<a name="l00004"></a>00004 <span class="comment"></span>
<a name="l00005"></a>00005 <span class="comment"> Copyright (c) 1989, 1993</span>
<a name="l00006"></a>00006 <span class="comment"> The Regents of the University of California. All rights reserved.</span>
<a name="l00007"></a>00007 <span class="comment"></span>
<a name="l00008"></a>00008 <span class="comment"> This library is free software; you can redistribute it and/or</span>
<a name="l00009"></a>00009 <span class="comment"> modify it under the terms of the GNU General Public License as</span>
<a name="l00010"></a>00010 <span class="comment"> published by the Free Software Foundation; either version 2 of the</span>
<a name="l00011"></a>00011 <span class="comment"> License, or (at your option) any later version.</span>
<a name="l00012"></a>00012 <span class="comment"></span>
<a name="l00013"></a>00013 <span class="comment"> This library is distributed in the hope that it will be useful,</span>
<a name="l00014"></a>00014 <span class="comment"> but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
<a name="l00015"></a>00015 <span class="comment"> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU</span>
<a name="l00016"></a>00016 <span class="comment"> General Public License for more details.</span>
<a name="l00017"></a>00017 <span class="comment"></span>
<a name="l00018"></a>00018 <span class="comment"> You should have received a copy of the GNU General Public License</span>
<a name="l00019"></a>00019 <span class="comment"> along with this program; if not, write to the Free Software</span>
<a name="l00020"></a>00020 <span class="comment"> Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA</span>
<a name="l00021"></a>00021 <span class="comment"></span>
<a name="l00022"></a>00022 <span class="comment"> In addition, as a special exception, the copyright holders give</span>
<a name="l00023"></a>00023 <span class="comment"> permission to link the code of this program with any edition of</span>
<a name="l00024"></a>00024 <span class="comment"> the Qt library by Trolltech AS, Norway (or with modified versions</span>
<a name="l00025"></a>00025 <span class="comment"> of Qt that use the same license as Qt), and distribute linked</span>
<a name="l00026"></a>00026 <span class="comment"> combinations including the two. You must obey the GNU General</span>
<a name="l00027"></a>00027 <span class="comment"> Public License in all respects for all of the code used other than</span>
<a name="l00028"></a>00028 <span class="comment"> Qt. If you modify this file, you may extend this exception to</span>
<a name="l00029"></a>00029 <span class="comment"> your version of the file, but you are not obligated to do so. If</span>
<a name="l00030"></a>00030 <span class="comment"> you do not wish to do so, delete this exception statement from</span>
<a name="l00031"></a>00031 <span class="comment"> your version.</span>
<a name="l00032"></a>00032 <span class="comment">*/</span>
<a name="l00033"></a>00033
<a name="l00034"></a>00034 <span class="preprocessor">#include "config.h"</span>
<a name="l00035"></a>00035
<a name="l00036"></a>00036 <span class="preprocessor">#include <kglobal.h></span>
<a name="l00037"></a>00037 <span class="preprocessor">#include <klocale.h></span>
<a name="l00038"></a>00038 <span class="preprocessor">#include <kdebug.h></span>
<a name="l00039"></a>00039
<a name="l00040"></a>00040 <span class="preprocessor">#include "lunarphase.h"</span>
<a name="l00041"></a>00041
<a name="l00042"></a>00042 LunarPhase::LunarPhase( Hemisphere hemisphere )
<a name="l00043"></a>00043 {
<a name="l00044"></a>00044 mHemisphere = hemisphere;
<a name="l00045"></a>00045 }
<a name="l00046"></a>00046
<a name="l00047"></a>00047 LunarPhase::~LunarPhase()
<a name="l00048"></a>00048 {
<a name="l00049"></a>00049 }
<a name="l00050"></a>00050
<a name="l00051"></a><a class="code" href="classLunarPhase.html#6c4096be9374cc33bb14adb770e817dd">00051</a> <span class="keywordtype">void</span> <a class="code" href="classLunarPhase.html#6c4096be9374cc33bb14adb770e817dd" title="Set the hemisphere.">LunarPhase::setHemisphere</a>( Hemisphere <a class="code" href="classLunarPhase.html#c63caa605e4c5a397d4d9fdb8ca970b4" title="Return the hemisphere.">hemisphere</a> )
<a name="l00052"></a>00052 {
<a name="l00053"></a>00053 mHemisphere = hemisphere;
<a name="l00054"></a>00054 }
<a name="l00055"></a>00055
<a name="l00056"></a><a class="code" href="classLunarPhase.html#c63caa605e4c5a397d4d9fdb8ca970b4">00056</a> LunarPhase::Hemisphere <a class="code" href="classLunarPhase.html#c63caa605e4c5a397d4d9fdb8ca970b4" title="Return the hemisphere.">LunarPhase::hemisphere</a>()<span class="keyword"> const</span>
<a name="l00057"></a>00057 <span class="keyword"></span>{
<a name="l00058"></a>00058 <span class="keywordflow">return</span>( mHemisphere );
<a name="l00059"></a>00059 }
<a name="l00060"></a>00060
<a name="l00061"></a><a class="code" href="classLunarPhase.html#4caa068aac19a40b136e2f119a9a886f">00061</a> QString <a class="code" href="classLunarPhase.html#4caa068aac19a40b136e2f119a9a886f" title="Return hemisphere as a clear text string.">LunarPhase::hemisphereStr</a>()<span class="keyword"> const</span>
<a name="l00062"></a>00062 <span class="keyword"></span>{
<a name="l00063"></a>00063 <span class="keywordflow">return</span> <a class="code" href="classLunarPhase.html#ea28e770ecd2312f5c81ba14f2cc1312" title="Return the string representation of hemisphere.">hemisphereName</a>( mHemisphere );
<a name="l00064"></a>00064 }
<a name="l00065"></a>00065
<a name="l00066"></a><a class="code" href="classLunarPhase.html#ea28e770ecd2312f5c81ba14f2cc1312">00066</a> QString <a class="code" href="classLunarPhase.html#ea28e770ecd2312f5c81ba14f2cc1312" title="Return the string representation of hemisphere.">LunarPhase::hemisphereName</a>( LunarPhase::Hemisphere <a class="code" href="classLunarPhase.html#c63caa605e4c5a397d4d9fdb8ca970b4" title="Return the hemisphere.">hemisphere</a> )
<a name="l00067"></a>00067 {
<a name="l00068"></a>00068 <span class="keywordflow">switch</span>( hemisphere ) {
<a name="l00069"></a>00069 <span class="keywordflow">case</span> Northern:
<a name="l00070"></a>00070 <span class="keywordflow">default</span>:
<a name="l00071"></a>00071 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"Northern"</span> ) );
<a name="l00072"></a>00072 <span class="keywordflow">break</span>;
<a name="l00073"></a>00073 <span class="keywordflow">case</span> Southern:
<a name="l00074"></a>00074 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"Southern"</span> ) );
<a name="l00075"></a>00075 <span class="keywordflow">break</span>;
<a name="l00076"></a>00076 }
<a name="l00077"></a>00077 }
<a name="l00078"></a>00078
<a name="l00079"></a><a class="code" href="classLunarPhase.html#363eacce10dfacbab85f6d7cc52a0fe6">00079</a> QString <a class="code" href="classLunarPhase.html#363eacce10dfacbab85f6d7cc52a0fe6" title="Return the lunar phase as a text string for the specified date.">LunarPhase::phaseStr</a>( <span class="keyword">const</span> QDate &date )<span class="keyword"> const</span>
<a name="l00080"></a>00080 <span class="keyword"></span>{
<a name="l00081"></a>00081 <span class="keywordflow">return</span> <a class="code" href="classLunarPhase.html#e2649d6c1b7138dc1ff82d8feb27eeeb" title="Return the string representation of phase.">phaseName</a>( <a class="code" href="classLunarPhase.html#0000f90da9f3213756589d82dc8ef528" title="Return the lunar phase for the specified Gregorian date.">phase</a>( date ) );
<a name="l00082"></a>00082 }
<a name="l00083"></a>00083
<a name="l00084"></a><a class="code" href="classLunarPhase.html#e2649d6c1b7138dc1ff82d8feb27eeeb">00084</a> QString <a class="code" href="classLunarPhase.html#e2649d6c1b7138dc1ff82d8feb27eeeb" title="Return the string representation of phase.">LunarPhase::phaseName</a>( LunarPhase::Phase <a class="code" href="classLunarPhase.html#0000f90da9f3213756589d82dc8ef528" title="Return the lunar phase for the specified Gregorian date.">phase</a> )
<a name="l00085"></a>00085 {
<a name="l00086"></a>00086 <span class="keywordflow">switch</span> ( phase ) {
<a name="l00087"></a>00087 <span class="keywordflow">case</span> New:
<a name="l00088"></a>00088 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"New Moon"</span> ) );
<a name="l00089"></a>00089 <span class="keywordflow">break</span>;
<a name="l00090"></a>00090 <span class="keywordflow">case</span> Full:
<a name="l00091"></a>00091 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"Full Moon"</span> ) );
<a name="l00092"></a>00092 <span class="keywordflow">break</span>;
<a name="l00093"></a>00093 <span class="keywordflow">case</span> FirstQ:
<a name="l00094"></a>00094 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"First Quarter Moon"</span> ) );
<a name="l00095"></a>00095 <span class="keywordflow">break</span>;
<a name="l00096"></a>00096 <span class="keywordflow">case</span> LastQ:
<a name="l00097"></a>00097 <span class="keywordflow">return</span>( i18n( <span class="stringliteral">"Last Quarter Moon"</span> ) );
<a name="l00098"></a>00098 <span class="keywordflow">break</span>;
<a name="l00099"></a>00099 <span class="keywordflow">default</span>:
<a name="l00100"></a>00100 <span class="keywordflow">case</span> None:
<a name="l00101"></a>00101 <span class="keywordflow">return</span>( QString::null );
<a name="l00102"></a>00102 <span class="keywordflow">break</span>;
<a name="l00103"></a>00103 }
<a name="l00104"></a>00104 }
<a name="l00105"></a>00105
<a name="l00106"></a><a class="code" href="classLunarPhase.html#0000f90da9f3213756589d82dc8ef528">00106</a> LunarPhase::Phase <a class="code" href="classLunarPhase.html#0000f90da9f3213756589d82dc8ef528" title="Return the lunar phase for the specified Gregorian date.">LunarPhase::phase</a>( <span class="keyword">const</span> QDate &date )<span class="keyword"> const</span>
<a name="l00107"></a>00107 <span class="keyword"></span>{
<a name="l00108"></a>00108 Phase retPhase = None;
<a name="l00109"></a>00109
<a name="l00110"></a>00110 <span class="comment">// compute percent-full for the middle of today and yesterday.</span>
<a name="l00111"></a>00111 QTime noontime( 12, 0, 0 );
<a name="l00112"></a>00112 QDateTime today( date, noontime );
<a name="l00113"></a>00113 <span class="keywordtype">double</span> todayPer = percentFull( today.toTime_t() );
<a name="l00114"></a>00114 QDateTime yesterday( date.addDays(-1), noontime );
<a name="l00115"></a>00115 <span class="keywordtype">double</span> yesterdayPer = percentFull( yesterday.toTime_t() );
<a name="l00116"></a>00116
<a name="l00117"></a>00117 <span class="keywordflow">if</span> ( ( todayPer < 0.50 ) && ( yesterdayPer > 0.50 ) ) {
<a name="l00118"></a>00118 retPhase = New;
<a name="l00119"></a>00119 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> ( ( todayPer > 99.50 ) && ( yesterdayPer < 99.50 ) ) {
<a name="l00120"></a>00120 retPhase = Full;
<a name="l00121"></a>00121 } <span class="keywordflow">else</span> {
<a name="l00122"></a>00122 <span class="comment">// compute percent-full for the start of today.</span>
<a name="l00123"></a>00123 QTime sqt( 0, 0, 0 );
<a name="l00124"></a>00124 QDateTime start( date, sqt );
<a name="l00125"></a>00125 <span class="keywordtype">double</span> startPer = percentFull( start.toTime_t() );
<a name="l00126"></a>00126 <span class="comment">// compute percent-full for the end of today.</span>
<a name="l00127"></a>00127 QTime eqt( 23, 59, 59 );
<a name="l00128"></a>00128 QDateTime end( date, eqt );
<a name="l00129"></a>00129 <span class="keywordtype">double</span> endPer = percentFull( end.toTime_t() );
<a name="l00130"></a>00130
<a name="l00131"></a>00131 <span class="keywordflow">if</span> ( ( startPer <= 50 ) && ( endPer > 50 ) ) {
<a name="l00132"></a>00132 <span class="keywordflow">if</span> ( mHemisphere == Northern ) {
<a name="l00133"></a>00133 retPhase = FirstQ;
<a name="l00134"></a>00134 } <span class="keywordflow">else</span> {
<a name="l00135"></a>00135 retPhase = LastQ;
<a name="l00136"></a>00136 }
<a name="l00137"></a>00137 }
<a name="l00138"></a>00138 <span class="keywordflow">if</span> ( ( endPer <= 50 ) && ( startPer > 50 ) ) {
<a name="l00139"></a>00139 <span class="keywordflow">if</span> ( mHemisphere == Northern ) {
<a name="l00140"></a>00140 retPhase = LastQ;
<a name="l00141"></a>00141 } <span class="keywordflow">else</span> {
<a name="l00142"></a>00142 retPhase = FirstQ;
<a name="l00143"></a>00143 }
<a name="l00144"></a>00144 }
<a name="l00145"></a>00145 <span class="comment">// Note: if you want to support crescent and gibbous phases then please</span>
<a name="l00146"></a>00146 <span class="comment">// read the source for the original BSD 'pom' program.</span>
<a name="l00147"></a>00147 }
<a name="l00148"></a>00148 <span class="keywordflow">return</span>( retPhase );
<a name="l00149"></a>00149 }
<a name="l00150"></a>00150
<a name="l00151"></a>00151 <span class="comment">/*</span>
<a name="l00152"></a>00152 <span class="comment"> * Copyright (c) 1989, 1993</span>
<a name="l00153"></a>00153 <span class="comment"> * The Regents of the University of California. All rights reserved.</span>
<a name="l00154"></a>00154 <span class="comment"> *</span>
<a name="l00155"></a>00155 <span class="comment"> * This code is derived from software posted to USENET.</span>
<a name="l00156"></a>00156 <span class="comment"> *</span>
<a name="l00157"></a>00157 <span class="comment"> * Redistribution and use in source and binary forms, with or without</span>
<a name="l00158"></a>00158 <span class="comment"> * modification, are permitted provided that the following conditions</span>
<a name="l00159"></a>00159 <span class="comment"> * are met:</span>
<a name="l00160"></a>00160 <span class="comment"> * 1. Redistributions of source code must retain the above copyright</span>
<a name="l00161"></a>00161 <span class="comment"> * notice, this list of conditions and the following disclaimer.</span>
<a name="l00162"></a>00162 <span class="comment"> * 2. Redistributions in binary form must reproduce the above copyright</span>
<a name="l00163"></a>00163 <span class="comment"> * notice, this list of conditions and the following disclaimer in the</span>
<a name="l00164"></a>00164 <span class="comment"> * documentation and/or other materials provided with the distribution.</span>
<a name="l00165"></a>00165 <span class="comment"> * 3. All advertising materials mentioning features or use of this software</span>
<a name="l00166"></a>00166 <span class="comment"> * must display the following acknowledgement:</span>
<a name="l00167"></a>00167 <span class="comment"> * This product includes software developed by the University of</span>
<a name="l00168"></a>00168 <span class="comment"> * California, Berkeley and its contributors.</span>
<a name="l00169"></a>00169 <span class="comment"> * 4. Neither the name of the University nor the names of its contributors</span>
<a name="l00170"></a>00170 <span class="comment"> * may be used to endorse or promote products derived from this software</span>
<a name="l00171"></a>00171 <span class="comment"> * without specific prior written permission.</span>
<a name="l00172"></a>00172 <span class="comment"> *</span>
<a name="l00173"></a>00173 <span class="comment"> * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND</span>
<a name="l00174"></a>00174 <span class="comment"> * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE</span>
<a name="l00175"></a>00175 <span class="comment"> * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE</span>
<a name="l00176"></a>00176 <span class="comment"> * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE</span>
<a name="l00177"></a>00177 <span class="comment"> * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL</span>
<a name="l00178"></a>00178 <span class="comment"> * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS</span>
<a name="l00179"></a>00179 <span class="comment"> * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)</span>
<a name="l00180"></a>00180 <span class="comment"> * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT</span>
<a name="l00181"></a>00181 <span class="comment"> * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY</span>
<a name="l00182"></a>00182 <span class="comment"> * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF</span>
<a name="l00183"></a>00183 <span class="comment"> * SUCH DAMAGE.</span>
<a name="l00184"></a>00184 <span class="comment"> */</span>
<a name="l00185"></a>00185
<a name="l00186"></a>00186 <span class="preprocessor">#if HAVE_SYS_CDEFS_H</span>
<a name="l00187"></a>00187 <span class="preprocessor"></span><span class="preprocessor">#include <sys/cdefs.h></span>
<a name="l00188"></a>00188 <span class="preprocessor">#endif</span>
<a name="l00189"></a>00189 <span class="preprocessor"></span>
<a name="l00190"></a>00190 <span class="comment">/*</span>
<a name="l00191"></a>00191 <span class="comment"> * Phase of the Moon. Calculates the current phase of the moon.</span>
<a name="l00192"></a>00192 <span class="comment"> * Based on routines from `Practical Astronomy with Your Calculator',</span>
<a name="l00193"></a>00193 <span class="comment"> * by Duffett-Smith. Comments give the section from the book that</span>
<a name="l00194"></a>00194 <span class="comment"> * particular piece of code was adapted from.</span>
<a name="l00195"></a>00195 <span class="comment"> *</span>
<a name="l00196"></a>00196 <span class="comment"> * -- Keith E. Brandt VIII 1984</span>
<a name="l00197"></a>00197 <span class="comment"> *</span>
<a name="l00198"></a>00198 <span class="comment"> * Updated to the Third Edition of Duffett-Smith's book, Paul Janzen, IX 1998</span>
<a name="l00199"></a>00199 <span class="comment"> *</span>
<a name="l00200"></a>00200 <span class="comment"> */</span>
<a name="l00201"></a>00201
<a name="l00202"></a>00202
<a name="l00203"></a>00203 <span class="preprocessor">#include <ctype.h></span>
<a name="l00204"></a>00204 <span class="preprocessor">#if HAVE_ERR_H</span>
<a name="l00205"></a>00205 <span class="preprocessor"></span><span class="preprocessor">#include <err.h></span>
<a name="l00206"></a>00206 <span class="preprocessor">#endif</span>
<a name="l00207"></a>00207 <span class="preprocessor"></span><span class="preprocessor">#include <math.h></span>
<a name="l00208"></a>00208 <span class="preprocessor">#include <string.h></span>
<a name="l00209"></a>00209 <span class="preprocessor">#include <stdlib.h></span>
<a name="l00210"></a>00210 <span class="preprocessor">#include <time.h></span>
<a name="l00211"></a>00211 <span class="preprocessor">#include <unistd.h></span>
<a name="l00212"></a>00212
<a name="l00213"></a>00213 <span class="preprocessor">#ifndef PI</span>
<a name="l00214"></a>00214 <span class="preprocessor"></span><span class="preprocessor">#define PI 3.14159265358979323846</span>
<a name="l00215"></a>00215 <span class="preprocessor"></span><span class="preprocessor">#endif</span>
<a name="l00216"></a>00216 <span class="preprocessor"></span>
<a name="l00217"></a>00217 <span class="comment">/*</span>
<a name="l00218"></a>00218 <span class="comment"> * The EPOCH in the third edition of the book is 1990 Jan 0.0 TDT.</span>
<a name="l00219"></a>00219 <span class="comment"> * In this program, we do not bother to correct for the differences</span>
<a name="l00220"></a>00220 <span class="comment"> * between UTC (as shown by the UNIX clock) and TDT. (TDT = TAI + 32.184s;</span>
<a name="l00221"></a>00221 <span class="comment"> * TAI-UTC = 32s in Jan 1999.)</span>
<a name="l00222"></a>00222 <span class="comment"> */</span>
<a name="l00223"></a>00223 <span class="preprocessor">#define EPOCH_MINUS_1970 (20 * 365 + 5 - 1) </span><span class="comment">/* 20 years, 5 leaps, back 1 day to Jan 0 */</span>
<a name="l00224"></a>00224 <span class="preprocessor">#define EPSILONg 279.403303 </span><span class="comment">/* solar ecliptic long at EPOCH */</span>
<a name="l00225"></a>00225 <span class="preprocessor">#define RHOg 282.768422 </span><span class="comment">/* solar ecliptic long of perigee at EPOCH */</span>
<a name="l00226"></a>00226 <span class="preprocessor">#define ECCEN 0.016713 </span><span class="comment">/* solar orbit eccentricity */</span>
<a name="l00227"></a>00227 <span class="preprocessor">#define lzero 318.351648 </span><span class="comment">/* lunar mean long at EPOCH */</span>
<a name="l00228"></a>00228 <span class="preprocessor">#define Pzero 36.340410 </span><span class="comment">/* lunar mean long of perigee at EPOCH */</span>
<a name="l00229"></a>00229 <span class="preprocessor">#define Nzero 318.510107 </span><span class="comment">/* lunar mean long of node at EPOCH */</span>
<a name="l00230"></a>00230
<a name="l00231"></a>00231 <span class="comment">/*</span>
<a name="l00232"></a>00232 <span class="comment"> * percentFull --</span>
<a name="l00233"></a>00233 <span class="comment"> * return phase of the moon as a percentage of full</span>
<a name="l00234"></a>00234 <span class="comment"> */</span>
<a name="l00235"></a>00235 <span class="keywordtype">double</span> LunarPhase::percentFull( uint tmpt )<span class="keyword"> const</span>
<a name="l00236"></a>00236 <span class="keyword"></span>{
<a name="l00237"></a>00237 <span class="keywordtype">double</span> N, Msol, Ec, LambdaSol, l, Mm, Ev, Ac, A3, Mmprime;
<a name="l00238"></a>00238 <span class="keywordtype">double</span> A4, lprime, V, ldprime, D, Nm;
<a name="l00239"></a>00239
<a name="l00240"></a>00240 <span class="keywordtype">double</span> days;
<a name="l00241"></a>00241 days = ( tmpt - EPOCH_MINUS_1970 * 86400 ) / 86400.0;
<a name="l00242"></a>00242
<a name="l00243"></a>00243 N = 360 * days / 365.242191; <span class="comment">/* sec 46 #3 */</span>
<a name="l00244"></a>00244 adj360(&N);
<a name="l00245"></a>00245 Msol = N + EPSILONg - RHOg; <span class="comment">/* sec 46 #4 */</span>
<a name="l00246"></a>00246 adj360(&Msol);
<a name="l00247"></a>00247 Ec = 360 / PI * ECCEN * sin(degreesToRadians(Msol)); <span class="comment">/* sec 46 #5 */</span>
<a name="l00248"></a>00248 LambdaSol = N + Ec + EPSILONg; <span class="comment">/* sec 46 #6 */</span>
<a name="l00249"></a>00249 adj360(&LambdaSol);
<a name="l00250"></a>00250 l = 13.1763966 * days + lzero; <span class="comment">/* sec 65 #4 */</span>
<a name="l00251"></a>00251 adj360(&l);
<a name="l00252"></a>00252 Mm = l - (0.1114041 * days) - Pzero; <span class="comment">/* sec 65 #5 */</span>
<a name="l00253"></a>00253 adj360(&Mm);
<a name="l00254"></a>00254 Nm = Nzero - (0.0529539 * days); <span class="comment">/* sec 65 #6 */</span>
<a name="l00255"></a>00255 adj360(&Nm);
<a name="l00256"></a>00256 Ev = 1.2739 * sin(degreesToRadians(2*(l - LambdaSol) - Mm)); <span class="comment">/* sec 65 #7 */</span>
<a name="l00257"></a>00257 Ac = 0.1858 * sin(degreesToRadians(Msol)); <span class="comment">/* sec 65 #8 */</span>
<a name="l00258"></a>00258 A3 = 0.37 * sin(degreesToRadians(Msol));
<a name="l00259"></a>00259 Mmprime = Mm + Ev - Ac - A3; <span class="comment">/* sec 65 #9 */</span>
<a name="l00260"></a>00260 Ec = 6.2886 * sin(degreesToRadians(Mmprime)); <span class="comment">/* sec 65 #10 */</span>
<a name="l00261"></a>00261 A4 = 0.214 * sin(degreesToRadians(2 * Mmprime)); <span class="comment">/* sec 65 #11 */</span>
<a name="l00262"></a>00262 lprime = l + Ev + Ec - Ac + A4; <span class="comment">/* sec 65 #12 */</span>
<a name="l00263"></a>00263 V = 0.6583 * sin(degreesToRadians(2 * (lprime - LambdaSol))); <span class="comment">/* sec 65 #13 */</span>
<a name="l00264"></a>00264 ldprime = lprime + V; <span class="comment">/* sec 65 #14 */</span>
<a name="l00265"></a>00265 D = ldprime - LambdaSol; <span class="comment">/* sec 67 #2 */</span>
<a name="l00266"></a>00266 <span class="keywordflow">return</span>(50.0 * (1 - cos(degreesToRadians(D)))); <span class="comment">/* sec 67 #3 */</span>
<a name="l00267"></a>00267 }
<a name="l00268"></a>00268
<a name="l00269"></a>00269 <span class="comment">/*</span>
<a name="l00270"></a>00270 <span class="comment"> * degreesToRadians --</span>
<a name="l00271"></a>00271 <span class="comment"> * convert degrees to radians</span>
<a name="l00272"></a>00272 <span class="comment"> */</span>
<a name="l00273"></a>00273 <span class="keywordtype">double</span> LunarPhase::degreesToRadians( <span class="keywordtype">double</span> degree )<span class="keyword"> const</span>
<a name="l00274"></a>00274 <span class="keyword"></span>{
<a name="l00275"></a>00275 <span class="keywordflow">return</span>( degree * PI / 180 );
<a name="l00276"></a>00276 }
<a name="l00277"></a>00277
<a name="l00278"></a>00278 <span class="comment">/*</span>
<a name="l00279"></a>00279 <span class="comment"> * adj360 --</span>
<a name="l00280"></a>00280 <span class="comment"> * adjust value so 0 <= degree <= 360</span>
<a name="l00281"></a>00281 <span class="comment"> */</span>
<a name="l00282"></a>00282 <span class="keywordtype">void</span> LunarPhase::adj360( <span class="keywordtype">double</span> *degree )<span class="keyword"> const</span>
<a name="l00283"></a>00283 <span class="keyword"></span>{
<a name="l00284"></a>00284 <span class="keywordflow">for</span>( ;; )
<a name="l00285"></a>00285 <span class="keywordflow">if</span>( *degree < 0 )
<a name="l00286"></a>00286 *degree += 360;
<a name="l00287"></a>00287 <span class="keywordflow">else</span> <span class="keywordflow">if</span>( *degree > 360 )
<a name="l00288"></a>00288 *degree -= 360;
<a name="l00289"></a>00289 <span class="keywordflow">else</span>
<a name="l00290"></a>00290 <span class="keywordflow">break</span>;
<a name="l00291"></a>00291 }
</pre></div></div>
</div></div>
</td>
</tr>
<tr>
<td valign="top" id="leftmenu" width="25%">
<a name="navigation"></a>
<div class="menu_box"><h2>libkholidays</h2>
<div class="nav_list">
<ul><li><a href="index.html">Main Page</a></li><li><a href="classes.html">Alphabetical List</a></li><li><a href="annotated.html">Class List</a></li><li><a href="files.html">File List</a></li><li><a href="functions.html">Class Members</a></li></ul>
<!--
<h2>Class Picker</h2>
<div style="text-align: center;">
<form name="guideform">
<select name="guidelinks" style="width:100%;" onChange="window.location=document.guideform.guidelinks.options[document.guideform.guidelinks.selectedIndex].value">
<option value="annotated.html">-- Choose --</option>
<option value="classLunarPhase.html">lunarphase</option>,
</select>
</form>
</div>
-->
</div></div>
<div class="menu_box"><h2>API Dox</h2>
<div class="nav_list">
<ul>
<li><a href="../../akregator/html/index.html">akregator</a></li><li> <a href="../../akregator/src/html/index.html">src</a></li><li> <a href="../../akregator/src/librss/html/index.html">librss</a></li><li><a href="../../certmanager/html/index.html">certmanager</a></li><li> <a href="../../certmanager/lib/html/index.html">lib</a></li><li><a href="../../kaddressbook/html/index.html">kaddressbook</a></li><li><a href="../../kalarm/html/index.html">kalarm</a></li><li> <a href="../../kalarm/lib/html/index.html">lib</a></li><li><a href="../../kandy/html/index.html">kandy</a></li><li><a href="../../karm/html/index.html">karm</a></li><li><a href="../../kdgantt/html/index.html">kdgantt</a></li><li><a href="../../kgantt/html/index.html">kgantt</a></li><li><a href="../../kioslaves/html/index.html">kioslaves</a></li><li> <a href="../../kioslaves/imap4/html/index.html">imap4</a></li><li> <a href="../../kioslaves/mbox/html/index.html">mbox</a></li><li><a href="../../kitchensync/html/index.html">kitchensync</a></li><li><a href="../../kmail/html/index.html">kmail</a></li><li><a href="../../knotes/html/index.html">knotes</a></li><li><a href="../../konsolekalendar/html/index.html">konsolekalendar</a></li><li><a href="../../kontact/html/index.html">kontact</a></li><li><a href="../../korganizer/html/index.html">korganizer</a></li><li><a href="../../kpilot/html/index.html">kpilot</a></li><li> <a href="../../kpilot/kpilot/html/index.html">kpilot</a></li><li> <a href="../../kpilot/lib/html/index.html">lib</a></li><li><a href="../../libemailfunctions/html/index.html">libemailfunctions</a></li><li><a href="../../libkcal/html/index.html">libkcal</a></li><li><a href="../../libkdepim/html/index.html">libkdepim</a></li><li><a href="../../libkholidays/html/index.html">libkholidays</a></li><li><a href="../../libkmime/html/index.html">libkmime</a></li><li><a href="../../libkpgp/html/index.html">libkpgp</a></li><li><a href="../../libkpimidentities/html/index.html">libkpimidentities</a></li>
</ul></div></div>
</td>
</tr>
</table>
<span class="doNotDisplay">
<a href="http://www.kde.org/" accesskey="8">KDE Home</a> |
<a href="http://accessibility.kde.org/" accesskey="9">KDE Accessibility Home</a> |
<a href="http://www.kde.org/media/accesskeys.php" accesskey="0">Description of Access Keys</a>
</span>
<div style="height: 8px"></div>
<div id="footer">
<div id="footer_left">
Maintained by <a href="mailto:groot@kde.org">Adriaan de Groot</a>
and
<a href="mailto:winter@kde.org">Allen Winter</a>.
<br/>
KDE and K Desktop Environment are trademarks of <a href="http://www.kde.org/areas/kde-ev/" title="Homepage of the KDE non-profit Organization">KDE e.V.</a> |
<a href="http://www.kde.org/contact/impressum.php">Legal</a>
</div>
<div id="footer_right"><img src="/media/images/footer_right.png" style="margin: 0px" alt="" /></div>
</div>
<!--
WARNING: DO NOT SEND MAIL TO THE FOLLOWING EMAIL ADDRESS! YOU WILL
BE BLOCKED INSTANTLY AND PERMANENTLY!
<a href="mailto:aaaatrap-425acc3b5374943f@kde.org">Block me</a>
WARNING END
-->
</body>
</html>
