wball/Circle-Packings
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

asdf

Browse source
This commit is contained in:
William Ball 2021-06-27 12:05:41 -04:00
parent e08e34a7c4
commit 9e00cffb28
7 changed files with 253843 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
fractal_dimension/bai_finch_rust/Cargo.toml
View file

@ -11,3 +11,6 @@ num = "0.3.1"
nalgebra = "0.27.1" nalgebra = "0.27.1"
num-complex = "0.3.1" num-complex = "0.3.1"
rayon = "1.5.1" rayon = "1.5.1"
[profile.release]
debug = true

BIN
fractal_dimension/bai_finch_rust/flamegraph.png Normal file
View file

Binary file not shown.
Side by side

After

Width:  |  Height:  |  Size: 180 KiB

419
fractal_dimension/bai_finch_rust/flamegraph.svg Normal file
View file

File diff suppressed because one or more lines are too long
Side by side

After

Width:  |  Height:  |  Size: 76 KiB

BIN
fractal_dimension/bai_finch_rust/perf.data Normal file
View file

Binary file not shown.

BIN
fractal_dimension/bai_finch_rust/perf.data.old Normal file
View file

Binary file not shown.

253413
fractal_dimension/bai_finch_rust/perf.data.txt Normal file
View file

File diff suppressed because it is too large Load diff

15
fractal_dimension/bai_finch_rust/src/main.rs
View file

@ -1,12 +1,13 @@
#![allow(dead_code)] #![allow(dead_code)]
use num_complex::Complex64; use num_complex::Complex64;
use nalgebra::{VecStorage, Const};
use rayon::prelude::*; use rayon::prelude::*;
type F = f64; type F = f64;
type Matrix2x2 = nalgebra::SMatrix<Complex64, 2, 2>; type Matrix2x2 = nalgebra::SMatrix<Complex64, 2, 2>;
type MatrixBig = nalgebra::SMatrix<Complex64, NC2, NC2>; type MatrixBig = nalgebra::Matrix<Complex64, Const<NC2>, Const<NC2>, VecStorage<Complex64, Const<NC2>, Const<NC2>>>;
type MatrixBigr = nalgebra::SMatrix<f64, NC2, NC2>; type MatrixBigr = nalgebra::Matrix<f64, Const<NC2>, Const<NC2>, VecStorage<f64, Const<NC2>, Const<NC2>>>;
const G: Matrix2x2 = Matrix2x2::new( const G: Matrix2x2 = Matrix2x2::new(
Complex64::new(2.0 / 12.0, -2.0 / 12.0), Complex64::new(2.0 / 12.0, -2.0 / 12.0),
@ -22,9 +23,9 @@ const R: Matrix2x2 = Matrix2x2::new(
Complex64::new(-6.0 / 12.0, -8.0 / 12.0), Complex64::new(-6.0 / 12.0, -8.0 / 12.0),
); );
const NC: usize = 15; const NC: usize = 10;
const K0: i32 = 100; const K0: i32 = 100;
const LC: usize = 7; const LC: usize = 5;
const NC2: usize = NC * NC; const NC2: usize = NC * NC;
fn fancy_l(q: f64) -> MatrixBigr { fn fancy_l(q: f64) -> MatrixBigr {
@ -58,7 +59,7 @@ fn normal_f(q: f64, n: i32, s: i32, l: i32) -> Complex64 {
let g21 = G[(1, 0)]; let g21 = G[(1, 0)];
let g22 = G[(1, 1)]; let g22 = G[(1, 1)];
(0..=s) (0..=s)
.into_par_iter() .into_iter()
.map(|j: i32| -> Complex64 { .map(|j: i32| -> Complex64 {
non_integral_choose(-n as f64 - q, j) non_integral_choose(-n as f64 - q, j)
* integral_choose(n, s - j) * integral_choose(n, s - j)
@ -101,7 +102,7 @@ fn normal_f(q: f64, n: i32, s: i32, l: i32) -> Complex64 {
fn zeta(s: f64, k0: i32) -> f64 { fn zeta(s: f64, k0: i32) -> f64 {
const UPPER_BOUND: i32 = 40_000; const UPPER_BOUND: i32 = 40_000;
(k0..=UPPER_BOUND) (k0..=UPPER_BOUND)
.into_par_iter() .into_iter()
.map(|j| (j as f64).powf(-s)) .map(|j| (j as f64).powf(-s))
.sum::<f64>() .sum::<f64>()
- (UPPER_BOUND as f64).powf(1.0 - s) / (1.0 - s) - (UPPER_BOUND as f64).powf(1.0 - s) / (1.0 - s)
@ -144,7 +145,7 @@ fn normal_m(k: i32, q: f64, n: i32, s: i32) -> Complex64 {
let a22 = ak[(1, 1)]; let a22 = ak[(1, 1)];
(0..=s) (0..=s)
.into_par_iter() .into_iter()
.map(|j| { .map(|j| {
non_integral_choose(-n as f64 - q, j) non_integral_choose(-n as f64 - q, j)
* integral_choose(n, s - j) * integral_choose(n, s - j)