/*****************************************************************************************
 *                                                                                       *
 * OpenSpace                                                                             *
 *                                                                                       *
 * Copyright (c) 2014-2020                                                               *
 *                                                                                       *
 * Permission is hereby granted, free of charge, to any person obtaining a copy of this  *
 * software and associated documentation files (the "Software"), to deal in the Software *
 * without restriction, including without limitation the rights to use, copy, modify,    *
 * merge, publish, distribute, sublicense, and/or sell copies of the Software, and to    *
 * permit persons to whom the Software is furnished to do so, subject to the following   *
 * conditions:                                                                           *
 *                                                                                       *
 * The above copyright notice and this permission notice shall be included in all copies *
 * or substantial portions of the Software.                                              *
 *                                                                                       *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,   *
 * INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A         *
 * PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT    *
 * HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF  *
 * CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE  *
 * OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.                                         *
 ****************************************************************************************/

#version __CONTEXT__

#include "floatoperations.glsl"

layout (location = 0) out vec4 outColor;

// Keep in sync with gaiaoptions.h:RenderOption enum
const int RENDEROPTION_STATIC = 0;
const int RENDEROPTION_COLOR = 1;
const int RENDEROPTION_MOTION = 2; 
const float ONE_PARSEC = 3.08567758e16; // 1 Parsec
const float FLT_MAX  = 3.402823466e38; // Max float constant in GLSL
const float LUM_LOWER_CAP = 0.01;

in vec2 ge_brightness;
in vec4 ge_gPosition;
in vec2 texCoord;
in float ge_starDistFromSun;
in float ge_cameraDistFromSun;
in float ge_observedDist;

uniform sampler2D psfTexture;
uniform sampler1D colorTexture;
uniform float luminosityMultiplier;
uniform float sharpness;
uniform int renderOption;

vec3 color2rgb(float color) {
    // BV is [-0.4, 2.0]
    float st = (color + 0.4) / (2.0 + 0.4);

    // Bp-Rp[-2.0, 6.5], Bp-G[-2.1, 5.0], G-Rp[-1.0, 3.0]
    //float st = (color + 1.0) / (5.0 + 1.0);

    return texture(colorTexture, st).rgb;
}

void main() {

    // Assume all stars has equal luminosity as the Sun when no magnitude is loaded.
    float luminosity = 0.05;
    vec3 color = vec3(luminosity);
    float ratioMultiplier = 0.05;

    vec4 textureColor = texture(psfTexture, texCoord);
    textureColor.a = pow(textureColor.a, sharpness);
    if (textureColor.a < 0.001) {
        discard;
    }

    // Calculate the color and luminosity if we have the magnitude and B-V color.
    if ( renderOption != RENDEROPTION_STATIC ) {
        color = color2rgb(ge_brightness.y);
        ratioMultiplier = 0.5;

        // Absolute magnitude is brightness a star would have at 10 pc away.
        float absoluteMagnitude = ge_brightness.x;

        // From formula: MagSun - MagStar = 2.5*log(LumStar / LumSun), it gives that:
        // LumStar = 10^(1.89 - 0.4*Magstar) , if LumSun = 1 and MagSun = 4.72
        luminosity = pow(10.0, 1.89 - 0.4 * absoluteMagnitude);

        // If luminosity is really really small then set it to a static low number.
        if (luminosity < LUM_LOWER_CAP) {
            luminosity = LUM_LOWER_CAP;
        }
    }

    // Luminosity decrease by {squared} distance [measured in Pc].
    float observedDistance = ge_observedDist / ONE_PARSEC;
    luminosity /= pow(observedDistance, 2.0);

    // Multiply our color with the luminosity as well as a user-controlled property.
    color *= luminosity * pow(luminosityMultiplier, 3.0);

    // Decrease contributing brightness for stars in central cluster.
    if ( ge_cameraDistFromSun > ge_starDistFromSun ) {
        float ratio = ge_starDistFromSun / ge_cameraDistFromSun;
        //color *= ratio * ratioMultiplier;
    }

    // Use truncating tonemapping here so we don't overexposure individual stars.
    //color = 1.0 - 1.0 * exp(-5.0 * color.rgb);
    float maxVal = max(max(color.r, color.g), color.b);
    if (maxVal > 1.0) {
        color /= maxVal;
    }

    outColor = vec4(color, textureColor.a);
}
