panchangam 0.1.0

High-precision Vedic Panchangam calculations using Swiss Ephemeris
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115

//! Tithi (Lunar Day) calculations
//! Tithi = (Moon longitude - Sun longitude) / 12°

use wasm_bindgen::prelude::*;
use serde::{Serialize, Deserialize};
use alloc::string::{String, ToString};
use crate::astronomy::planets::{sun_longitude, moon_longitude};

/// Tithi names (1-30)
pub const TITHI_NAMES: [&str; 30] = [
    "Pratipada", "Dwitiya", "Tritiya", "Chaturthi", "Panchami",
    "Shashthi", "Saptami", "Ashtami", "Navami", "Dashami",
    "Ekadashi", "Dwadashi", "Trayodashi", "Chaturdashi", "Purnima",
    "Pratipada", "Dwitiya", "Tritiya", "Chaturthi", "Panchami",
    "Shashthi", "Saptami", "Ashtami", "Navami", "Dashami",
    "Ekadashi", "Dwadashi", "Trayodashi", "Chaturdashi", "Amavasya",
];

/// Paksha (lunar fortnight)
#[derive(Debug, Clone, Copy, PartialEq, Serialize, Deserialize)]
#[wasm_bindgen]
pub enum Paksha {
    Shukla = 0, // Waxing (bright fortnight)
    Krishna = 1, // Waning (dark fortnight)
}

/// Tithi information
#[derive(Debug, Clone, Serialize, Deserialize)]
#[wasm_bindgen(getter_with_clone)]
pub struct TithiInfo {
    /// Tithi index (1-30)
    pub index: u8,
    /// Tithi name
    pub name: String,
    /// Lunar fortnight
    pub paksha: Paksha,
    /// Completion percentage (0.0 to 1.0)
    pub completion: f64,
}

#[wasm_bindgen]
impl TithiInfo {
    #[wasm_bindgen(getter)]
    pub fn paksha_name(&self) -> String {
        match self.paksha {
            Paksha::Shukla => "Shukla".to_string(),
            Paksha::Krishna => "Krishna".to_string(),
        }
    }
}

/// Calculate Tithi for a given Julian Day
/// Returns TithiInfo with index, name, paksha, and completion percentage
#[wasm_bindgen]
pub fn calculate_tithi(jd: f64) -> TithiInfo {
    let sun_long = sun_longitude(jd);
    let moon_long = moon_longitude(jd);
    
    // Angular distance: Moon - Sun
    let mut diff = moon_long - sun_long;
    if diff < 0.0 {
        diff += 360.0;
    }
    
    // Each Tithi spans 12 degrees
    let tithi_float = diff / 12.0;
    let tithi_index = (tithi_float.floor() as u8) + 1;
    
    // Completion within current Tithi
    let completion = tithi_float - tithi_float.floor();
    
    // Paksha: Shukla (1-15), Krishna (16-30)
    let paksha = if tithi_index <= 15 {
        Paksha::Shukla
    } else {
        Paksha::Krishna
    };
    
    // Adjust index for display (1-15 for both fortnights)
    let _display_index = if tithi_index <= 15 { tithi_index } else { tithi_index - 15 };
    
    TithiInfo {
        index: tithi_index,
        name: TITHI_NAMES[(tithi_index - 1) as usize].to_string(),
        paksha,
        completion,
    }
}

/// Calculate the Julian Day when the current Tithi ends
pub fn tithi_end_time(jd: f64) -> f64 {
    let current = calculate_tithi(jd);
    let target_angle = current.index as f64 * 12.0; // End of current Tithi
    
    // Iterative solver (simple Newton-Raphson)
    let mut search_jd = jd;
    for _ in 0..20 {
        let sun_long = sun_longitude(search_jd);
        let moon_long = moon_longitude(search_jd);
        let mut diff = moon_long - sun_long;
        if diff < 0.0 { diff += 360.0; }
        
        let error = target_angle - diff;
        if error.abs() < 0.001 {
            break;
        }
        
        // Moon moves ~13.2 deg/day, Sun moves ~1 deg/day
        // Net motion ~12.2 deg/day
        let step = error / 12.2;
        search_jd += step;
    }
    
    search_jd
}