-Copyright (c) 2016-2017, Mor-Absa Loum ; 2017, Benjamin Auder
-
-Permission is hereby granted, free of charge, to any person obtaining
-a copy of this software and associated documentation files (the
-"Software"), to deal in the Software without restriction, including
-without limitation the rights to use, copy, modify, merge, publish,
-distribute, sublicense, and/or sell copies of the Software, and to
-permit persons to whom the Software is furnished to do so, subject to
-the following conditions:
-
-The above copyright notice and this permission notice shall be
-included in all copies or substantial portions of the Software.
-
-THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
-EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
-MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
-NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
-LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
-OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
-WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+YEAR: 2016-2018
+COPYRIGHT HOLDER: Mor-Absa Loum, Benjamin Auder
--- /dev/null
+// RegisteringDynamic Symbols
+
+#include <R.h>
+#include <Rinternals.h>
+#include <R_ext/Rdynload.h>
+
+void R_init_morpheus(DllInfo* info) {
+ R_registerRoutines(info, NULL, NULL, NULL, NULL);
+ R_useDynamicSymbols(info, TRUE);
+}
test_that("naive computation provides the same result as vectorized computations",
{
h <- 1e-7 #for finite-difference tests
- tol <- .25 * sqrt(h) #about 7.9 e-5
+ tol <- 1e-3 #large tolerance, necessary in some cases... (generally 1e-6 is OK)
for (dK in list( c(2,2), c(5,3)))
{
d = dK[1]