alignment/SPK: complete pointing model with altaz convention and polar axis mapping

Implement the full 6-term Pmfit least-squares fitting loop with hot-path
incremental normal-equation accumulation. Apply the altaz roll convention
(pi - Az, matching vtel Note 7). Add call to MathPlugin::Initialise() so
the base class sanitizes sync points before Pmfit runs. Activate EQUATORIAL_PE
output in the CCD/guide simulators only when tracking state is Ok, so plate
solves see a stable star field.

Map ME (polar elevation error, pmv[2]) to pan/AN (along-meridian tilt) and
MA (polar azimuth error, pmv[3]) to paw/AW (across-meridian tilt) in
ParsePmfitCoefficients, following Wallace PM struct conventions.

Use RA-based (ANTI_CLOCKWISE) TDV encoding throughout: telescope_simulator
Sync/Goto, MathPlugin polar sanitizer, and SPKMathPlugin boundary functions.
Store LST in UpdateAstrometry so DirectionVectorToRollPitch and
RollPitchToDirectionVector can convert between RA-based TDVs and the
HA-based roll/pitch that the Wallace kernel expects internally.

Add altaz/EQ polar degeneracy and regression tests. Update ARCHITECTURE.md
to document the Wallace error injection pipeline.

Author: ckemper67

drivers/ccd/ccd_simulator.cpp

@@ -658,6 +658,12 @@ int CCDSim::DrawCcdFrame(INDI::CCDChip * targetChip)
         }
         else
         {
+            static bool s_warnedNoSnoop = false;
+            if (!s_warnedNoSnoop)
+            {
+                LOG_WARN("EQUATORIAL_PE not yet snooped — rendering at raw EQUATORIAL_EOD_COORD (mount errors not active)");
+                s_warnedNoSnoop = true;
+            }
 #endif
             currentRA  = RA;
             currentDE = Dec;
@@ -1382,7 +1388,15 @@ bool CCDSim::ISSnoopDevice(XMLEle * root)
             INDI::ObservedToJ2000(&epochPos, ln_get_julian_from_sys(), &J2000Pos);
             raPE  = J2000Pos.rightascension;
             decPE = J2000Pos.declination;
-            usePE = true;
+
+            // Only activate PE rendering once the telescope reports a valid (Ok) state.
+            // The initial defNumberVector arrives with state Idle and values 0/0.
+            const char * state = findXMLAttValu(root, "state");
+            if (!usePE && strcmp(state, "Ok") == 0)
+            {
+                LOGF_INFO("EQUATORIAL_PE snoop active: JNow RA %.4f Dec %.4f -> J2000 RA %.4f Dec %.4f", newra, newdec, raPE, decPE);
+                usePE = true;
+            }
 
             EqPENP[AXIS_RA].setValue(newra);
             EqPENP[AXIS_DE].setValue(newdec);

drivers/ccd/guide_simulator.cpp

@@ -519,6 +519,12 @@ int GuideSim::DrawCcdFrame(INDI::CCDChip * targetChip)
         }
         else
         {
+            static bool s_warnedNoSnoop = false;
+            if (!s_warnedNoSnoop)
+            {
+                LOG_WARN("EQUATORIAL_PE not yet snooped — rendering at raw EQUATORIAL_EOD_COORD (mount errors not active)");
+                s_warnedNoSnoop = true;
+            }
 #endif
             m_CurrentRA  = RA;
             m_CurrentDEC = Dec;
@@ -1168,7 +1174,15 @@ bool GuideSim::ISSnoopDevice(XMLEle * root)
             INDI::ObservedToJ2000(&epochPos, ln_get_julian_from_sys(), &J2000Pos);
             raPE  = J2000Pos.rightascension;
             decPE = J2000Pos.declination;
-            m_UsePE = true;
+
+            // Only activate PE rendering once the telescope reports a valid (Ok) state.
+            // The initial defNumberVector arrives with state Idle and values 0/0.
+            const char * state = findXMLAttValu(root, "state");
+            if (!m_UsePE && strcmp(state, "Ok") == 0)
+            {
+                LOGF_INFO("EQUATORIAL_PE snoop active: JNow RA %.4f Dec %.4f -> J2000 RA %.4f Dec %.4f", newra, newdec, raPE, decPE);
+                m_UsePE = true;
+            }
 
             EqPENP[AXIS_RA].setValue(newra);
             EqPENP[AXIS_DE].setValue(newdec);

drivers/telescope/telescope_simulator.cpp

@@ -579,6 +579,12 @@ bool ScopeSim::Goto(double r, double d)
             EquatorialCoordinatesFromTelescopeDirectionVector(tdv, raCoords);
             targetRA  = raCoords.rightascension;
             targetDec = raCoords.declination;
+            LOGF_INFO("EQ Goto: alignment corrected RA %.4fh Dec %.4f -> RA %.4fh Dec %.4f",
+                      r, d, targetRA, targetDec);
+        }
+        else
+        {
+            LOG_WARN("EQ Goto: TransformCelestialToTelescope FAILED, using raw coordinates");
         }
     }