dnorm.hpp

Go to the documentation of this file.

// Copyright (C) 2013-2015 Kasper Kristensen
// License: GPL-2

template<class Type>
Type dnorm(Type x, Type mean, Type sd, int give_log=0)
{
  Type resid = (x - mean) / sd;
  Type logans = Type(-log(sqrt(2*M_PI))) - log(sd) - Type(.5) * resid * resid;
  if(give_log) return logans; else return exp(logans);
}
VECTORIZE4_ttti(dnorm)