utils.js

const THIEN_CAN_LIST = ['Giáp', 'Ất', 'Bính', 'Đinh', 'Mậu', 'Kỷ', 'Canh', 'Tân', 'Nhâm', 'Quý'];
const DIA_CHI_LIST = ['Tí', 'Sửu', 'Dần', 'Mão', 'Thìn', 'Tị', 'Ngọ', 'Mùi', 'Thân', 'Dậu', 'Tuất', 'Hợi'];
const TIET_KHI_LIST = [
    "Lập Xuân", "Vũ Thủy", "Kinh Trập", "Xuân Phân",
    "Thanh Minh", "Cốc Vũ", "Lập Hạ", "Tiểu Mãn",
    "Mang Chủng", "Hạ Chí", "Tiểu Thử", "Đại Thử",
    "Lập Thu", "Xử Thử", "Bạch Lộ", "Thu Phân",
    "Hàn Lộ", "Sương Giáng", "Lập Đông", "Tiểu Tuyết",
    "Đại Tuyết", "Đông Chí", "Tiểu Hàn", "Đại Hàn"
];

function jdFromDate(dd, mm, yy) {
    let a, y, m, jd;
    a = Math.floor((14 - mm) / 12);
    y = yy + 4800 - a;
    m = mm + 12*a - 3;
    jd = dd + Math.floor((153*m+2)/5) + 365*y + Math.floor(y/4) - Math.floor(y/100) + Math.floor(y/400) - 32045;
    if (jd < 2299161) {
        jd = dd + Math.floor((153*m+2)/5) + 365*y + Math.floor(y/4) - 32083;
    }
    return jd;
}

function jdToDate(jd) {
    let a, b, c, d, e, m, day, month, year;
    if (jd > 2299160) { // After 5/10/1582, Gregorian calendar
        a = jd + 32044;
        b = Math.floor((4*a+3)/146097);
        c = a - Math.floor((b*146097)/4);
    } else {
        b = 0;
        c = jd + 32082;
    }
    d = Math.floor((4*c+3)/1461);
    e = c - Math.floor((1461*d)/4);
    m = Math.floor((5*e+2)/153);
    day = e - Math.floor((153*m+2)/5) + 1;
    month = m + 3 - 12*Math.floor(m/10);
    year = b*100 + d - 4800 + Math.floor(m/10);
    return {day, month, year};
}

function getNewMoonDay(k, timeZone) {
    let T, T2, T3, dr, Jd1, M, Mpr, F, C1, deltat, JdNew;
    T = k/1236.85; // Time in Julian centuries from 1900 January 0.5
    T2 = T * T;
    T3 = T2 * T;
    dr = Math.PI/180;
    Jd1 = 2415020.75933 + 29.53058868*k + 0.0001178*T2 - 0.000000155*T3;
    Jd1 = Jd1 + 0.00033*Math.sin((166.56 + 132.87*T - 0.009173*T2)*dr); // Mean new moon
    M = 359.2242 + 29.10535608*k - 0.0000333*T2 - 0.00000347*T3; // Sun's mean anomaly
    Mpr = 306.0253 + 385.81691806*k + 0.0107306*T2 + 0.00001236*T3; // Moon's mean anomaly
    F = 21.2964 + 390.67050646*k - 0.0016528*T2 - 0.00000239*T3; // Moon's argument of latitude
    C1=(0.1734 - 0.000393*T)*Math.sin(M*dr) + 0.0021*Math.sin(2*dr*M);
    C1 = C1 - 0.4068*Math.sin(Mpr*dr) + 0.0161*Math.sin(dr*2*Mpr);
    C1 = C1 - 0.0004*Math.sin(dr*3*Mpr);
    C1 = C1 + 0.0104*Math.sin(dr*2*F) - 0.0051*Math.sin(dr*(M+Mpr));
    C1 = C1 - 0.0074*Math.sin(dr*(M-Mpr)) + 0.0004*Math.sin(dr*(2*F+M));
    C1 = C1 - 0.0004*Math.sin(dr*(2*F-M)) - 0.0006*Math.sin(dr*(2*F+Mpr));
    C1 = C1 + 0.0010*Math.sin(dr*(2*F-Mpr)) + 0.0005*Math.sin(dr*(2*Mpr+M));
    if (T < -11) {
        deltat= 0.001 + 0.000839*T + 0.0002261*T2 - 0.00000845*T3 - 0.000000081*T*T3;
    } else {
        deltat= -0.000278 + 0.000265*T + 0.000262*T2;
    };
    JdNew = Jd1 + C1 - deltat;
    return Math.floor(JdNew + 0.5 + timeZone/24);
}

