panchangam 0.1.2

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

//! Nakshatra (Lunar Mansion) calculations
//! Moon's longitude divided into 27 segments of 13°20' each

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

/// Nakshatra names (27 lunar mansions)
pub const NAKSHATRA_NAMES: [&str; 27] = [
    "Ashwini", "Bharani", "Krittika", "Rohini", "Mrigashirsha",
    "Ardra", "Punarvasu", "Pushya", "Ashlesha", "Magha",
    "Purva Phalguni", "Uttara Phalguni", "Hasta", "Chitra", "Swati",
    "Vishakha", "Anuradha", "Jyeshtha", "Mula", "Purva Ashadha",
    "Uttara Ashadha", "Shravana", "Dhanishta", "Shatabhisha", "Purva Bhadrapada",
    "Uttara Bhadrapada", "Revati",
];

/// Nakshatra lords (planetary rulers)
pub const NAKSHATRA_LORDS: [&str; 27] = [
    "Ketu", "Venus", "Sun", "Moon", "Mars",
    "Rahu", "Jupiter", "Saturn", "Mercury", "Ketu",
    "Venus", "Sun", "Moon", "Mars", "Rahu",
    "Jupiter", "Saturn", "Mercury", "Ketu", "Venus",
    "Sun", "Moon", "Mars", "Rahu", "Jupiter",
    "Saturn", "Mercury",
];

/// Nakshatra qualities
pub const NAKSHATRA_QUALITIES: [&str; 27] = [
    "Swift/Light", "Fierce/Severe", "Mixed", "Fixed/Permanent", "Soft/Gentle",
    "Sharp/Dreadful", "Movable/Ephemeral", "Light/Swift", "Sharp/Dreadful", "Fierce/Severe",
    "Fierce/Severe", "Fixed/Permanent", "Light/Swift", "Soft/Gentle", "Movable/Ephemeral",
    "Mixed", "Soft/Gentle", "Sharp/Dreadful", "Sharp/Dreadful", "Fierce/Severe",
    "Fixed/Permanent", "Movable/Ephemeral", "Movable/Ephemeral", "Movable/Ephemeral", "Fierce/Severe",
    "Fixed/Permanent", "Soft/Gentle",
];

/// Nakshatra information
#[derive(Debug, Clone, Serialize, Deserialize)]
#[wasm_bindgen(getter_with_clone)]
pub struct NakshatraInfo {
    /// Nakshatra index (1-27)
    pub index: u8,
    /// Nakshatra name
    pub name: String,
    /// Planetary ruler
    pub ruler: String,
    /// Quality/nature
    pub quality: String,
    /// Pada (quarter, 1-4)
    pub pada: u8,
}

/// Calculate Nakshatra for a given Julian Day
/// Each Nakshatra spans 13°20' (13.333... degrees)
#[wasm_bindgen]
pub fn calculate_nakshatra(jd: f64, ayanamsha_mode: AyanamshaMode) -> NakshatraInfo {
    let moon_tropical = moon_longitude(jd);
    let ayanamsha = get_ayanamsha(ayanamsha_mode, jd);
    let moon_sidereal = tropical_to_sidereal(moon_tropical, ayanamsha);
    
    // Each Nakshatra = 360/27 = 13.333... degrees
    let nakshatra_span = 360.0 / 27.0;
    let nakshatra_float = moon_sidereal / nakshatra_span;
    let nakshatra_index = (nakshatra_float.floor() as u8) + 1;
    
    // Pada (quarter within Nakshatra)
    let within_nakshatra = nakshatra_float - nakshatra_float.floor();
    let pada = (within_nakshatra * 4.0).floor() as u8 + 1;
    
    let idx = (nakshatra_index - 1) as usize;
    
    NakshatraInfo {
        index: nakshatra_index,
        name: NAKSHATRA_NAMES[idx].to_string(),
        ruler: NAKSHATRA_LORDS[idx].to_string(),
        quality: NAKSHATRA_QUALITIES[idx].to_string(),
        pada,
    }
}

/// Calculate Julian Day when Moon transitions to next Nakshatra
pub fn nakshatra_end_time(jd: f64, ayanamsha_mode: AyanamshaMode) -> f64 {
    let current = calculate_nakshatra(jd, ayanamsha_mode);
    let nakshatra_span = 360.0 / 27.0;
    let target_long = current.index as f64 * nakshatra_span; // End of current Nakshatra
    
    // Iterative solver
    let mut search_jd = jd;
    for _ in 0..20 {
        let moon_tropical = moon_longitude(search_jd);
        let ayanamsha = get_ayanamsha(ayanamsha_mode, search_jd);
        let moon_sidereal = tropical_to_sidereal(moon_tropical, ayanamsha);
        
        let error = target_long - moon_sidereal;
        let error = if error < -180.0 { error + 360.0 } else if error > 180.0 { error - 360.0 } else { error };
        
        if error.abs() < 0.001 {
            break;
        }
        
        // Moon moves ~13.2 deg/day
        let step = error / 13.2;
        search_jd += step;
    }
    
    search_jd
}