---
title: "Parallel Processing"
output:
    rmarkdown::html_vignette:
        toc: true
description: >
  Leverage parallel processing to speed up the metasnf pipeline.
vignette: >
  %\VignetteIndexEntry{Parallel Processing}
  %\VignetteEngine{knitr::rmarkdown}
  %\VignetteEncoding{UTF-8}
---

<style>
div.aside { background-color:#fff2e6; }
</style>

```{r, include = FALSE}
# Default chunk options
knitr::opts_chunk$set(
    collapse = TRUE,
    comment = "#>",
    fig.width = 6,
    fig.height = 4.5,
    fig.align = "center"
)
```

Download a copy of the vignette to follow along here: [parallel_processing.Rmd](https://raw.githubusercontent.com/BRANCHlab/metasnf/main/vignettes/parallel_processing.Rmd)

metasnf can use [future apply](https://future.apply.futureverse.org/) to run the rows of a settings matrix in parallel.

Parallel processing has some overhead cost, so it's possible that your particular situation won't greatly benefit (or may even run slower) with parallel processing enabled.

However, if you have access to multiple cores and each integration you are running is quite intensive, the parallel processing will likely help speed things up quite a lot.

## Basic usage

```{r eval = FALSE}
# Load the package
library(metasnf)

# Setting up the data
data_list <- generate_data_list(
    list(cort_t, "cortical_thickness", "neuroimaging", "continuous"),
    list(cort_sa, "cortical_surface_area", "neuroimaging", "continuous"),
    list(subc_v, "subcortical_volume", "neuroimaging", "continuous"),
    list(income, "household_income", "demographics", "continuous"),
    list(pubertal, "pubertal_status", "demographics", "continuous"),
    uid = "unique_id"
)

# Specifying 5 different sets of settings for SNF
set.seed(42)
settings_matrix <- generate_settings_matrix(
    data_list,
    nrow = 10,
    max_k = 40
)

solutions_matrix <- batch_snf(
    data_list,
    settings_matrix,
    processes = "max" # Can also be a specific integer
)
```

## Including a progress bar

Use the progressr package to visualize the progress of parallel batch_snf.

```{r eval = FALSE}
progressr::with_progress({
    solutions_matrix <- batch_snf(
        data_list,
        settings_matrix,
        processes = "max"
    )
})
```

## Number of processes

Setting processes to "max" will make use of as many cores as R can find.

If you want to dial things back a little, you can specify more precisely the number of processes you want:

```{r eval = FALSE}
solutions_matrix <- batch_snf(
    data_list,
    settings_matrix,
    processes = 4
)
```

To find out how many processes you have access to (or at least, how many metasnf will think you have access to), use the `availableCores()` function from the future package:

```{r eval = FALSE}
library(future)

availableCores()
```