mihrab-methods 0.0.0-stage → 1.0.0

package/dist/esm/sun/preciseSun.js

@@ -0,0 +1,181 @@
+/**
+ * A precise apparent Sun: the Earth's VSOP87 series (truncated), IAU 2006
+ * precession, IAU 2000B nutation and light time, after Don Cross's
+ * Astronomy Engine (MIT licence). It agrees with JPL DE421 to about 0.05 s
+ * in the equation of time and 0.5″ in declination, where Meeus's short
+ * series (`meeusApparentSun`) drifts by about a second over the year —
+ * enough to move a published minute.
+ */
+const ASEC2RAD = 4.848136811095359935899141e-6;
+const ASEC360 = 1296000;
+const DEG = Math.PI / 180;
+const C_AUDAY = 173.1446326846693;
+const DAYS_PER_MILLENNIUM = 365250;
+/** Heliocentric ecliptic (J2000) longitude, latitude and radius of the Earth: [ampl, phase, freq] per power of t. */
+const EARTH = [
+    [
+        [
+            [1.75347045673, 0, 0], [0.03341656453, 4.66925680415, 6283.0758499914], [0.00034894275, 4.62610242189, 12566.1516999828],
+            [0.00003417572, 2.82886579754, 3.523118349], [0.00003497056, 2.74411783405, 5753.3848848968], [0.00003135899, 3.62767041756, 77713.7714681205],
+            [0.00002676218, 4.41808345438, 7860.4193924392], [0.00002342691, 6.13516214446, 3930.2096962196], [0.00001273165, 2.03709657878, 529.6909650946],
+            [0.00001324294, 0.74246341673, 11506.7697697936], [0.00000901854, 2.04505446477, 26.2983197998], [0.00001199167, 1.10962946234, 1577.3435424478],
+            [0.00000857223, 3.50849152283, 398.1490034082], [0.00000779786, 1.17882681962, 5223.6939198022], [0.0000099025, 5.23268072088, 5884.9268465832],
+            [0.00000753141, 2.53339052847, 5507.5532386674], [0.00000505267, 4.58292599973, 18849.2275499742], [0.00000492392, 4.20505711826, 775.522611324],
+            [0.00000356672, 2.91954114478, 0.0673103028], [0.00000284125, 1.89869240932, 796.2980068164], [0.00000242879, 0.34481445893, 5486.777843175],
+            [0.00000317087, 5.84901948512, 11790.6290886588], [0.00000271112, 0.31486255375, 10977.078804699], [0.00000206217, 4.80646631478, 2544.3144198834],
+            [0.00000205478, 1.86953770281, 5573.1428014331], [0.00000202318, 2.45767790232, 6069.7767545534], [0.00000126225, 1.08295459501, 20.7753954924],
+            [0.00000155516, 0.83306084617, 213.299095438],
+        ],
+        [[6283.0758499914, 0, 0], [0.00206058863, 2.67823455808, 6283.0758499914], [0.00004303419, 2.63512233481, 12566.1516999828]],
+        [[0.00008721859, 1.07253635559, 6283.0758499914]],
+    ],
+    [[], [[0.00227777722, 3.4137662053, 6283.0758499914], [0.00003805678, 3.37063423795, 12566.1516999828]]],
+    [
+        [
+            [1.00013988784, 0, 0], [0.01670699632, 3.09846350258, 6283.0758499914], [0.00013956024, 3.05524609456, 12566.1516999828],
+            [0.0000308372, 5.19846674381, 77713.7714681205], [0.00001628463, 1.17387558054, 5753.3848848968], [0.00001575572, 2.84685214877, 7860.4193924392],
+            [0.00000924799, 5.45292236722, 11506.7697697936], [0.00000542439, 4.56409151453, 3930.2096962196], [0.0000047211, 3.66100022149, 5884.9268465832],
+            [0.00000085831, 1.27079125277, 161000.6857376741], [0.00000057056, 2.01374292245, 83996.84731811189], [0.00000055736, 5.2415979917, 71430.69561812909],
+            [0.00000174844, 3.01193636733, 18849.2275499742], [0.00000243181, 4.2734953079, 11790.6290886588],
+        ],
+        [[0.00103018607, 1.10748968172, 6283.0758499914], [0.00001721238, 1.06442300386, 12566.1516999828]],
+        [[0.00004359385, 5.78455133808, 6283.0758499914]],
+    ],
+];
+function vsop(formula, t, clamp) {
+    let tp = 1;
+    let out = 0;
+    for (const series of formula) {
+        let sum = 0;
+        for (const [a, p, f] of series)
+            sum += a * Math.cos(p + t * f);
+        let incr = tp * sum;
+        if (clamp)
+            incr %= 2 * Math.PI;
+        out += incr;
+        tp *= t;
+    }
+    return out;
+}
+/** Earth−J2000 equatorial position (AU) at `tt` days from J2000. */
+function earthJ2000(tt) {
+    const t = tt / DAYS_PER_MILLENNIUM;
+    const lon = vsop(EARTH[0], t, true);
+    const lat = vsop(EARTH[1], t, false);
+    const r = vsop(EARTH[2], t, false);
+    const x = r * Math.cos(lat) * Math.cos(lon);
+    const y = r * Math.cos(lat) * Math.sin(lon);
+    const z = r * Math.sin(lat);
+    return [
+        x + 0.00000044036 * y - 0.000000190919 * z,
+        -0.000000479966 * x + 0.917482137087 * y - 0.397776982902 * z,
+        0.397776982902 * y + 0.917482137087 * z,
+    ];
+}
+/** ΔT (s), Espenak & Meeus, for UT days from J2000 — held at its 2017 value onward (about 69 s, as observed), as JPL Horizons does. */
+export function deltaT(ut) {
+    const y = 2000 + (Math.min(ut, 17 * 365.24217) - 14) / 365.24217;
+    if (y >= 2005 && y < 2050) {
+        const u = y - 2000;
+        return 62.92 + 0.32217 * u + 0.005589 * u * u;
+    }
+    if (y >= 2050 && y < 2150) {
+        const u = (y - 1820) / 100;
+        return -20 + 32 * u * u - 0.5628 * (2150 - y);
+    }
+    if (y >= 1986 && y < 2005) {
+        const u = y - 2000;
+        return 63.86 + 0.3345 * u - 0.060374 * u ** 2 + 0.0017275 * u ** 3 + 0.000651814 * u ** 4 + 0.00002373599 * u ** 5;
+    }
+    const u = (y - 1820) / 100;
+    return -20 + 32 * u * u;
+}
+function nutation(tt) {
+    const mod = (x) => (x % ASEC360) * ASEC2RAD;
+    const t = tt / 36525;
+    const elp = mod(1287104.79305 + t * 129596581.0481);
+    const f = mod(335779.526232 + t * 1739527262.8478);
+    const d = mod(1072260.70369 + t * 1602961601.209);
+    const om = mod(450160.398036 - t * 6962890.5431);
+    let dp = (-172064161.0 - 174666.0 * t) * Math.sin(om) + 33386.0 * Math.cos(om);
+    let de = (92052331.0 + 9086.0 * t) * Math.cos(om) + 15377.0 * Math.sin(om);
+    let a = 2 * (f - d + om);
+    dp += (-13170906.0 - 1675.0 * t) * Math.sin(a) - 13696.0 * Math.cos(a);
+    de += (5730336.0 - 3015.0 * t) * Math.cos(a) - 4587.0 * Math.sin(a);
+    a = 2 * (f + om);
+    dp += (-2276413.0 - 234.0 * t) * Math.sin(a) + 2796.0 * Math.cos(a);
+    de += (978459.0 - 485.0 * t) * Math.cos(a) + 1374.0 * Math.sin(a);
+    a = 2 * om;
+    dp += (2074554.0 + 207.0 * t) * Math.sin(a) - 698.0 * Math.cos(a);
+    de += (-897492.0 + 470.0 * t) * Math.cos(a) - 291.0 * Math.sin(a);
+    dp += (1475877.0 - 3633.0 * t) * Math.sin(elp) + 11817.0 * Math.cos(elp);
+    de += (73871.0 - 184.0 * t) * Math.cos(elp) - 1924.0 * Math.sin(elp);
+    return { dpsi: -0.000135 + dp * 1e-7, deps: 0.000388 + de * 1e-7 };
+}
+function meanObliquity(tt) {
+    const t = tt / 36525;
+    return ((((-0.0000000434 * t - 0.000000576) * t + 0.0020034) * t - 0.0001831) * t - 46.836769) * t + 84381.406;
+}
+/** J2000 equatorial → mean equator and equinox of date (IAU 2006). */
+function precess(v, tt) {
+    const t = tt / 36525;
+    const eps0 = 84381.406 * ASEC2RAD;
+    const psia = (((((-0.0000000951 * t + 0.000132851) * t - 0.00114045) * t - 1.0790069) * t + 5038.481507) * t) * ASEC2RAD;
+    const omegaa = (((((0.0000003337 * t - 0.000000467) * t - 0.00772503) * t + 0.0512623) * t - 0.025754) * t + 84381.406) * ASEC2RAD;
+    const chia = (((((-0.000000056 * t + 0.000170663) * t - 0.00121197) * t - 2.3814292) * t + 10.556403) * t) * ASEC2RAD;
+    const sa = Math.sin(eps0), ca = Math.cos(eps0);
+    const sb = Math.sin(-psia), cb = Math.cos(-psia);
+    const sc = Math.sin(-omegaa), cc = Math.cos(-omegaa);
+    const sd = Math.sin(chia), cd = Math.cos(chia);
+    const xx = cd * cb - sb * sd * cc;
+    const yx = cd * sb * ca + sd * cc * cb * ca - sa * sd * sc;
+    const zx = cd * sb * sa + sd * cc * cb * sa + ca * sd * sc;
+    const xy = -sd * cb - sb * cd * cc;
+    const yy = -sd * sb * ca + cd * cc * cb * ca - sa * cd * sc;
+    const zy = -sd * sb * sa + cd * cc * cb * sa + ca * cd * sc;
+    const xz = sb * sc;
+    const yz = -sc * cb * ca - sa * cc;
+    const zz = -sc * cb * sa + cc * ca;
+    return [xx * v[0] + yx * v[1] + zx * v[2], xy * v[0] + yy * v[1] + zy * v[2], xz * v[0] + yz * v[1] + zz * v[2]];
+}
+/** Mean → true equator and equinox of date. */
+function nutate(v, oblm, oblt, psi) {
+    const cobm = Math.cos(oblm), sobm = Math.sin(oblm);
+    const cobt = Math.cos(oblt), sobt = Math.sin(oblt);
+    const cpsi = Math.cos(psi), spsi = Math.sin(psi);
+    const xx = cpsi;
+    const yx = -spsi * cobm;
+    const zx = -spsi * sobm;
+    const xy = spsi * cobt;
+    const yy = cpsi * cobm * cobt + sobm * sobt;
+    const zy = cpsi * sobm * cobt - cobm * sobt;
+    const xz = spsi * sobt;
+    const yz = cpsi * cobm * sobt - sobm * cobt;
+    const zz = cpsi * sobm * sobt + cobm * cobt;
+    return [xx * v[0] + yx * v[1] + zx * v[2], xy * v[0] + yy * v[1] + zy * v[2], xz * v[0] + yz * v[1] + zz * v[2]];
+}
+/**
+ * The apparent Sun at `d` days from J2000.0 (UT): declination (rad) and
+ * equation of time (minutes), such that the transit at longitude λ (east) is
+ * 720 − 4λ − eot minutes UT. Same contract as `meeusApparentSun`.
+ */
+export function preciseApparentSun(d) {
+    const tt = d + deltaT(d) / 86400;
+    const e = earthJ2000(tt - 1 / C_AUDAY);
+    const nut = nutation(tt);
+    const mobl = meanObliquity(tt);
+    const oblm = (mobl / 3600) * DEG;
+    const oblt = ((mobl + nut.deps) / 3600) * DEG;
+    const v = nutate(precess([-e[0], -e[1], -e[2]], tt), oblm, oblt, nut.dpsi * ASEC2RAD);
+    const ra = Math.atan2(v[1], v[0]) / DEG;
+    const dec = Math.atan2(v[2], Math.hypot(v[0], v[1]));
+    // Greenwich apparent sidereal time (degrees): Earth rotation angle + the precession and equation of the equinoxes.
+    const t = tt / 36525;
+    const era = 360 * (((0.779057273264 + 0.00273781191135448 * d + (d % 1)) % 1) + 1);
+    const eqeq = nut.dpsi * Math.cos(oblm);
+    const st = eqeq + 0.014506 + ((((-0.0000000368 * t - 0.000029956) * t - 0.00000044) * t + 1.3915817) * t + 4612.156534) * t;
+    const gast = era + st / 3600;
+    const tmin = (d - Math.floor(d + 0.5) + 0.5) * 1440; // minutes past 0h UT
+    const eotDeg = ((((gast - ra - (tmin - 720) / 4 + 180) % 360) + 360) % 360) - 180;
+    return { dec, eot: 4 * eotDeg };
+}

