From ce53f1f45d32a68656d1cab94726dbea3dfa9aa3 Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <w.l.netto@medisin.uio.no>
Date: Fri, 23 Dec 2022 13:11:55 +0100
Subject: [PATCH 1/3] Translated ind2subv() (#3)

---
 R/ind2subv.R | 67 ++++++++++++++++++++++++++--------------------------
 1 file changed, 33 insertions(+), 34 deletions(-)

diff --git a/R/ind2subv.R b/R/ind2subv.R
index 6fef574..92b4943 100644
--- a/R/ind2subv.R
+++ b/R/ind2subv.R
@@ -1,38 +1,37 @@
-ind2subv <- function(siz, ndx) stop("Needs translation")
-# function sub = ind2subv(siz, ndx)
-# % IND2SUBV Like the built-in ind2sub, but returns the answer as a row vector.
-# % sub = ind2subv(siz, ndx)
-# %
-# % siz and ndx can be row or column vectors.
-# % sub will be of size length(ndx) * length(siz).
-# %
-# % Example
-# % ind2subv([2 2 2], 1:8) returns
-# %  [1 1 1
-# %   2 1 1
-# %   ...
-# %   2 2 2]
-# % That is, the leftmost digit toggle fastest.
-# %
-# % See also SUBV2IND
+ind2subv <- function(siz, ndx) {
+  #  IND2SUBV Like the built - in ind2sub, but returns the answer as a row vector.
+  #  sub <- ind2subv(siz, ndx)
+  #  siz and ndx can be row or column vectors.
+  #  sub will be of size length(ndx) * length(siz).
+  #  Example
+  #  ind2subv([2 2 2], 1:8) returns
+  #   [1 1 1
+  #    2 1 1
+  #     ...
+  #    2 2 2]
+  #  That is, the leftmost digit toggle fastest.
+  #
+  #  See also SUBV2IND
 
-# n = length(siz);
+  n <- length(siz)
 
-# if n==0
-#   sub = ndx;
-#   return;
-# end
+  if (n == 0) {
+    sub <- ndx
+    return(sub)
+  }
 
-# if all(siz==2)
-#   sub = dec2bitv(ndx-1, n);
-#   sub = sub(:,n:-1:1)+1;
-#   return;
-# end
+  if (all(siz == 2)) {
+    sub <- dec2bitv(ndx - 1, n)
+    sub <- sub[, seq(n, 1, - 1)] + 1
+    return(sub)
+  }
 
-# cp = [1 cumprod(siz(:)')];
-# ndx = ndx(:) - 1;
-# sub = zeros(length(ndx), n);
-# for i = n:-1:1 % i'th digit
-#   sub(:,i) = floor(ndx/cp(i))+1;
-#   ndx = rem(ndx,cp(i));
-# end
+  cp <- c(1, cumprod(t(siz[])))
+  ndx <- ndx[] - 1
+  sub <- zeros(length(ndx), n)
+  for (i in seq(n, 1, -1)) {# i'th digit
+    sub[, i] <- floor(ndx / cp[i]) + 1
+    ndx <- ndx %% cp(i)
+  }
+  return(sub)
+}

From 98cde6b23e67a881faa6077111675891c0e1caff Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <w.l.netto@medisin.uio.no>
Date: Fri, 23 Dec 2022 13:41:50 +0100
Subject: [PATCH 2/3] Translated dec2bitv() (#3)

---
 R/dec2bitv.R      | 16 ++++++++++++++++
 R/rBAPS-package.R |  2 +-
 2 files changed, 17 insertions(+), 1 deletion(-)
 create mode 100644 R/dec2bitv.R

diff --git a/R/dec2bitv.R b/R/dec2bitv.R
new file mode 100644
index 0000000..2124936
--- /dev/null
+++ b/R/dec2bitv.R
@@ -0,0 +1,16 @@
+dec2bitv <- function(d, n) {
+  # DEC2BITV Convert a decimal integer to a bit vector.
+  # bits <- dec2bitv(d, n) is just like the built - in dec2bin, except the answer is a vector, not a as.character.
+  # n is an optional minimum length on the bit vector.
+  # If d is a vector,  each row of the output array will be a bit vector.
+
+  if (nargin() < 2) {
+    n <- 1 # Need at least one digit even for 0.
+  }
+  d <- d[]
+
+  f <- e <- NA
+  c(f, e) <- log2(max(d)) # How many digits do we need to represent the numbers?
+    bits <- floor(d * 2 ^ (seq(1 - max(n, e), 0))) %% 2
+  return(bits)
+}
diff --git a/R/rBAPS-package.R b/R/rBAPS-package.R
index 25e3934..810a2f1 100644
--- a/R/rBAPS-package.R
+++ b/R/rBAPS-package.R
@@ -10,5 +10,5 @@
 #' size sortrows squeeze strcmp times zeros disp
 #' @importFrom stats runif
 #' @importFrom zeallot %<-%
-#' @importFrom matlab2r nargin
+#' @importFrom matlab2r nargin log2
 NULL

From 8dc8a198a398999e0f56db466a3f0b42f37dec84 Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <w.l.netto@medisin.uio.no>
Date: Fri, 23 Dec 2022 13:42:35 +0100
Subject: [PATCH 3/3] Increment version number to 0.0.0.9014

---
 DESCRIPTION | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/DESCRIPTION b/DESCRIPTION
index 6146e4c..38cca8d 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: rBAPS
 Title: Bayesian Analysis of Population Structure
-Version: 0.0.0.9013
+Version: 0.0.0.9014
 Date: 2020-11-09
 Authors@R:
     c(