function getSunLongitude(jdn, timeZone) {
    let T, T2, dr, M, L0, DL, L;
    T = (jdn - 2451545.5 - timeZone/24) / 36525; // Time in Julian centuries from 2000-01-01 12:00:00 GMT
    T2 = T*T;
    dr = Math.PI/180; // degree to radian
    M = 357.52910 + 35999.05030*T - 0.0001559*T2 - 0.00000048*T*T2; // mean anomaly, degree
    L0 = 280.46645 + 36000.76983*T + 0.0003032*T2; // mean longitude, degree
    DL = (1.914600 - 0.004817*T - 0.000014*T2)*Math.sin(dr*M);
    DL = DL + (0.019993 - 0.000101*T)*Math.sin(dr*2*M) + 0.000290*Math.sin(dr*3*M);
    L = L0 + DL; // true longitude, degree
    L = L*dr;
    L = L - Math.PI*2*(Math.floor(L/(Math.PI*2))); // Normalize to (0, 2*PI)
    return Math.floor(L / Math.PI * 6);
}

function getLunarMonth11(yy, timeZone) {
    let k, off, nm, sunLong;
    off = jdFromDate(31, 12, yy) - 2415021;
    k = Math.floor(off / 29.530588853);
    nm = getNewMoonDay(k, timeZone);
    sunLong = getSunLongitude(nm, timeZone); // sun longitude at local midnight
    if (sunLong >= 9) {
        nm = getNewMoonDay(k-1, timeZone);
    }
    return nm;
}

function getLeapMonthOffset(a11, timeZone) {
    let k, last, arc, i;
    k = Math.floor((a11 - 2415021.076998695) / 29.530588853 + 0.5);
    last = 0;
    i = 1; // We start with the month following lunar month 11
    arc = getSunLongitude(getNewMoonDay(k+i, timeZone), timeZone);
    do {
        last = arc;
        i++;
        arc = getSunLongitude(getNewMoonDay(k+i, timeZone), timeZone);
    } while (arc != last && i < 14);
    return i-1;
}

function convertSolar2Lunar(dd, mm, yy, timeZone) {
    let k, dayNumber, monthStart, a11, b11, lunarDay, lunarMonth, lunarYear, lunarLeap;
    dayNumber = jdFromDate(dd, mm, yy);
    k = Math.floor((dayNumber - 2415021.076998695) / 29.530588853);
    monthStart = getNewMoonDay(k+1, timeZone);
    if (monthStart > dayNumber) {
        monthStart = getNewMoonDay(k, timeZone);
    }
    a11 = getLunarMonth11(yy, timeZone);
    b11 = a11;
    if (a11 >= monthStart) {
        lunarYear = yy;
        a11 = getLunarMonth11(yy-1, timeZone);
    } else {
        lunarYear = yy+1;
        b11 = getLunarMonth11(yy+1, timeZone);
    }
    lunarDay = dayNumber-monthStart+1;
    diff = Math.floor((monthStart - a11)/29);
    lunarLeap = 0;
    lunarMonth = diff+11;
    if (b11 - a11 > 365) {
        leapMonthDiff = getLeapMonthOffset(a11, timeZone);
        if (diff >= leapMonthDiff) {
            lunarMonth = diff + 10;
            if (diff == leapMonthDiff) {
                lunarLeap = 1;
            }
        }
    }
    if (lunarMonth > 12) {
        lunarMonth = lunarMonth - 12;
    }
    if (lunarMonth >= 11 && diff < 4) {
        lunarYear -= 1;
    }
    const canChiDate = getCanChiDateByJd(dayNumber);
    const canChiMonth = getCanChiMonth(lunarMonth, lunarYear)
    const canChiYear = getCanChiYear(lunarYear);
    return {lunarDay, lunarMonth, lunarYear, lunarLeap, canChiDate, canChiMonth, canChiYear};
}