package/dist/esm/sweden.d.ts

@@ -0,0 +1,4 @@
+/** mihrab-methods/sweden — the sweden method with its places and correction list, and nothing else's. */
+import './authorities/sweden.js';
+export * from './core/index.js';
+export { places } from './authorities/sweden.js';

package/dist/esm/sweden.js

@@ -0,0 +1,4 @@
+/** mihrab-methods/sweden — the sweden method with its places and correction list, and nothing else's. */
+import './authorities/sweden.js';
+export * from './core/index.js';
+export { places } from './authorities/sweden.js';

package/docs/methods.md

@@ -0,0 +1,190 @@
+# The methods, in words and formulas
+
+What each authority's program does, as read off its published tables, so
+that a port in another language can do the same. The TypeScript under
+`src/methods/` is the reference: where this page and the code differ, the
+code is what the certification holds. `test/vectors.json` has a small sample
+of places, dates and the times to expect, with the clock they print in.
+
+Conventions throughout:
+
+- Angles in degrees unless marked; φ the place's latitude, λ its longitude
+  (east positive), δ the sun's declination, E the equation of time in
+  minutes (apparent solar time less mean), T Julian centuries since J2000.0.
+- The **hour angle** at which the sun's centre is at altitude *a*:
+  `cos H = (sin a − sin φ · sin δ) / (cos φ · cos δ)`, H in degrees, 4 minutes per degree.
+- **Transit** (local apparent noon) in minutes past 0h UT: `720 − 4λ − E`.
+- A time that does not form (the sun never reaches the altitude that day)
+  is null unless the method says otherwise.
+- The **printed minute** is the instant on the table's clock. A table prints
+  the whole day in the clock kept at noon of that day, so on the morning the
+  clocks go forward a Fajr before the change is printed in the new time.
+- Each method is followed by its **correction list**
+  (`data/corrections/<authority>.json`): the few times the authority's own
+  arithmetic puts a fraction of a second on the other side of a minute's
+  rounding point. An entry moves one printed time of one point on one date by
+  one minute. The lists are not part of the method; they are what remains
+  after it.
+
+## Sweden and Diyanet
+
+`src/methods/sweden.ts`. The International Prayer Times Congress method
+(Istanbul, 2021 — the Diyanet and the European Council for Fatwa and
+Research) as Diyanet's program realises it. Islamiska Förbundet's Swedish
+tables are Diyanet's, row for row. Diyanet's worldwide tables are the same
+program with Isha's angle by country.
+
+**Published rules** (awqatsalah.com, "Tespit Kriterleri"):
+
+- Dhuhr five minutes after transit. Asr at one shadow length, four minutes
+  late; where the shadow never forms, Asr is printed at Dhuhr's time.
+- Sunrise seven minutes early, Maghrib seven minutes late ("temkin"), at
+  the standard −0.833°.
+- True Fajr at −18°, true Isha at −16° (the Congress) or −17° (Diyanet's own,
+  outside the countries listed below).
+- Neither the day nor the night is shorter than five hours.
+- Where the 18° dawn no longer forms, the share of the night that the sharia
+  night's third took on the last day it did is carried: Isha that share of
+  the night after Maghrib, Fajr the same share scaled by 18/16 before sunrise.
+- The move from the true times to the carried ones is gradual, starting
+  twenty minutes before (Isha) or after (Fajr) the takdir.
+
+**What the tables do with them**, read off Diyanet's own tables for 160
+Swedish places and 8 000 places worldwide, 2024–2027:
+
+- The **sun is taken once a day**, at 0h UT less 45 seconds, from Meeus's
+  mean elements:
+  - `L0 = 280.46646 + 36000.76983 T + 0.0003032 T²` (mean longitude)
+  - `M = 357.52911 + 35999.05029 T − 0.0001537 T²` (mean anomaly)
+  - `C = (1.914602 − 0.004817 T − 0.000014 T²) sin M + (0.019993 − 0.000101 T) sin 2M + 0.000289 sin 3M`
+  - true longitude `λ☉ = L0 + C` (the *true*, not the apparent, longitude:
+    no aberration, no nutation in longitude);
+  - obliquity `ε = 23°26′21.448″ − 46.815″ T − 0.00059″ T² + 0.001813″ T³`
+    plus the nutation in obliquity `9.2 cos Ω + 0.57 cos 2L☉ + 0.1 cos 2L☾ − 0.09 cos 2Ω` arc-seconds,
+    with `Ω = 125.04 − 1934.136 T`, `L☉ = 280.4665 + 36000.7698 T`, `L☾ = 218.3165 + 481267.8813 T`;
+  - right ascension `α = atan2(cos ε · sin λ☉, cos λ☉)`, declination `δ = asin(sin ε · sin λ☉)`;
+  - equation of time `E = 4 (L0 − α)` minutes, wrapped to ±720, plus 0.00176 min;
+  - and two small waves in the mean longitude the formula leaves out, read
+    off the tables (`eotResidualS`, `decResidualDeg` in the code): a few
+    hundredths of a second on the transit, a tenth of an arc-second on the
+    declination. A port that omits them will be a minute off on a few times
+    in ten thousand.
+- Coordinates are kept to **three decimals** of a degree.
+- The night a share is taken of is the day's own: 24 h less the printed day
+  (sunrise to Maghrib, after temkin).
+- The **five-hour floors** hold each end separately: 150 minutes from the
+  printed Dhuhr, or from the midnight after it. Where midsummer's night is
+  floored at both ends (to within a minute), that night anchors the season,
+  and the floored times never step back: before midsummer no sunrise earlier
+  than midsummer's and a floored Maghrib no earlier than the day before;
+  after it no sunrise earlier nor Maghrib later than on any day since
+  midsummer — and Fajr and Isha carried the same way.
+- The **gradual transition** is a straight line across the days next to the
+  gap on which the true time lies within twenty minutes of the takdir at the
+  gap's edge: twenty minutes away at the far end, the takdir at the edge.
+- Where the dawn forms all summer but midsummer's true Isha (at 16°) falls
+  later than the share carried from that night, **midsummer is the takdir's
+  day** — a gap of one day, with its ramps (from about 44.4° north).
+- **Isha by country:** 16° in most of Europe, Russia, Ukraine, Belarus,
+  Kazakhstan, Mongolia, China, the United States and Canada; 17° elsewhere.
+  The high-latitude rules are worked at 16° either way.
+- **Türkiye's districts** (and Northern Cyprus's): each prayer is computed at
+  a point of its own in the district, and the points move by a few
+  thousandths of a degree from one year's tables to the next
+  (`Place.prayerPoints`: a point per prayer per year; the nearest year's
+  beyond the tables).
+- **Past local midnight:** where the first true dawn after the summer gap
+  falls before local midnight on the table's clock (far east of a zone's
+  meridian), there is no transition: the gap's last day is printed twenty
+  minutes after the takdir and the true times follow.
+- **Tables on their own clocks:** some of Diyanet's tables are printed on a
+  clock other than their place's zone (Alberta keeping summer time, Xinjiang
+  on Beijing time, and so on); `Place.clock` gives that clock, read off the
+  table. The instants are the method's; only the labels differ.
+- Every time is rounded to the nearest minute after the temkin.
+
+The full algorithm — the search for the gap, the anchoring, the ramps — is
+`CongressYear` in the code, about four hundred lines; there is no shorter
+statement of it that gives the tables.
+
+## Morocco
+
+`src/methods/morocco.ts`. The Ministry of Habous and Islamic Affairs, read
+off its tables for 164 cities.
+
+- The sun is **Meeus's apparent sun** with the full nutation (the mean
+  elements above; apparent longitude `L0 + C − 0.00569 + Δψ`, with
+  `Δψ = −17.2 sin Ω − 1.32 sin 2L☉ − 0.23 sin 2L☾ + 0.21 sin 2Ω` arc-seconds)
+  and his chapter-28 equation of time
+  `E = L0 − 0.0057183 − α + Δψ cos ε`, taken at **each prayer's own
+  moment** (iterated), declination and equation of time alike.
+- Fajr −19°, Isha −17°, Asr one shadow length (Shafiʿi).
+- Dhuhr five minutes after the transit.
+- Sunrise 2.85 minutes before, and Maghrib 3.7 minutes after, the sun's
+  upper limb on the horizon: altitude `−(34′ + s) − dip`, with s the sun's
+  semi-diameter for its distance that day and dip the city's own horizon
+  dip — eastward at sunrise (`Place.dipDeg`), westward at Maghrib
+  (`Place.dipMaghribDeg`), each fitted from the tables with the city's
+  point; where not fitted, `0.02°·√height(m)`, and 0.86 of it at Maghrib.
+- Fajr rounded **down** to the minute, every other time to the nearest.
+- The ministry's clock is GMT+1 until 17 September 2026 and GMT from then
+  on (`Place.clock`); the zone databases still move Morocco's clock, so the
+  table's clock is not a zone's.
+
+## Algeria
+
+`src/methods/algeria.ts`. The Ministry of Religious Affairs and Wakfs, read
+off its 1448 AH table for 68 cities.
+
+- The sun from **NOAA's solar-calculator formulas**: the mean elements above
+  with the eccentricity `e = 0.016708634 − 0.000042037 T − 0.0000001267 T²`,
+  NOAA's series for the equation of time
+  `E = 4 [ y sin 2L0 − 2e sin M + 4e y sin M cos 2L0 − ½ y² sin 4L0 − 1¼ e² sin 2M ]`
+  with `y = tan²(ε₀/2)` (the mean obliquity), and its declination from the
+  apparent longitude `L0 + C − 0.00569 − 0.00478 sin Ω` and the obliquity
+  `ε₀ + 0.00256 cos Ω`. They fit the table far better than a precise
+  ephemeris, so this is what the ministry's program uses.
+- Dhuhr at the sun's transit, the equation of time taken there.
+- Every other time is the transit plus or minus the hour angle, with the
+  declination taken at that prayer's own moment (three iterations) and the
+  equation of time kept at noon's value (itself iterated three times).
+- Fajr −18°, Isha −17°, sunrise and sunset at −0.833°.
+- Asr from the noon shadow `tan(φ − δ)` **with its sign**, so where the sun
+  culminates north of the zenith (the far south around June) Asr comes
+  earlier, up to thirteen minutes.
+- Maghrib three minutes after sunset.
+- Every time rounded **up** to the minute.
+
+## Indonesia
+
+`src/methods/indonesia.ts`. Kemenag's tables for all 517 regencies and
+cities.
+
+- The **precise apparent sun** (`src/sun/preciseSun.ts`: a VSOP87 series
+  after Astronomy Engine, within 0.05 s and 0.5″ of JPL's DE421), taken at
+  each prayer's own moment, iterated three times. Meeus's shorter series
+  drifts by a second and leaves a minute in two hundred off.
+- Fajr −20°, Isha −18°.
+- Sunrise and Maghrib at −1°, and at −2° in the regencies the ministry
+  classes as highland (`Place.highland`); in six highland regencies of Papua
+  the degree is taken off Fajr and Isha too (`Place.twilightDip`).
+- Asr when the shadow is its length plus the noon shadow, the noon shadow
+  from the sun's declination at 00:00 UT of the day (as an almanac's daily
+  entry gives it), the hour angle from the sun at Asr (four iterations).
+- Every time rounded **up**, after a margin: +2 minutes on Fajr, Asr,
+  Maghrib and Isha, +3 on Dhuhr, −3 on sunrise.
+- Each regency at its own fitted point; Karawang's Asr, Maghrib and Isha at a
+  second point about 2 km north-east (`Place.afternoon`).
+- Each regency is printed in its province's zone (WIB, WITA or WIT).
+
+## Checking a port
+
+1. Take `test/vectors.json`. For each case, compute the six times at the
+   given point (with the place's own fields: dips, highland class, Isha
+   angle, per-prayer points, clock) and compare the instants and the printed
+   minutes.
+2. The vectors include the days the clocks change and a table on its own
+   clock; a port that prints the wrong day's offset shows up there.
+3. Then apply the correction list and compare against the authority's own
+   published table for a place and a month. The list is small on purpose: if
+   a port needs more than it, the port is not the method.

