adic 0.5.1

Arithmetic and rootfinding for p-adic numbers
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63

//! Create and manipulate multi-prime adics, e.g. 10-adics

use clap::Parser;
use itertools::Itertools;
use adic::{
    divisible::Composite,
    error::AdicResult,
    num_adic::MAdic,
    traits::HasDigits,
    ZAdic,
};


#[derive(Parser)]
struct Args {

    #[arg(short)]
    mp: u32,

    #[arg(short, allow_hyphen_values=true)]
    a: i32,

    #[arg(long)]
    precision: usize,

}


fn main() -> AdicResult<()> {

    let args = Args::parse();
    let ma = MAdic::approx_from_i32(args.mp, args.a, args.precision)?;
    println!("{ma}");

    let digit_vec = ma.digits().collect::<Vec<_>>();
    println!("Digits:");
    println!("{}", digit_vec.into_iter().rev().map(|d|
        if d < 10 {
            d.to_string()
        } else {
            format!("[{d}]")
        }
    ).collect::<String>());

    let c = Composite::from(args.mp);

    println!("Prime idempotents:");
    let primes = c.primes();
    let prime_idemps = MAdic::<ZAdic>::prime_idempotents(args.mp, args.precision)?;
    for (prime, idemp) in primes.zip(prime_idemps) {
        println!("Prime {prime}, idemp {}", idemp.to_str_radix(args.mp));
    }

    println!("All idempotents:");
    let primes = c.primes();
    let all_idemps = MAdic::<ZAdic>::all_idempotents(args.mp, args.precision)?;
    for (primes, idemp) in primes.powerset().zip(all_idemps) {
        println!("Primes {}, idemp {}", primes.into_iter().join(" x "), idemp.to_str_radix(args.mp));
    }

    Ok(())

}