bschneidr

library(rextendr)

rust_source(
  code = readLines("r-ndarray-eigendecomposition.rs"),
  dependencies = list("ndarray-linalg" = "0.17.0"),
  features = "ndarray",
  profile = "dev"
)

X_dim <- 5

bschneidr

library(rextendr)

rust_source(
  code         = readLines("r-faer-eigendecomposition.rs"),
  dependencies = list("faer" = "0.22.6"),
  features     = "faer",
  profile      = "release"
)

X <- matrix(rnorm(n = 25), 5, 5)

bschneidr

library(rextendr) # Version 0.4.0

rust_code <- r"(
use faer::prelude::*;
use extendr_api::prelude::*;

#[extendr]
fn return_matrix(X: RMatrix<f64>) -> RMatrix<f64> {
  let A: Mat<f64> = Mat::from_fn(X.nrows(), X.ncols(), |row, col| X[[row, col]]);
  let result: RMatrix<f64> = RMatrix::new_matrix(A.nrows(), A.ncols(), |row, col| A[(row, col)]);

bschneidr

import polars as pl
import pandas as pd

# Read input data as a data frame with two columns
input = pl.from_pandas(
  pd.read_table("day-1/input", sep = "   ", header = None, 
                names = ["list_1", "list_2"])
)

# Sort each column separately, 

bschneidr

read.table("day-1/input", col.names = c("list_1", "list_2")) |>
  apply(MARGIN = 2, sort) |>
  apply(MARGIN = 1, FUN = \(x) abs(x[2] - x[1])) |>
  sum()

bschneidr

sudo apt-get install libopenblas-dev
sudo apt-get install liblapacke-dev
sudo apt-get install liblapacke
sudo apt-get install libopenblas-serial-dev
sudo apt-get install libopenblas0-serial
sudo apt-get install libopenblas0

bschneidr

library(dplyr)
library(tidyr)
library(survey)

# Load example data
# 'apipop' has one row per school

  data('api', package = 'survey')

# Draw a sample of school districts

bschneidr

library(dplyr)
library(broom)
library(tidyr)

# Make example data
  my_data <- data.frame(x = sample(1:5, size = 10, replace = TRUE),
                        y = sample(1:5, size = 10, replace = TRUE))

# Generate every pair of variables
  var_list <- colnames(my_data)

bschneidr

describe_design(
  srs_stage(id = "PSU_ID", method = "SRS",
            stratum = "PSU_STRATUM_ID"),
  nonresponse_stage(response_indicator = "PSU_RESPONDENT"),
  pps_stage(id = "SSU_ID", method = "PPS",
            stratum = "SSU_STRATUM_ID"),
  nonresponse_stage(response_indicator = "SSU_RESPONDENT")
)

bschneidr

n = 5

# Establish quadratic form

quad_form_matrix <- matrix(nrow = n, ncol = n)
for (i in seq(n)) {
  for (j in seq(n)) {
    if (i == j) {
      quad_form_matrix[i,j] <- (1/n)
    } else {