eddy/mountcontrol
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

...

Browse Source
This commit is contained in:
Timur A. Fatkhullin
2025-08-15 12:18:19 +03:00
parent e934bd3932
commit 61e41b1a1d
2 changed files with 19 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
cxx/mcc_mount_pz.h
View File

@@ -214,6 +214,7 @@ public:
_coord2coord(target.coordPairKind, target.x, target.y, target.time_point,
MccCoordPairKind::COORDS_KIND_HADEC_APP, ha, dec);
// _altLimit = 0.001_degs;
// compute HA for intersection point
double cos_ha =
(std::sin(_altLimit) - std::sin(dec) * std::sin(_latitude)) / std::cos(dec) / std::cos(_latitude);
@@ -223,25 +224,36 @@ public:
std::cout << "HA(85.0) = " << MccAngle(acos(cos_ha)).sexagesimal(true) << "\n";
std::cout << "HA(85.0) = " << MccAngle(-acos(cos_ha)).sexagesimal(true) << "\n";
std::cout << "PHI = " << MccAngle(_latitude).sexagesimal() << "\n";
std::cout << "COS_HA = " << cos_ha << "\n";
if (cos_ha > 1.0) { // no intersection
// compute culmination points?
return false;
}
double cosA =
(-std::sin(dec) * std::cos(_latitude) + std::cos(dec) * std::sin(_latitude) * cos_ha) / cos(_altLimit);
// cosA /= std::cos(_altLimit);
std::cout << "COS_A = " << cosA << "\n";
double sinA = std::cos(dec) * sqrt(1.0 - cos_ha * cos_ha) / cos(_altLimit);
double tgA = sqrt(1.0 - cos_ha * cos_ha) / (cos(_latitude) * tan(dec) - sin(_latitude) * cos_ha);
auto z = (-std::sin(dec) * std::cos(_latitude) + std::cos(dec) * std::sin(_latitude) * cos_ha) / cosA;
// (-std::sin(dec) * std::cos(_latitude) + std::cos(dec) * std::sin(_latitude) * cos_ha) / cos(133.75_degs);
std::cout << "Z = " << MccAngle(asin(z)).sexagesimal() << "\n";
if constexpr (KIND ==
MccAltLimitKind::MIN_ALT_LIMIT) { // the closest time point is one after upper culmination
az = -std::acos(cosA);
az = std::acos(cosA);
// az = atan2(sinA, cosA);
} else if constexpr (KIND == MccAltLimitKind::MAX_ALT_LIMIT) { // the closest time point is one before upper
// culmination
az = std::acos(cosA);
az = -std::acos(cosA) + std::numbers::pi; // to system of azimuth started from the North!!!
// az = atan(tgA);
// az = std::asin(sinA);
// az = atan2(sinA, cosA);
}
@@ -249,7 +261,6 @@ public:
_coord2coord(MccCoordPairKind::COORDS_KIND_AZALT, az, _altLimit, target.time_point, int_point.coordPairKind,
int_point.x, int_point.y);
int_point.x = az;
return true;
}