Francesc Alted FrancescAlted

CEO of https://ironarray.io Project leader of the Blosc project. Compression and efficient data I/O to efficient storage and computation.

FrancescAlted / circle_func.py

Created September 12, 2022 07:18

	@udf.jit()
	def filter_func(out: udf.Array(udf.float32, 2),
	vals: udf.Array(udf.float32, 2),
	nrows: udf.int64, ncols: udf.int64,
	iscircle: udf.bool) -> udf.int32:
	n = out.window_shape[0]
	m = out.window_shape[1]
	row_start = out.window_start[0]
	col_start = out.window_start[1]
	for i in range(n):

FrancescAlted / create_square_board.py

Created September 12, 2022 07:23

	expr = ia.expr_from_udf(filter_func,
	[rand_data],
	[shape[0], shape[1], False])
	square = expr.eval()

FrancescAlted / create_circle_board.py

Created September 12, 2022 07:25

	expr = ia.expr_from_udf(filter_func,
	[rand_data],
	[shape[0], shape[1], True])
	circle = expr.eval()

FrancescAlted / print_pi_estimate.py

Created September 12, 2022 07:27

	area_circle = ia.nansum(circle)
	area_square = ia.nansum(square)
	print(f"PI estimate: {4 * area_circle / area_square}")

FrancescAlted / blosc2-jit-microbench.ipynb

Created March 10, 2025 18:08

Micro-benchmark for blosc2, using numpy functions, reductions and broadcasting

Sorry, something went wrong. Reload?

Sorry, we cannot display this file.

Sorry, this file is invalid so it cannot be displayed.

FrancescAlted / big-calculation2.py

Last active September 23, 2025 05:11

A calculation of a 1000x1000x1000 cube using NumPy and Blosc2

	import os
	import psutil
	from time import time
	import numpy as np
	import blosc2


	# --- Experiment Setup ---
	n_frames = 1000 # Raise this for more frames
	width, height = np.array((n_frames, n_frames)) # Size of the grid