package/package.json

@@ -1,6 +1,180 @@
 {
   "name": "mihrab-methods",
-  "version": "0.0.0-stage",
-  "stub": true,
-  "description": "Temporary package placeholder for staged publishing"
-}
+  "version": "1.0.0",
+  "description": "The prayer-time authorities' own calculation methods, each reproducing its published timetable to the minute: Sweden (Diyanet, as Islamiska Förbundet prints it), Diyanet worldwide, Morocco (Habous), Algeria (Ministry of Religious Affairs), Indonesia (Kemenag).",
+  "license": "MPL-2.0",
+  "author": "Hassan El Ghamri",
+  "repository": {
+    "type": "git",
+    "url": "git+https://github.com/MihrabHQ/mihrab-methods.git"
+  },
+  "homepage": "https://github.com/MihrabHQ/mihrab-methods#readme",
+  "bugs": "https://github.com/MihrabHQ/mihrab-methods/issues",
+  "keywords": [
+    "prayer times",
+    "salah",
+    "namaz",
+    "diyanet",
+    "habous",
+    "kemenag",
+    "islamiska förbundet",
+    "adhan",
+    "mihrab"
+  ],
+  "sideEffects": [
+    "./dist/*/index.js",
+    "./dist/*/sweden.js",
+    "./dist/*/diyanet.js",
+    "./dist/*/morocco.js",
+    "./dist/*/algeria.js",
+    "./dist/*/indonesia.js",
+    "./dist/*/authorities/*.js"
+  ],
+  "engines": {
+    "node": ">=18"
+  },
+  "files": [
+    "dist",
+    "schemas",
+    "docs/methods.md",
+    "certification.json",
+    "data/corrections"
+  ],
+  "main": "./dist/cjs/index.js",
+  "module": "./dist/esm/index.js",
+  "types": "./dist/cjs/index.d.ts",
+  "exports": {
+    ".": {
+      "import": {
+        "types": "./dist/esm/index.d.ts",
+        "default": "./dist/esm/index.js"
+      },
+      "require": {
+        "types": "./dist/cjs/index.d.ts",
+        "default": "./dist/cjs/index.js"
+      }
+    },
+    "./core": {
+      "import": {
+        "types": "./dist/esm/core.d.ts",
+        "default": "./dist/esm/core.js"
+      },
+      "require": {
+        "types": "./dist/cjs/core.d.ts",
+        "default": "./dist/cjs/core.js"
+      }
+    },
+    "./sweden": {
+      "import": {
+        "types": "./dist/esm/sweden.d.ts",
+        "default": "./dist/esm/sweden.js"
+      },
+      "require": {
+        "types": "./dist/cjs/sweden.d.ts",
+        "default": "./dist/cjs/sweden.js"
+      }
+    },
+    "./diyanet": {
+      "import": {
+        "types": "./dist/esm/diyanet.d.ts",
+        "default": "./dist/esm/diyanet.js"
+      },
+      "require": {
+        "types": "./dist/cjs/diyanet.d.ts",
+        "default": "./dist/cjs/diyanet.js"
+      }
+    },
+    "./morocco": {
+      "import": {
+        "types": "./dist/esm/morocco.d.ts",
+        "default": "./dist/esm/morocco.js"
+      },
+      "require": {
+        "types": "./dist/cjs/morocco.d.ts",
+        "default": "./dist/cjs/morocco.js"
+      }
+    },
+    "./algeria": {
+      "import": {
+        "types": "./dist/esm/algeria.d.ts",
+        "default": "./dist/esm/algeria.js"
+      },
+      "require": {
+        "types": "./dist/cjs/algeria.d.ts",
+        "default": "./dist/cjs/algeria.js"
+      }
+    },
+    "./indonesia": {
+      "import": {
+        "types": "./dist/esm/indonesia.d.ts",
+        "default": "./dist/esm/indonesia.js"
+      },
+      "require": {
+        "types": "./dist/cjs/indonesia.d.ts",
+        "default": "./dist/cjs/indonesia.js"
+      }
+    },
+    "./locate": {
+      "import": {
+        "types": "./dist/esm/locate.d.ts",
+        "default": "./dist/esm/locate.js"
+      },
+      "require": {
+        "types": "./dist/cjs/locate.d.ts",
+        "default": "./dist/cjs/locate.js"
+      }
+    },
+    "./certification.json": "./certification.json",
+    "./data/corrections/*.json": "./data/corrections/*.json",
+    "./schemas/*.json": "./schemas/*.json",
+    "./package.json": "./package.json"
+  },
+  "typesVersions": {
+    "*": {
+      "core": [
+        "./dist/cjs/core.d.ts"
+      ],
+      "sweden": [
+        "./dist/cjs/sweden.d.ts"
+      ],
+      "diyanet": [
+        "./dist/cjs/diyanet.d.ts"
+      ],
+      "morocco": [
+        "./dist/cjs/morocco.d.ts"
+      ],
+      "algeria": [
+        "./dist/cjs/algeria.d.ts"
+      ],
+      "indonesia": [
+        "./dist/cjs/indonesia.d.ts"
+      ],
+      "locate": [
+        "./dist/cjs/locate.d.ts"
+      ]
+    }
+  },
+  "scripts": {
+    "embed": "node scripts/embed.mjs",
+    "build": "npm run embed && rm -rf dist && tsc -p tsconfig.build.cjs.json && tsc -p tsconfig.build.esm.json && echo '{\"type\":\"module\"}' > dist/esm/package.json",
+    "prepack": "npm run build",
+    "pretest": "npm run embed",
+    "test": "tsx --test test/*.test.ts",
+    "pretypecheck": "npm run embed",
+    "typecheck": "tsc --noEmit",
+    "preverify": "npm run embed",
+    "verify": "tsx tools/verify.ts",
+    "precorrections": "npm run embed",
+    "corrections": "tsx tools/corrections.ts",
+    "precertify": "npm run embed",
+    "certify": "tsx tools/certify.ts",
+    "vectors": "npm run embed && tsx tools/vectors.ts",
+    "serve": "npm run embed && tsx tools/serve.ts"
+  },
+  "devDependencies": {
+    "@types/node": "^22.10.0",
+    "ajv": "^8.20.0",
+    "tsx": "^4.19.2",
+    "typescript": "^5.7.2"
+  }
+}

