/* Convert radian argument to DMS ascii format */

#include <math.h>
#include <stddef.h>
#include <stdio.h>
#include <string.h>

#include "proj.h"
#include "proj_internal.h"

/*
** RES is fractional second figures
** RES60 = 60 * RES
** CONV = 180 * 3600 * RES / PI (radians to RES seconds)
*/
static double RES = 1000., RES60 = 60000., CONV = 206264806.24709635516;
static char format[50] = "%dd%d'%.3f\"%c";
static int dolong = 0;
void set_rtodms(int fract, int con_w) {
    int i;

    if (fract >= 0 && fract < 9) {
        RES = 1.;
        /* following not very elegant, but used infrequently */
        for (i = 0; i < fract; ++i)
            RES *= 10.;
        RES60 = RES * 60.;
        CONV = 180. * 3600. * RES / M_PI;
        if (!con_w)
            (void)snprintf(format, sizeof(format), "%%dd%%d'%%.%df\"%%c",
                           fract);
        else
            (void)snprintf(format, sizeof(format), "%%dd%%02d'%%0%d.%df\"%%c",
                           fract + 2 + (fract ? 1 : 0), fract);
        dolong = con_w;
    }
}
char *rtodms(char *s, size_t sizeof_s, double r, int pos, int neg) {
    int deg, min, sign;
    char *ss = s;
    double sec;
    size_t sizeof_ss = sizeof_s;

    if (r < 0) {
        r = -r;
        if (!pos) {
            if (sizeof_s == 1) {
                *s = 0;
                return s;
            }
            sizeof_ss--;
            *ss++ = '-';
            sign = 0;
        } else
            sign = neg;
    } else
        sign = pos;
    r = floor(r * CONV + .5);
    sec = fmod(r / RES, 60.);
    r = floor(r / RES60);
    min = (int)fmod(r, 60.);
    r = floor(r / 60.);
    deg = (int)r;

    if (dolong)
        (void)snprintf(ss, sizeof_ss, format, deg, min, sec, sign);
    else if (sec != 0.0) {
        char *p, *q;
        /* double prime + pos/neg suffix (if included) + NUL */
        size_t suffix_len = sign ? 3 : 2;

        (void)snprintf(ss, sizeof_ss, format, deg, min, sec, sign);
        /* Replace potential decimal comma by decimal point for non C locale */
        for (p = ss; *p != '\0'; ++p) {
            if (*p == ',') {
                *p = '.';
                break;
            }
        }
        if (suffix_len > strlen(ss))
            return s;
        for (q = p = ss + strlen(ss) - suffix_len; *p == '0'; --p)
            ;
        if (*p != '.')
            ++p;
        if (++q != p)
            (void)memmove(p, q, suffix_len);
    } else if (min)
        (void)snprintf(ss, sizeof_ss, "%dd%d'%c", deg, min, sign);
    else
        (void)snprintf(ss, sizeof_ss, "%dd%c", deg, sign);
    return s;
}