Fix PugsFunctionAdapter for R^p -> R^dxd functions

src/language/utils/PugsFunctionAdapter.hpp

@@ -236,8 +236,8 @@ class PugsFunctionAdapter<OutputType(InputType...)>
           AggregateDataVariant& v = std::get<AggregateDataVariant>(result);
           OutputType x;
 
-          for (size_t i = 0, l = 0; i < x.dimension(); ++i) {
-            for (size_t j = 0; j < x.dimension(); ++j, ++l) {
+          for (size_t i = 0, l = 0; i < x.nbRows(); ++i) {
+            for (size_t j = 0; j < x.nbColumns(); ++j, ++l) {
               std::visit(
                 [&](auto&& Aij) {
                   using Aij_T = std::decay_t<decltype(Aij)>;