ealdent-lda-ruby 0.1.1

data/utils.c

@@ -0,0 +1,111 @@
+#include "utils.h"
+
+/*
+ * given log(a) and log(b), return log(a + b)
+ *
+ */
+
+double log_sum(double log_a, double log_b)
+{
+  double v;
+
+  if (log_a < log_b)
+  {
+      v = log_b+log(1 + exp(log_a-log_b));
+  }
+  else
+  {
+      v = log_a+log(1 + exp(log_b-log_a));
+  }
+  return(v);
+}
+
+ /**
+   * Proc to calculate the value of the trigamma, the second
+   * derivative of the loggamma function. Accepts positive matrices.
+   * From Abromowitz and Stegun.  Uses formulas 6.4.11 and 6.4.12 with
+   * recurrence formula 6.4.6.  Each requires workspace at least 5
+   * times the size of X.
+   *
+   **/
+
+double trigamma(double x)
+{
+    double p;
+    int i;
+
+    x=x+6;
+    p=1/(x*x);
+    p=(((((0.075757575757576*p-0.033333333333333)*p+0.0238095238095238)
+         *p-0.033333333333333)*p+0.166666666666667)*p+1)/x+0.5*p;
+    for (i=0; i<6 ;i++)
+    {
+        x=x-1;
+        p=1/(x*x)+p;
+    }
+    return(p);
+}
+
+
+/*
+ * taylor approximation of first derivative of the log gamma function
+ *
+ */
+
+double digamma(double x)
+{
+    double p;
+    x=x+6;
+    p=1/(x*x);
+    p=(((0.004166666666667*p-0.003968253986254)*p+
+	0.008333333333333)*p-0.083333333333333)*p;
+    p=p+log(x)-0.5/x-1/(x-1)-1/(x-2)-1/(x-3)-1/(x-4)-1/(x-5)-1/(x-6);
+    return p;
+}
+
+
+double log_gamma(double x)
+{
+     double z=1/(x*x);
+
+    x=x+6;
+    z=(((-0.000595238095238*z+0.000793650793651)
+	*z-0.002777777777778)*z+0.083333333333333)/x;
+    z=(x-0.5)*log(x)-x+0.918938533204673+z-log(x-1)-
+	log(x-2)-log(x-3)-log(x-4)-log(x-5)-log(x-6);
+    return z;
+}
+
+
+
+/*
+ * make directory
+ *
+ */
+
+void make_directory(char* name)
+{
+    mkdir(name, S_IRUSR|S_IWUSR|S_IXUSR);
+}
+
+
+/*
+ * argmax
+ *
+ */
+
+int argmax(double* x, int n)
+{
+    int i;
+    double max = x[0];
+    int argmax = 0;
+    for (i = 1; i < n; i++)
+    {
+        if (x[i] > max)
+        {
+            max = x[i];
+            argmax = i;
+        }
+    }
+    return(argmax);
+}

data/utils.h

@@ -0,0 +1,18 @@
+#ifndef UTILS_H
+#define UTILS_H
+
+#include <stdio.h>
+#include <math.h>
+#include <float.h>
+#include <stdlib.h>
+#include <sys/stat.h>
+#include <sys/types.h>
+
+double log_sum(double log_a, double log_b);
+double trigamma(double x);
+double digamma(double x);
+double log_gamma(double x);
+void make_directory(char* name);
+int argmax(double* x, int n);
+
+#endif

metadata

@@ -0,0 +1,70 @@
+--- !ruby/object:Gem::Specification 
+name: ealdent-lda-ruby
+version: !ruby/object:Gem::Version 
+  version: 0.1.1
+platform: ruby
+authors: 
+- Jason M. Adams
+- David M. Blei
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2008-11-14 00:00:00 -08:00
+default_executable: 
+dependencies: []
+
+description: 
+email: jasonmadams@gmail.com
+executables: []
+
+extensions: []
+
+extra_rdoc_files: []
+
+files: 
+- README
+- license.txt
+- cokus.c
+- cokus.h
+- extconf.rb
+- lda-alpha.c
+- lda-alpha.h
+- lda-data.c
+- lda-data.h
+- lda-inference.c
+- lda-inference.h
+- lda-model.c
+- lda-model.h
+- lda.h
+- lda.rb
+- utils.c
+- utils.h
+has_rdoc: true
+homepage: http://github.com/ealdent/lda-ruby
+post_install_message: 
+rdoc_options: []
+
+require_paths: 
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+requirements: []
+
+rubyforge_project: 
+rubygems_version: 1.2.0
+signing_key: 
+specification_version: 2
+summary: Ruby port of Latent Dirichlet Allocation by David M. Blei.
+test_files: []
+