highpmstarlist.cpp

Go to the documentation of this file.

/***************************************************************************
                        highpmstarlist.cpp  -  K Desktop Planetarium
                             -------------------
    begin                : 2007-08-06
    copyright            : (C) 2007 by James B. Bowlin
    email                : bowlin@mindspring.com
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#include "highpmstarlist.h"
#include "skyobjects/starobject.h"
#include "kstarsdatetime.h"
#include "skymesh.h"


typedef struct HighPMStar
{
    HighPMStar( Trixel t, StarObject* s ) : trixel(t), star(s) {}
    Trixel      trixel;
    StarObject  *star;

} HighPMStar;


HighPMStarList::HighPMStarList( double threshold )
        : m_reindexNum(J2000), m_threshold(threshold), m_maxPM(0.0)
{
    m_skyMesh = SkyMesh::Instance();
}

HighPMStarList::~HighPMStarList()
{
    for (int i = 0; i < size(); i++ ) {
        delete m_stars[ i ];
    }
}

bool HighPMStarList::append( Trixel trixel, StarObject* star, double pm )
{
    if ( pm < m_threshold ) return false;

    if( trixel >= (unsigned int)m_skyMesh->size() )
        kDebug() << "### Trixel ID out of range for the Mesh currently in use!" << trixel;

    m_stars.append( new HighPMStar( trixel, star ) );
    if ( m_maxPM >= pm ) return true;

    m_maxPM = pm;
    m_reindexInterval = StarObject::reindexInterval( pm );

    return true;
}

void HighPMStarList::setIndexTime( KSNumbers *num )
{
    m_reindexNum = KSNumbers( *num );
}

void HighPMStarList::reindex( KSNumbers *num, StarIndex* starIndex )
{
    if ( fabs( num->julianCenturies() -
               m_reindexNum.julianCenturies() ) < m_reindexInterval ) return;

    m_reindexNum = KSNumbers( *num );
    m_skyMesh->setKSNumbers( num );

    int cnt(0);

    for ( int i = 0; i < m_stars.size(); i++ ) {
        HighPMStar* HPStar = m_stars.at( i );
        Trixel trixel = m_skyMesh->indexStar( HPStar->star );
        if ( trixel == HPStar->trixel ) continue;
        cnt++;
        StarObject* star = HPStar->star;

        // out with the old ...
        if( HPStar->trixel >= (unsigned int)m_skyMesh->size() ) {
            kDebug() << "### Expect an Index out-of-range error. star->trixel =" << HPStar->trixel;
        }
            
        StarList *old = starIndex->at( HPStar->trixel );
        old->removeAt( old->indexOf( star ) );
            
        float mag = star->mag();
        //printf("\n mag = %4.2f  trixel %d -> %d\n", mag, HPStar->trixel, trixel );

        // in with the new ...
        HPStar->trixel = trixel;
        if( trixel >= (unsigned int)m_skyMesh->size() )
            kDebug() << "### Expect an Index out-of-range error. trixel =" << trixel;
                
        StarList *list = starIndex->at( trixel );
        int j;
        for ( j = 0; j < list->size(); j++ ) {
            if ( list->at(j)->mag() < mag ) continue;
            list->insert( j, star );
            break;
        }
        if ( j == list->size() ) list->append( star );

        //for ( j = 0; j < list->size(); j++ ) {
        //    printf("    %4.2f\n", list->at(j)->mag() );
        //}
    }
    //printf("Re-indexed %d stars at interval %6.1f\n", cnt, 100.0 * m_reindexInterval );
}


void HighPMStarList::stats()
{
    printf("\n");
    printf("maxPM: %6.1f  threshold %5.1f\n", m_maxPM, m_threshold );
    printf("stars: %d\n", size() );
    printf("Update Interval: %6.1f years\n",  100.0 * m_reindexInterval );
    printf("Last Update: %6.1f\n", 2000.0 + 100.0 * (int) m_reindexNum.julianCenturies() );
}