add transform

simplex_transforms/stan/transforms/Exponential/Exponential_functions.stan

@@ -0,0 +1,21 @@
+vector exponential_simplex_constrain_lp(vector y) {
+  int N = rows(y);
+  vector[N] x_pos = y - log1p_exp(-y);
+  real r = log_sum_exp(x_pos);
+  vector[N] x = exp(x_pos - r);
+  target += x_pos;
+  target += std_normal_lpdf(r - log(N));
+  target += sum(y) - N * r;
+  return x;
+}
+
+vector exponential_log_simplex_constrain_lp(vector y) {
+  int N = rows(y);
+  vector[N] x = y - log1p_exp(-y);
+  real r = log_sum_exp(x);
+  vector[N] log_x = x - r;
+  target += x;
+  target += std_normal_lpdf(r - log(N));
+  target += log_x[N];
+  return log_x;
+}

simplex_transforms/stan/transforms/Exponential/Exponential_parameters_log_simplex.stan

@@ -0,0 +1,7 @@
+parameters {
+  vector[N] y;
+}
+transformed parameters {
+  vector<upper=0>[N] log_x = exponential_log_simplex_constrain_lp(y);
+  simplex[N] x = exp(log_x);
+}

simplex_transforms/stan/transforms/Exponential/Exponential_parameters_simplex.stan

@@ -0,0 +1,6 @@
+parameters {
+  vector[N] y;
+}
+transformed parameters {
+  simplex[N] x = exponential_simplex_constrain_lp(y);
+}