package/schemas/certification.schema.json

@@ -0,0 +1,59 @@
+{
+  "$schema": "https://json-schema.org/draft/2020-12/schema",
+  "$id": "https://raw.githubusercontent.com/MihrabHQ/mihrab-methods/main/schemas/certification.schema.json",
+  "title": "mihrab-methods certification",
+  "description": "What may be used: for each authority, whether its method reproduces every time of the authority's tables on the printed minute, over which dates, at which method revision, and which correction list that holds with. Written by `npm run certify`; read by apps to keep current (fetchLatestCorrections).",
+  "type": "object",
+  "propertyNames": { "enum": ["sweden", "diyanet", "morocco", "algeria", "indonesia"] },
+  "additionalProperties": { "$ref": "#/$defs/entry" },
+  "$defs": {
+    "date": { "type": "string", "pattern": "^[0-9]{4}-[0-9]{2}-[0-9]{2}$" },
+    "entry": {
+      "type": "object",
+      "required": ["certified", "checkedAt", "checks"],
+      "additionalProperties": false,
+      "properties": {
+        "certified": { "type": "boolean", "description": "True only when every check below has every time on the printed minute. Not trusted alone: apps re-derive it from the checks." },
+        "checkedAt": { "$ref": "#/$defs/date" },
+        "checks": { "type": "array", "minItems": 1, "items": { "$ref": "#/$defs/check" } },
+        "corrections": {
+          "type": "object",
+          "description": "The correction list this certification is for: its identity, and the SHA-256 of data/corrections/<authority>.json as published.",
+          "required": ["from", "to", "entries"],
+          "additionalProperties": false,
+          "properties": {
+            "from": { "$ref": "#/$defs/date" },
+            "to": { "$ref": "#/$defs/date" },
+            "entries": { "type": "integer", "minimum": 0 },
+            "sha256": { "type": "string", "pattern": "^[0-9a-f]{64}$" }
+          }
+        },
+        "revision": { "type": "integer", "minimum": 1, "description": "The method's revision the list is for. An app computing another revision must not apply the list." }
+      }
+    },
+    "check": {
+      "type": "object",
+      "required": ["tables", "places", "times", "off", "from", "to"],
+      "additionalProperties": false,
+      "properties": {
+        "tables": { "type": "string", "description": "Whose tables were checked." },
+        "places": { "type": "integer", "minimum": 0 },
+        "times": { "type": "integer", "minimum": 0 },
+        "off": { "type": "integer", "minimum": 0, "description": "Times not on the printed minute. Certified only at 0." },
+        "from": { "$ref": "#/$defs/date" },
+        "to": { "$ref": "#/$defs/date" },
+        "setAside": {
+          "type": "object",
+          "description": "Rows set aside as the source's fault, not the method's, so each can be judged.",
+          "required": ["misplaced", "sourceFaults", "copiedRows"],
+          "additionalProperties": false,
+          "properties": {
+            "misplaced": { "type": "array", "items": { "type": "string" } },
+            "sourceFaults": { "type": "integer", "minimum": 0 },
+            "copiedRows": { "type": "integer", "minimum": 0 }
+          }
+        }
+      }
+    }
+  }
+}

