diff git a/DESCRIPTION b/DESCRIPTION
index 3a5505b2991bb40b0659f0afdeeb289aba31a4f8..3f6fffe11e7d7162bcb355baf17678cd9f992ddd 100644
 a/DESCRIPTION
+++ b/DESCRIPTION
@@ 1,6 +1,6 @@
Package: multiApply
Title: Apply Functions to Multiple Multidimensional Arrays or Vectors
Version: 2.0.1
+Version: 2.1.0
Authors@R: c(
person("BSCCNS", role = c("aut", "cph")),
person("Nicolau", "Manubens", , "nicolau.manubens@bsc.es", role = "aut"),
diff git a/man/Apply.Rd b/man/Apply.Rd
index b93c173352638ae4b7e37aa65e42b1cc72f41b63..a368553101b7307567bf436584696e72824ef6d1 100644
 a/man/Apply.Rd
+++ b/man/Apply.Rd
@@ 9,11 +9,11 @@ Apply(data, target_dims = NULL, fun, ..., output_dims = NULL,
guess_dim_names = TRUE, ncores = NULL, split_factor = 1)
}
\arguments{
\item{data}{One or a list of numeric object (vector, matrix or array). They must be in the same order as expected by the function provided in the parameter 'fun'. The dimensions do not necessarily have to be ordered. If the 'target_dims' require a different order than the provided, \code{Apply} will automatically reorder the dimensions as needed.}
+\item{data}{One or a list of vectors, matrices or arrays. They must be in the same order as expected by the function provided in the parameter 'fun'. The dimensions do not necessarily have to be ordered. If the 'target_dims' require a different order than the provided, \code{Apply} will automatically reorder the dimensions as needed.}
\item{target_dims}{One or a list of vectors (or NULLs) containing the dimensions to be input into fun for each of the objects in the data. If a single vector of target dimensions is specified and multiple inputs are provided in 'data, then the single set of target dimensions is reused for all of the inputs. These vectors can contain either integers specifying the position of the dimensions, or character strings corresponding to the dimension names. This parameter is mandatory if 'margins' are not specified. If both 'margins' and 'target_dims' are specified, 'margins' takes priority.}
\item{fun}{Function to be applied to the arrays. Must receive as many inputs as provided in 'data', each with as many dimensions as specified in 'target_dims' or as the total number of dimensions in 'data' minus the ones specified in 'margins'. The function can receive other additional fixed parameters (see parameter '...' of \code{Apply}). The function can return one or a list of numeric vectors or multidimensional arrays, optionally with dimension names which will be propagated to the final result. The returned list can optionally be named, with a name for each output, which will be propagated to the resulting array. The function can optionally be provided with the attributes 'target_dims' and 'output_dims'. In that case, the corresponding parameters of \code{Apply} do not need to be provided. The function can expect named dimensions for each of its inputs, in the same order as specified in 'target_dims' or, if no 'target_dims' have been provided, in the same order as provided in 'data'.}
+\item{fun}{Function to be applied to the arrays. Must receive as many inputs as provided in 'data', each with as many dimensions as specified in 'target_dims' or as the total number of dimensions in 'data' minus the ones specified in 'margins'. The function can receive other additional fixed parameters (see parameter '...' of \code{Apply}). The function can return one or a list of vectors or multidimensional arrays, optionally with dimension names which will be propagated to the final result. The returned list can optionally be named, with a name for each output, which will be propagated to the resulting array. The function can optionally be provided with the attributes 'target_dims' and 'output_dims'. In that case, the corresponding parameters of \code{Apply} do not need to be provided. The function can expect named dimensions for each of its inputs, in the same order as specified in 'target_dims' or, if no 'target_dims' have been provided, in the same order as provided in 'data'. The function can access the variable \code{.margin_indices}, a named numeric vector that provides the indices of the current iteration over the margins, as well as any other variables specified in the parameter \code{extra_info} or input attributes specified in the parameter \code{use_attributes}.}
\item{...}{Additional fixed arguments expected by the function provided in the parameter 'fun'.}
@@ 35,7 +35,7 @@ Apply(data, target_dims = NULL, fun, ..., output_dims = NULL,
List of arrays or matrices or vectors resulting from applying 'fun' to 'data'.
}
\description{
This function efficiently applies a given function, which takes N vectors or multidimensional arrays as inputs (which may have different numbers of dimensions and dimension lengths), and applies it to a list of N vectors or multidimensional arrays with at least as many dimensions as expected by the given function. The user can specify which dimensions of each array the function is to be applied over with the \code{margins} or \code{target_dims} parameters. The function to be applied can receive other helper parameters and return any number of numeric vectors or multidimensional arrays. The target dimensions or margins can be specified by their names, as long as the inputs are provided with dimension names (recommended). This function can also use multicore in a transparent way if requested via the \code{ncores} parameter.\cr\cr The following steps help to understand how \code{Apply} works:\cr\cr  The function receives N arrays with Dn dimensions each.\cr  The user specifies, for each of the arrays, which of its dimensions are 'target' dimensions (dimensions which the function provided in 'fun' operates with) and which are 'margins' (dimensions to be looped over).\cr  \code{Apply} will generate an array with as many dimensions as margins in all of the input arrays. If a margin is repeated across different inputs, it will appear only once in the resulting array.\cr  For each element of this resulting array, the function provided in the parameter'fun' is applied to the corresponding subarrays in 'data'.\cr  If the function returns a vector or a multidimensional array, the additional dimensions will be prepended to the resulting array (in leftmost positions).\cr  If the provided function returns more than one vector or array, the process above is carried out for each of the outputs, resulting in a list with multiple arrays, each with the combination of all target dimensions (at the rightmost positions) and resulting dimensions (at the leftmost positions).
+This function efficiently applies a given function, which takes N vectors or multidimensional arrays as inputs (which may have different numbers of dimensions and dimension lengths), and applies it to a list of N vectors or multidimensional arrays with at least as many dimensions as expected by the given function. The user can specify which dimensions of each array the function is to be applied over with the \code{margins} or \code{target_dims} parameters. The function to be applied can receive other helper parameters and return any number of vectors or multidimensional arrays. The target dimensions or margins can be specified by their names, as long as the inputs are provided with dimension names (recommended). This function can also use multicore in a transparent way if requested via the \code{ncores} parameter.\cr\cr The following steps help to understand how \code{Apply} works:\cr\cr  The function receives N arrays with Dn dimensions each.\cr  The user specifies, for each of the arrays, which of its dimensions are 'target' dimensions (dimensions which the function provided in 'fun' operates with) and which are 'margins' (dimensions to be looped over).\cr  \code{Apply} will generate an array with as many dimensions as margins in all of the input arrays. If a margin is repeated across different inputs, it will appear only once in the resulting array.\cr  For each element of this resulting array, the function provided in the parameter'fun' is applied to the corresponding subarrays in 'data'.\cr  If the function returns a vector or a multidimensional array, the additional dimensions will be prepended to the resulting array (in leftmost positions).\cr  If the provided function returns more than one vector or array, the process above is carried out for each of the outputs, resulting in a list with multiple arrays, each with the combination of all target dimensions (at the rightmost positions) and resulting dimensions (at the leftmost positions).
}
\details{
When using a single object as input, Apply is almost identical to the apply function (as fast or slightly slower in some cases; with equal or improved smaller memory footprint).
diff git a/multiApplymanual.pdf b/multiApplymanual.pdf
index 8fabe141a52cc70c0e91997c670591579bd8489b..f5ae17b78b4c01d8c7637076bcd2c2412a9db99d 100644
Binary files a/multiApplymanual.pdf and b/multiApplymanual.pdf differ
diff git a/tests/testthat/testusecases.R b/tests/testthat/testusecases.R
index 6dd8d87a2605b5f2306e23f747ca8d6b28757e0f..b2533c83f2898adc93233eb9553a4d2291ca8394 100644
 a/tests/testthat/testusecases.R
+++ b/tests/testthat/testusecases.R
@@ 1346,21 +1346,15 @@ test_that("Dimension names are propagated correctly.", {
a < array(1:prod(1:6), dim = c(a = 1, b = 2, c = 3, d = 4, e = 5, f = 6))
b < array(1:prod(c(1, 2, 3, 5, 6)), dim = c(a = 1, b = 2, c = 3, e = 5, f = 6))
 x < ""
 y < ""

f < function(a, b) {
 x << names(dim(a))
 y << names(dim(b))
+ stopifnot(identical(names(dim(a)), c('b', 'c', 'd')))
+ stopifnot(identical(names(dim(b)), c('b', 'c')))
}
r < multiApply::Apply(list(a, b),
list(c('b', 'c', 'd'),
c('b', 'c')),
f)

 expect_equal(x, c('b', 'c', 'd'))
 expect_equal(y, c('b', 'c'))
})
# TODOS: