| 1 | /* | 
|---|
| 2 | INDI Developers Manual | 
|---|
| 3 | Tutorial #2 | 
|---|
| 4 |  | 
|---|
| 5 | "Simple Telescope Simulator" | 
|---|
| 6 |  | 
|---|
| 7 | In this tutorial, we create a simple simulator. We will use a few handy utility functions provided by INDI | 
|---|
| 8 | like timers, string <---> number conversion, and more. | 
|---|
| 9 |  | 
|---|
| 10 | Refer to README, which contains instruction on how to build this driver, and use it | 
|---|
| 11 | with an INDI-compatible client. | 
|---|
| 12 |  | 
|---|
| 13 | */ | 
|---|
| 14 |  | 
|---|
| 15 | /** \file tutorial_two.c | 
|---|
| 16 | \brief Implement a simple telescope simulator using more complex INDI concepts. | 
|---|
| 17 | \author Jasem Mutlaq | 
|---|
| 18 | */ | 
|---|
| 19 |  | 
|---|
| 20 | /* Standard headers */ | 
|---|
| 21 | #include <stdio.h> | 
|---|
| 22 | #include <stdlib.h> | 
|---|
| 23 | #include <string.h> | 
|---|
| 24 | #include <math.h> | 
|---|
| 25 | #include <unistd.h> | 
|---|
| 26 | #include <sys/time.h> | 
|---|
| 27 |  | 
|---|
| 28 | #include <sys/time.h> | 
|---|
| 29 | #include <time.h> | 
|---|
| 30 |  | 
|---|
| 31 | void show_runtime(int state) { | 
|---|
| 32 | static struct timeval tv; | 
|---|
| 33 | struct timeval tv1; | 
|---|
| 34 | double x, y; | 
|---|
| 35 |  | 
|---|
| 36 | if (state) { | 
|---|
| 37 | gettimeofday(&tv, NULL); | 
|---|
| 38 | } else { | 
|---|
| 39 | gettimeofday(&tv1, NULL); | 
|---|
| 40 | fprintf(stderr, "Ran for: %fmsec\n", (double)(tv1.tv_sec * 1000.0 + tv1.tv_usec / 1000.0) - (double)(tv.tv_sec * 1000.0 + tv.tv_usec / 1000.0)); | 
|---|
| 41 | } | 
|---|
| 42 | } | 
|---|
| 43 |  | 
|---|
| 44 | /* INDI Core headers */ | 
|---|
| 45 |  | 
|---|
| 46 | /* indidevapi.h contains API declerations */ | 
|---|
| 47 | #include "indidevapi.h" | 
|---|
| 48 |  | 
|---|
| 49 | /* INDI Eventloop mechanism */ | 
|---|
| 50 | #include "eventloop.h" | 
|---|
| 51 |  | 
|---|
| 52 | /* INDI Common Routines */ | 
|---|
| 53 | #include "indicom.h" | 
|---|
| 54 |  | 
|---|
| 55 | /* Definitions */ | 
|---|
| 56 |  | 
|---|
| 57 | #define mydev           "Telescope Simulator"           /* Device name */ | 
|---|
| 58 | #define MAIN_GROUP      "Main Control"                  /* Group name */ | 
|---|
| 59 | #define SLEWRATE        1                               /* slew rate, degrees/s */ | 
|---|
| 60 | #define POLLMS          250                             /* poll period, ms */ | 
|---|
| 61 | #define SIDRATE         0.004178                        /* sidereal rate, degrees/s */ | 
|---|
| 62 |  | 
|---|
| 63 | /* Function protptypes */ | 
|---|
| 64 | static void connectTelescope (void); | 
|---|
| 65 | static void mountSim (void *); | 
|---|
| 66 |  | 
|---|
| 67 | /* operational info */ | 
|---|
| 68 | static double targetRA; | 
|---|
| 69 | static double targetDEC; | 
|---|
| 70 |  | 
|---|
| 71 | /* main connection switch | 
|---|
| 72 | Note that the switch will appear to the user as On and Off (versus Connect and Disconnect in tutorial one) | 
|---|
| 73 | Nevertheless, the members _names_ are still CONNECT and DISCONNECT and therefore this is a perfectly legal standard property decleration | 
|---|
| 74 | */ | 
|---|
| 75 | static ISwitch connectS[] = { | 
|---|
| 76 | {"CONNECT",  "On",  ISS_OFF, 0, 0}, {"DISCONNECT", "Off", ISS_ON, 0, 0}}; | 
|---|
| 77 |  | 
|---|
| 78 | static ISwitchVectorProperty connectSP = { mydev, "CONNECTION", "Connection",  MAIN_GROUP, IP_RW, ISR_1OFMANY, 0, IPS_IDLE,  connectS, NARRAY(connectS), "", 0 }; | 
|---|
| 79 |  | 
|---|
| 80 | /* Equatorial position. EQUATORIAL_EOD_COORD is one of INDI's reserved Standard  Properties */ | 
|---|
| 81 | static INumber eqN[] = { | 
|---|
| 82 | /* 1st member is Right ascension */ | 
|---|
| 83 | {"RA"                           /* 1st Number name */ | 
|---|
| 84 | ,"RA  H:M:S"                    /* Number label */ | 
|---|
| 85 | , "%10.6m"                      /* Format. Refer to the indiapi.h for details on number formats */ | 
|---|
| 86 | ,0.                                     /* Minimum value */ | 
|---|
| 87 | , 24.                                   /* Maximum value */ | 
|---|
| 88 | , 0.                                    /* Steps */ | 
|---|
| 89 | , 0.                                    /* Initial value */ | 
|---|
| 90 | , 0                                     /* Pointer to parent, we don't use it, so set it to 0 */ | 
|---|
| 91 | , 0                                     /* Auxiluary member, set it to 0 */ | 
|---|
| 92 | , 0},                                   /* Autxiluar member, set it to 0 */ | 
|---|
| 93 |  | 
|---|
| 94 | /* 2nd member is Declination */ | 
|---|
| 95 | {"DEC", "Dec D:M:S", "%10.6m", -90., 90., 0., 0., 0, 0, 0} | 
|---|
| 96 | }; | 
|---|
| 97 |  | 
|---|
| 98 | static INumberVectorProperty eqNP = {  mydev, "EQUATORIAL_EOD_COORD", "Equatorial JNow",  MAIN_GROUP , IP_RO, 0, IPS_IDLE,  eqN, NARRAY(eqN), "", 0}; | 
|---|
| 99 |  | 
|---|
| 100 |  | 
|---|
| 101 | /* Equatorial EOD Coord Request. This property is for requesting changes to target equatorial coordinates. However, the CURRENT coordinates are reported in EQUATORIAL_EOD_COORDS above.*/ | 
|---|
| 102 | static INumber eqNR[] = {{"RA" ,"RA  H:M:S" , "%10.6m" ,0. , 24., 0., 0., 0, 0, 0}, | 
|---|
| 103 | {"DEC", "Dec D:M:S", "%10.6m", -90., 90., 0., 0., 0, 0, 0}}; | 
|---|
| 104 | static INumberVectorProperty eqNPR = {  mydev, "EQUATORIAL_EOD_COORD_REQUEST", "Equatorial Request",  MAIN_GROUP , IP_RW, 0, IPS_IDLE,  eqNR, NARRAY(eqNR), "", 0}; | 
|---|
| 105 |  | 
|---|
| 106 | /* Property naming convention. All property names are lower case with a postfix to indicate their type. connectS is a switch, | 
|---|
| 107 | * connectSP is a switch vector. eqN is a number, eqNP is a number property, and so on. While this is not strictly required, it makes the code easier to read. */ | 
|---|
| 108 |  | 
|---|
| 109 | #define currentRA       eqN[0].value            /* scope's current simulated RA, rads. Handy macro to right ascension from eqN[] */ | 
|---|
| 110 | #define currentDec      eqN[1].value            /* scope's current simulated Dec, rads. Handy macro to declination from eqN[] */ | 
|---|
| 111 |  | 
|---|
| 112 | /******************************************** | 
|---|
| 113 | Property: Movement (Arrow keys on handset). North/South | 
|---|
| 114 | *********************************************/ | 
|---|
| 115 | static ISwitch MovementNSS[]       = {{"MOTION_NORTH", "North", ISS_OFF, 0, 0}, {"MOTION_SOUTH", "South", ISS_OFF, 0, 0}}; | 
|---|
| 116 | ISwitchVectorProperty MovementNSSP      = { mydev, "TELESCOPE_MOTION_NS", "North/South", MAIN_GROUP, IP_RW, ISR_ATMOST1, 0, IPS_IDLE, MovementNSS, NARRAY(MovementNSS), "", 0}; | 
|---|
| 117 |  | 
|---|
| 118 | /******************************************** | 
|---|
| 119 | Property: Movement (Arrow keys on handset). West/East | 
|---|
| 120 | *********************************************/ | 
|---|
| 121 | static ISwitch MovementWES[]       = {{"MOTION_WEST", "West", ISS_OFF, 0, 0}, {"MOTION_EAST", "East", ISS_OFF, 0, 0}}; | 
|---|
| 122 | ISwitchVectorProperty MovementWESP      = { mydev, "TELESCOPE_MOTION_WE", "West/East", MAIN_GROUP, IP_RW, ISR_ATMOST1, 0, IPS_IDLE, MovementWES, NARRAY(MovementWES), "", 0}; | 
|---|
| 123 |  | 
|---|
| 124 | /* Initlization routine */ | 
|---|
| 125 | static void mountInit() | 
|---|
| 126 | { | 
|---|
| 127 | static int inited;              /* set once mountInit is called */ | 
|---|
| 128 |  | 
|---|
| 129 | if (inited) | 
|---|
| 130 | return; | 
|---|
| 131 |  | 
|---|
| 132 | /* start timer to simulate mount motion | 
|---|
| 133 | The timer will call function mountSim after POLLMS milliseconds */ | 
|---|
| 134 | IEAddTimer (POLLMS, mountSim, NULL); | 
|---|
| 135 |  | 
|---|
| 136 | inited = 1; | 
|---|
| 137 |  | 
|---|
| 138 | } | 
|---|
| 139 |  | 
|---|
| 140 | /* send client definitions of all properties */ | 
|---|
| 141 | void ISGetProperties (const char *dev) | 
|---|
| 142 | { | 
|---|
| 143 | if (dev && strcmp (mydev, dev)) | 
|---|
| 144 | return; | 
|---|
| 145 |  | 
|---|
| 146 | IDDefSwitch (&connectSP, NULL); | 
|---|
| 147 | IDDefNumber (&eqNP, NULL); | 
|---|
| 148 | IDDefNumber (&eqNPR, NULL); | 
|---|
| 149 | IDDefSwitch (&MovementNSSP, NULL); | 
|---|
| 150 | IDDefSwitch (&MovementWESP, NULL); | 
|---|
| 151 |  | 
|---|
| 152 | } | 
|---|
| 153 |  | 
|---|
| 154 | void ISNewText (const char *dev, const char *name, char *texts[], char *names[], int n) | 
|---|
| 155 | { | 
|---|
| 156 | return; | 
|---|
| 157 | } | 
|---|
| 158 |  | 
|---|
| 159 | /* client is sending us a new value for a Numeric vector property */ | 
|---|
| 160 | void ISNewNumber (const char *dev, const char *name, double values[], char *names[], int n) | 
|---|
| 161 | { | 
|---|
| 162 | /* Make sure to initalize */ | 
|---|
| 163 | mountInit(); | 
|---|
| 164 |  | 
|---|
| 165 | /* ignore if not ours */ | 
|---|
| 166 | if (strcmp (dev, mydev)) | 
|---|
| 167 | return; | 
|---|
| 168 |  | 
|---|
| 169 | if (!strcmp (name, eqNPR.name)) { | 
|---|
| 170 | /* new equatorial target coords */ | 
|---|
| 171 | double newra = 0, newdec = 0; | 
|---|
| 172 | int i, nset; | 
|---|
| 173 |  | 
|---|
| 174 | /* Check connectSP, if it is idle, then return */ | 
|---|
| 175 | if (connectSP.s == IPS_IDLE) | 
|---|
| 176 | { | 
|---|
| 177 | eqNPR.s = IPS_IDLE; | 
|---|
| 178 | IDSetNumber(&eqNP, "Telescope is offline."); | 
|---|
| 179 | return; | 
|---|
| 180 | } | 
|---|
| 181 |  | 
|---|
| 182 | for (nset = i = 0; i < n; i++) | 
|---|
| 183 | { | 
|---|
| 184 | /* Find numbers with the passed names in the eqNP property */ | 
|---|
| 185 | INumber *eqp = IUFindNumber (&eqNPR, names[i]); | 
|---|
| 186 |  | 
|---|
| 187 | /* If the number found is Right ascension (eqN[0]) then process it */ | 
|---|
| 188 | if (eqp == &eqNR[0]) | 
|---|
| 189 | { | 
|---|
| 190 | newra = (values[i]); | 
|---|
| 191 | nset += newra >= 0 && newra <= 24; | 
|---|
| 192 | } | 
|---|
| 193 | /* Otherwise, if the number found is Declination (eqN[1]) then process it */ | 
|---|
| 194 | else if (eqp == &eqNR[1]) { | 
|---|
| 195 | newdec = (values[i]); | 
|---|
| 196 | nset += newdec >= -90 && newdec <= 90; | 
|---|
| 197 | } | 
|---|
| 198 | } /* end for */ | 
|---|
| 199 |  | 
|---|
| 200 | /* Did we process the two numbers? */ | 
|---|
| 201 | if (nset == 2) | 
|---|
| 202 | { | 
|---|
| 203 | char r[32], d[32]; | 
|---|
| 204 |  | 
|---|
| 205 | /* Set the mount state to BUSY */ | 
|---|
| 206 | eqNP.s = IPS_BUSY; | 
|---|
| 207 | eqNPR.s = IPS_BUSY; | 
|---|
| 208 |  | 
|---|
| 209 | /* Set the new target coordinates */ | 
|---|
| 210 | targetRA = newra; | 
|---|
| 211 | targetDEC = newdec; | 
|---|
| 212 |  | 
|---|
| 213 | /* Convert the numeric coordinates to a sexagesmal string (H:M:S) */ | 
|---|
| 214 | fs_sexa (r, targetRA, 2, 3600); | 
|---|
| 215 | fs_sexa (d, targetDEC, 3, 3600); | 
|---|
| 216 |  | 
|---|
| 217 | IDSetNumber(&eqNP, NULL); | 
|---|
| 218 | IDSetNumber(&eqNPR, "Moving to RA Dec %s %s", r, d); | 
|---|
| 219 | } | 
|---|
| 220 | /* We didn't process the two number correctly, report an error */ | 
|---|
| 221 | else | 
|---|
| 222 | { | 
|---|
| 223 | /* Set property state to ALERT */ | 
|---|
| 224 | eqNPR.s = IPS_ALERT; | 
|---|
| 225 |  | 
|---|
| 226 | IDSetNumber(&eqNP, "RA or Dec absent or bogus."); | 
|---|
| 227 | } | 
|---|
| 228 |  | 
|---|
| 229 | return; | 
|---|
| 230 | } | 
|---|
| 231 | } | 
|---|
| 232 |  | 
|---|
| 233 | /* client is sending us a new value for a Switch property */ | 
|---|
| 234 | void ISNewSwitch (const char *dev, const char *name, ISState *states, char *names[], int n) | 
|---|
| 235 | { | 
|---|
| 236 | ISwitch *sp; | 
|---|
| 237 |  | 
|---|
| 238 | mountInit(); | 
|---|
| 239 |  | 
|---|
| 240 | /* ignore if not ours */ | 
|---|
| 241 | if (strcmp (dev, mydev)) | 
|---|
| 242 | return; | 
|---|
| 243 |  | 
|---|
| 244 |  | 
|---|
| 245 | if (!strcmp(name, connectSP.name)) | 
|---|
| 246 | { | 
|---|
| 247 | /* We update switches. This is different from the way we used to update switches in tutorial 1. This is | 
|---|
| 248 | * to illustrate that there are several ways to update the switches. Here, we try to find the switch with names[0], | 
|---|
| 249 | * and if found, we update its state to states[0] and call connectTelescope(). We must call IUResetSwitches to erase any previous history */ | 
|---|
| 250 |  | 
|---|
| 251 | sp = IUFindSwitch (&connectSP, names[0]); | 
|---|
| 252 |  | 
|---|
| 253 | if (sp) | 
|---|
| 254 | { | 
|---|
| 255 | IUResetSwitch(&connectSP); | 
|---|
| 256 | sp->s = states[0]; | 
|---|
| 257 | connectTelescope(); | 
|---|
| 258 | } | 
|---|
| 259 | } | 
|---|
| 260 | else if (! strcmp(name, MovementNSSP.name)) { | 
|---|
| 261 | sp = IUFindSwitch (&MovementNSSP, names[0]); | 
|---|
| 262 | if (sp) { | 
|---|
| 263 | IUResetSwitch(&MovementNSSP); | 
|---|
| 264 | sp->s = states[0]; | 
|---|
| 265 | show_runtime(sp->s); | 
|---|
| 266 | IDSetSwitch (&MovementNSSP, "Toggle North/South."); | 
|---|
| 267 | } | 
|---|
| 268 | } | 
|---|
| 269 | else if (! strcmp(name, MovementWESP.name)) { | 
|---|
| 270 | sp = IUFindSwitch (&MovementWESP, names[0]); | 
|---|
| 271 | if (sp) { | 
|---|
| 272 | IUResetSwitch(&MovementWESP); | 
|---|
| 273 | sp->s = states[0]; | 
|---|
| 274 | show_runtime(sp->s); | 
|---|
| 275 | IDSetSwitch (&MovementWESP, "Toggle West/East."); | 
|---|
| 276 | } | 
|---|
| 277 | } | 
|---|
| 278 |  | 
|---|
| 279 |  | 
|---|
| 280 | } | 
|---|
| 281 |  | 
|---|
| 282 | /* update the "mount" over time */ | 
|---|
| 283 | void mountSim (void *p) | 
|---|
| 284 | { | 
|---|
| 285 | static struct timeval ltv; | 
|---|
| 286 | struct timeval tv; | 
|---|
| 287 | double dt, da, dx; | 
|---|
| 288 | int nlocked; | 
|---|
| 289 |  | 
|---|
| 290 | /* If telescope is not on, do not simulate */ | 
|---|
| 291 | if (connectSP.s == IPS_IDLE) | 
|---|
| 292 | { | 
|---|
| 293 | IEAddTimer (POLLMS, mountSim, NULL); | 
|---|
| 294 | return; | 
|---|
| 295 | } | 
|---|
| 296 |  | 
|---|
| 297 | /* update elapsed time since last poll, don't presume exactly POLLMS */ | 
|---|
| 298 | gettimeofday (&tv, NULL); | 
|---|
| 299 |  | 
|---|
| 300 | if (ltv.tv_sec == 0 && ltv.tv_usec == 0) | 
|---|
| 301 | ltv = tv; | 
|---|
| 302 |  | 
|---|
| 303 | dt = tv.tv_sec - ltv.tv_sec + (tv.tv_usec - ltv.tv_usec)/1e6; | 
|---|
| 304 | ltv = tv; | 
|---|
| 305 | da = SLEWRATE*dt; | 
|---|
| 306 |  | 
|---|
| 307 | /* Process per current state. We check the state of EQUATORIAL_EOD_COORDS_REQUEST and act acoordingly */ | 
|---|
| 308 | switch (eqNPR.s) | 
|---|
| 309 | { | 
|---|
| 310 |  | 
|---|
| 311 | /* #1 State is idle, update telesocpe at sidereal rate */ | 
|---|
| 312 | case IPS_IDLE: | 
|---|
| 313 | /* RA moves at sidereal, Dec stands still */ | 
|---|
| 314 | currentRA += (SIDRATE*dt/15.); | 
|---|
| 315 | IDSetNumber(&eqNP, NULL); | 
|---|
| 316 | break; | 
|---|
| 317 |  | 
|---|
| 318 | case IPS_BUSY: | 
|---|
| 319 | /* slewing - nail it when both within one pulse @ SLEWRATE */ | 
|---|
| 320 | nlocked = 0; | 
|---|
| 321 |  | 
|---|
| 322 | dx = targetRA - currentRA; | 
|---|
| 323 |  | 
|---|
| 324 | if (fabs(dx) <= da) | 
|---|
| 325 | { | 
|---|
| 326 | currentRA = targetRA; | 
|---|
| 327 | nlocked++; | 
|---|
| 328 | } | 
|---|
| 329 | else if (dx > 0) | 
|---|
| 330 | currentRA += da/15.; | 
|---|
| 331 | else | 
|---|
| 332 | currentRA -= da/15.; | 
|---|
| 333 |  | 
|---|
| 334 |  | 
|---|
| 335 | dx = targetDEC - currentDec; | 
|---|
| 336 | if (fabs(dx) <= da) | 
|---|
| 337 | { | 
|---|
| 338 | currentDec = targetDEC; | 
|---|
| 339 | nlocked++; | 
|---|
| 340 | } | 
|---|
| 341 | else if (dx > 0) | 
|---|
| 342 | currentDec += da; | 
|---|
| 343 | else | 
|---|
| 344 | currentDec -= da; | 
|---|
| 345 |  | 
|---|
| 346 | if (nlocked == 2) | 
|---|
| 347 | { | 
|---|
| 348 | eqNP.s = IPS_OK; | 
|---|
| 349 | eqNPR.s = IPS_OK; | 
|---|
| 350 | IDSetNumber(&eqNP, NULL); | 
|---|
| 351 | IDSetNumber(&eqNPR, "Now tracking"); | 
|---|
| 352 | } else | 
|---|
| 353 | IDSetNumber(&eqNP, NULL); | 
|---|
| 354 |  | 
|---|
| 355 | break; | 
|---|
| 356 |  | 
|---|
| 357 | case IPS_OK: | 
|---|
| 358 | /* tracking */ | 
|---|
| 359 | IDSetNumber(&eqNP, NULL); | 
|---|
| 360 | break; | 
|---|
| 361 |  | 
|---|
| 362 | case IPS_ALERT: | 
|---|
| 363 | break; | 
|---|
| 364 | } | 
|---|
| 365 |  | 
|---|
| 366 | /* again */ | 
|---|
| 367 | IEAddTimer (POLLMS, mountSim, NULL); | 
|---|
| 368 | } | 
|---|
| 369 |  | 
|---|
| 370 | /* Note that we must define ISNewBLOB and ISSnoopDevice even if we don't use them, otherwise, the driver will NOT compile */ | 
|---|
| 371 | void ISNewBLOB (const char *dev, const char *name, int sizes[], int blobsizes[], char *blobs[], char *formats[], char *names[], int n) {} | 
|---|
| 372 | void ISSnoopDevice (XMLEle *root) {} | 
|---|
| 373 |  | 
|---|
| 374 |  | 
|---|
| 375 | static void connectTelescope () | 
|---|
| 376 | { | 
|---|
| 377 | if (connectS[0].s == ISS_ON) | 
|---|
| 378 | { | 
|---|
| 379 | connectSP.s   = IPS_OK; | 
|---|
| 380 | IDSetSwitch (&connectSP, "Telescope is connected."); | 
|---|
| 381 | } | 
|---|
| 382 | else | 
|---|
| 383 | { | 
|---|
| 384 | connectSP.s   = IPS_IDLE; | 
|---|
| 385 | IDSetSwitch (&connectSP, "Telescope is disconnected."); | 
|---|
| 386 | } | 
|---|
| 387 | } | 
|---|