planck.cc

Go to the documentation of this file.

#include "planck.h"

using namespace blitz;
using namespace FullPhysics;

// Inform C about the Fortran function we will use
extern "C" {
    void planckfunction(const double* wavenumber, const double* temperature, double* bbfunc, int* smallv, bool* fail, const int* message_len, char* message);
    void planckfunction_plus(const double* wavenumber, const double* temperature, double* bbfunc, double* deriv_bbfunc, int* smallv, bool* fail, const int* message_len, char* message);
}

double FullPhysics::planck(double wn, double temperature)
{
    double bbfunc;
    int smallv;
    bool fail;
    const int message_len = 50;
    char message[message_len];

    planckfunction(&wn, &temperature, &bbfunc, &smallv, &fail, &message_len, message);

    if (fail) {
        Exception err;
        err << "plankfunction_plus failed with error message: " << message;
        throw err;
    }

    return bbfunc;
}

AutoDerivative<double> FullPhysics::planck(double wn, double temperature, Array<double, 1>& gradient)
{
    double bbfunc;
    double deriv_bbfunc;
    int smallv;
    bool fail;
    const int message_len = 50;
    char message[message_len];

    planckfunction_plus(&wn, &temperature, &bbfunc, &deriv_bbfunc, &smallv, &fail, &message_len, message);

    if (fail) {
        Exception err;
        err << "plankfunction_plus failed with error message: " << message;
        throw err;
    }

    Array<double, 1> out_gradient(gradient.rows());
    out_gradient = gradient * deriv_bbfunc;

    AutoDerivative<double> bbody(bbfunc, out_gradient);

    return bbody;
}