function convertLunar2Solar({lunarDay, lunarMonth, lunarYear, lunarLeap, timeZone}) {
    let k, a11, b11, off, leapOff, leapMonth, monthStart;
    if (lunarMonth < 11) {
        a11 = getLunarMonth11(lunarYear-1, timeZone);
        b11 = getLunarMonth11(lunarYear, timeZone);
    } else {
        a11 = getLunarMonth11(lunarYear, timeZone);
        b11 = getLunarMonth11(lunarYear+1, timeZone);
    }
    off = lunarMonth - 11;
    if (off < 0) {
        off += 12;
    }
    if (b11 - a11 > 365) {
        leapOff = getLeapMonthOffset(a11, timeZone);
        leapMonth = leapOff - 2;
        if (leapMonth < 0) {
            leapMonth += 12;
        }
        if (lunarLeap != 0 && lunarMonth != leapMonth) {
            return new Array(0, 0, 0);
        } else if (lunarLeap != 0 || off >= leapOff) {
            off += 1;
        }
    }
    k = Math.floor(0.5 + (a11 - 2415021.076998695) / 29.530588853);
    monthStart = getNewMoonDay(k+off, timeZone);
    return jdToDate(monthStart+lunarDay-1);
}

function getCanChiYear(year) {
    return `${THIEN_CAN_LIST[(year + 6) % 10]} ${DIA_CHI_LIST[(year + 8) % 12]}`;
}

function getCanChiDateByJd(jd) {
    const can = THIEN_CAN_LIST[(jd + 9) % 10];
    const chi = DIA_CHI_LIST[(jd + 1) % 12];
    return `${can} ${chi}`;
}

function getCanChiMonth(lunarMonth, lunarYear) {
    const can = THIEN_CAN_LIST[(lunarYear * 12 + lunarMonth + 3) % 10];
    const chi = DIA_CHI_LIST[(lunarMonth + 1) % 12];
    return `${can} ${chi}`;
}

function calculateLichLapXuan(year) {
    // Dữ liệu cơ bản:
    const baseYear = 1900;
    const baseDate = new Date(baseYear, 1, 4); // Lập Xuân của năm 1900 là 4 tháng 2
    const daysPerYear = 365.2422; // Số ngày trung bình trong năm theo quỹ đạo Trái Đất
  
    // Tính toán số năm cách từ năm 1900
    const differenceInYears = year - baseYear;
  
    // Tính tổng số ngày đã qua kể từ năm 1900
    const totalDaysPassed = differenceInYears * daysPerYear;
  
    // Tạo đối tượng ngày mới dựa trên số ngày đã qua
    const lapXuanDate = new Date(baseDate.getTime() + totalDaysPassed * 24 * 60 * 60 * 1000);
  
    // Trả về ngày Lập Xuân (kết quả là một đối tượng Date)
    return lapXuanDate;
}

function tietKhiIndexToDiaChiMonthIndex(index) {
    if (index < 0 || index >= TIET_KHI_LIST.length) {
        throw new RangeError("index out of range");
    }

    return Math.floor((index + 2) / 2) + 1;
}

function calculateTietKhi(inputDate) {
    const year = inputDate.getFullYear();
    const lichLapXuan = calculateLichLapXuan(year);

    // Tính số ngày đã trôi qua kể từ ngày Lập Xuân
    const daysSinceLapXuan = Math.floor((inputDate - lichLapXuan) / (1000 * 60 * 60 * 24));

    // Xác định tiết khí dựa trên số ngày đã trôi qua
    const tietKhiIndex = Math.floor(daysSinceLapXuan / 15.2184);
    let realTietKhiIndex = tietKhiIndex;
    if (realTietKhiIndex < 0) {
        realTietKhiIndex = TIET_KHI_LIST.length + tietKhiIndex;
    } else if (tietKhiIndex >= TIET_KHI_LIST.length - 1) {
        realTietKhiIndex = tietKhiIndex % TIET_KHI_LIST.length;
    }

    const tietKhi = TIET_KHI_LIST[realTietKhiIndex];

    const diaChiMonthIndex = tietKhiIndexToDiaChiMonthIndex(realTietKhiIndex);
    const virtualLunarMonth = diaChiMonthIndex - 1;
    let virtualLunarYear = year;
    if (inputDate < lichLapXuan) {
        virtualLunarYear -= 1;
    }
    const canChiYear = getCanChiYear(virtualLunarYear);
    const canChiMonth = getCanChiMonth(virtualLunarMonth, virtualLunarYear);

    return {
        'tietKhi': tietKhi,
        'canChiYear': canChiYear,
        'canChiMonth': canChiMonth
    };
}

module.exports = {
    convertSolar2Lunar,
    convertLunar2Solar,
    calculateLichLapXuan,
    calculateTietKhi,
}