package/schemas/corrections.schema.json

@@ -0,0 +1,29 @@
+{
+  "$schema": "https://json-schema.org/draft/2020-12/schema",
+  "$id": "https://raw.githubusercontent.com/MihrabHQ/mihrab-methods/main/schemas/corrections.schema.json",
+  "title": "mihrab-methods correction list",
+  "description": "A method's correction list: the few times its authority's own program puts on the other side of a minute's rounding point. Each entry moves one printed time of one point on one date by a minute. Published at data/corrections/<authority>.json; certification.json names the list in force (its dates, its size and its SHA-256).",
+  "type": "object",
+  "required": ["from", "to", "points"],
+  "additionalProperties": false,
+  "properties": {
+    "from": { "$ref": "#/$defs/date", "description": "The first date the list was made over." },
+    "to": { "$ref": "#/$defs/date", "description": "The last date the list was made over; past it the method alone applies." },
+    "points": {
+      "type": "object",
+      "description": "The point the method is computed at, as 'lat,lon' to three decimals → its corrections.",
+      "propertyNames": { "pattern": "^-?[0-9]+\\.[0-9]{3},-?[0-9]+\\.[0-9]{3}$" },
+      "additionalProperties": {
+        "type": "array",
+        "items": {
+          "type": "string",
+          "description": "'YYYY-MM-DD prayer ±1': the minutes to add to that prayer's printed time on that date.",
+          "pattern": "^[0-9]{4}-[0-9]{2}-[0-9]{2} (fajr|sunrise|dhuhr|asr|maghrib|isha) [+-]1$"
+        }
+      }
+    }
+  },
+  "$defs": {
+    "date": { "type": "string", "pattern": "^[0-9]{4}-[0-9]{2}-[0-9]{2}$" }
+  }
+}