modular_forms 0.0.1 → 0.0.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: c2b8536ed2846f289afd967456cdd892858ad3af312acffab3e312c57b367a38
-  data.tar.gz: ea900a9133578ce1278141bda388d85fc222419ac2de207c50a10c243176166f
+  metadata.gz: f9b49161da093d436c9afd089b0d2935a43fafb71d1b1eece1af75cc3ef08833
+  data.tar.gz: 54672d7ef27bc439d12747a919a024e9c9fa58d76a3d4b3ae1fa3f807226be6b
 SHA512:
-  metadata.gz: 7542e5fa3105f3008bdf192c2471e950898567aca9c133ebb50c7e0d65ec1bd0c3b1765c944257c59782bd26f2f849f15d7653652d14844d844a91ab1174b7b9
-  data.tar.gz: 6878c24ff277884ff52aeb3ad54073ae50001b3f0f792b317b1d4ed7caf40151a867124a394eef0731846f1d3732b096f741e4e8c6e2473f9b172eddaed0d237
+  metadata.gz: d8ad361a6edc83f140ac815bb4b04a6a627100204f523e2192edfaa4cdc8e7b1022cd1c4e6651526582d287818231193d546dd2fb99b859d39401cb5002b3861
+  data.tar.gz: 2a7223689108395ebc2f389d1e1ae9b6e61a5165e9da9885624b7758500131ca01f983c433dc045928c433608d773b6ea22410d2196b43b2261cb88c40933b27

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 Edgar Armando Delgado Vega
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/README.md

@@ -2,44 +2,56 @@
 
 A creative toolkit for exploring modular forms and elliptic curves through [Sonic Pi](https://sonic-pi.net/).
 
-![Modular Forms - Image](modular_forms.png)
+<p align="center">
+  <img src="https://img.shields.io/gem/v/modular_forms" alt="Gem Version">
+  <img src="https://img.shields.io/gem/dt/modular_forms" alt="Gem Total Downloads">
+  <img src="https://img.shields.io/github/stars/edelveart/modular_forms" alt="GitHub Repo stars">
+  <img src="https://img.shields.io/github/license/edelveart/modular_forms" alt="GitHub License">
+</p>
 
-## Project Status
+[![Modular Forms - Image](https://raw.githubusercontent.com/edelveart/modular_forms/main/modular_forms.png)](https://rubygems.org/gems/modular_forms)
+
+## ⚙️ Project Status
 
 This is a pre-alpha release of `modular_forms`. At this stage, only a subset of core mathematical definitions and operations is implemented.
-Future updates might include a DSL, depending on how the library is used and the interest from the community.
 
-Contributions are welcome! Feel free to fork, make changes, and submit a pull request.
+Future updates might include a DSL, depending on how the library is used and the interest from the community. A key challenge lies in creating musical mappings that stay true to the underlying mathematics while also sounding intentional, expressive, and naturally fitting within a musical structure.
 
-## Features
+## 🧊 Features
 
 - **Accessible to both musicians and coders**: No math expertise required. Create musical patterns, rhythms, timbres, and harmonies by experimenting with mathematical ideas and turning them into sound and effects intuitively.
 - **Interactive Educational Resource**: Use **Sonic Pi** to discover introductory number theory in a hands-on, immersive way, gaining insights into abstract concepts through math in action.
 
-## Purpose and Scope
+## 💡 Purpose and Scope
+
+### Overview
 
-Given the vastness of the field, this tool intentionally focuses on a limited subset of definitions, without covering all aspects of each. Below is a list of the implemented modules:
+Given the vastness of the field, this tool intentionally focuses on a limited subset of definitions, without covering all aspects of each.
 
-- [Eisenstein Series](#eisenstein-series)
-- [Eta Functions and Eta Quotients](#eta-functions-and-eta-quotients)
-- [Theta functions](#theta-functions)
-- [Ramanujan Tau Function](#ramanujan-tau-function)
-- [J-Function](#j-function)
-- [Hecke Operators](#hecke-operators)
-- [SL(2,Z) Group](#sl2z-group)
-- [Dirichlet Characters](#dirichlet-characters)
-- [Elliptic Curves over Rationals](#elliptic-curves-over-rationals)
-- [Elliptic Curves over Finite Fields](#elliptic-curves-over-finite-fields)
-- [Newforms Invariants](#newforms-invariants)
+- 🧩 List of implemented modules:
+  - [Eisenstein Series](#eisenstein-series)
+  - [Eta Functions and Eta Quotients](#eta-functions-and-eta-quotients)
+  - [Theta Functions](#theta-functions)
+  - [Ramanujan Tau Function](#ramanujan-tau-function)
+  - [J-Function](#j-function)
+  - [Hecke Operators](#hecke-operators)
+  - [SL(2,Z) Group](#sl2z-group)
+  - [Dirichlet Characters](#dirichlet-characters)
+  - [Elliptic Curves over Rationals](#elliptic-curves-over-rationals)
+  - [Elliptic Curves over Finite Fields](#elliptic-curves-over-finite-fields)
+  - [Newform Invariants](#newform-invariants)
+  - [L-functions](#l-functions)
+  - [p-adic Fields](#p-adic-fields)
 
 ### Not Optimized for Computational Efficiency
 
-This library is designed for creative exploration rather than maximum computational efficiency. It is not intended to replace advanced mathematical software. Instead, it draws inspiration from tools like SageMath, Pari/GP, and the LMFDB database.
+This library is designed for creative exploration rather than maximum computational efficiency. It is **not intended** to replace specialized mathematical software optimized for heavy or large-scale computations. Instead, it draws inspiration from tools like **SageMath**, **Pari/GP**, and the **LMFDB database**.
 
 ### Goal
+
 The goal is simple: to provide an accessible and creative starting point for those who wish to explore, learn, and uncover new ideas, regardless of their mathematical background.
 
-## Installation
+## 💎 Installation
 
 You can install the `modular_forms` gem directly from **RubyGems** or clone it from GitHub.
 
@@ -47,31 +59,100 @@ You can install the `modular_forms` gem directly from **RubyGems** or clone it f
 gem install modular_forms
 ```
 
-## How to use?
+If you are **use Ruby**, then import via
+```rb
+require 'modular_forms'
+```
 
-You can dive into the beauty of math, both in Ruby and Sonic Pi, creating music in real-time. Here is a simple example of how to use **modular_forms** to generate a basic musical pattern:
+If you are **using Sonic Pi**, replace `<PATH>` with the directory path to the `modular_forms.rb` file inside the gem installation on your system. You can find the full path by running:
+
+```bash
+gem which modular_forms
+```
+
+Then require the file like this:
 
 ```rb
-# If you are using Ruby, simply load the gem with its name
-# after installing it via 'gem install modular_forms'.
-require 'modular_forms'
+require "<PATH>/modular_forms.rb"
+```
+
+## 🧰 How to use?
+
+### 🎶 Explore Fermat’s Last Theorem in Sonic Pi
+
+The Taniyama–Shimura Conjecture (now the Modularity Theorem), proven for semistable cases by Andrew Wiles, connects elliptic curves over the rationals to modular forms. This profound result was crucial in proving **Fermat’s Last Theorem**.
+
+For example, we can explore the relationship between the elliptic curve `144.a3` and the newform orbit `144.2.a.a` through their L-function,
+`L(E, s) = L(f, s)`, in a musical sense.
 
-# If you are using Sonic Pi, replace <PATH>
-# with the full path to the 'modular_forms.rb' file
-# inside the gem installation on your system.
+```rb
 require "<PATH>/modular_forms.rb"
 
-# Calculate the Eisenstein series of weight k = 4
-eisenstein_melody = ModularForms.eisenstein_serie(4)
+# Set precision for the Fourier q-expansion
+prec = 20
+
+# Construct a weight 2 newform as an eta quotient
+n = ModularForms.dedekind_eta_pow(12, prec, 12)
+d1 = ModularForms.dedekind_eta_pow(4, prec, 6)
+d2 = ModularForms.dedekind_eta_pow(4, prec, 24)
+eta_prod = ModularForms.eta_product(d1, d2)
+
+newform = ModularForms.eta_quotient(n, eta_prod, prec)
+
+# Define the elliptic curve E over F_p
+p = 13
+ellc = ModularForms.elliptic_curve_fp(p, [0, -1]) # y^2 = x^3 - 1
+points = ModularForms.cardinality_fp(ellc) # count points on E mod p
+
+# Compute the modular coefficient a_p of L-function(E, s)
+a_p = ModularForms.a_p(p, points)
 
-# Play the melody in a loop with a mathematical transformation
-120.times do
-  play eisenstein_melody.next % 12 * 7
-  sleep 0.5
+# Sonify the modular relationship
+live_loop :modularity_music do
+  play (chord (:a3 + a_p), :m11)[newform.ring.tick],
+    amp: a_p, release: 0.125
+  sleep 0.125
 end
 ```
 
-## Implemented Modular Forms, Elliptic Curves, and Related Definitions
+### 🎶 Merging Machine and Organism - Warp Up
+
+Unlike the structured Fermat example, this finite loop explores more abstract territory, using multiple concepts to create a sonic landscape that, while grounded in mathematical principles, goes beyond conventional musical forms.
+
+```rb
+require "<PATH>/modular_forms.rb"
+
+eisenstein_melody = ModularForms.eisenstein_series(8)
+p = 3
+ellc = ModularForms.elliptic_curve_q([2, 5])
+disc = ModularForms.discriminant_q(ellc)
+e8 = eisenstein_melody.take(58)
+hecke_op = ModularForms.hecke_operator_prime_non_cusp(e8, p, 8, 20)
+j_func = ModularForms.j_function(40)
+newform_ac = ModularForms.analytic_conductor(15, 2)
+pol = ModularForms.def_pol_2deg(41)
+
+(disc * -1).times do
+  synth :subpulse, note: ModularForms.padic_valuation(eisenstein_melody.next, p) % 7 + 70,
+    release: 0.25 if (spread (disc % 6), 7).tick(:d)
+
+  synth :chiplead, note: hecke_op.tick(:tp) % 7 + 50,
+    release: newform_ac, attack: newform_ac
+
+  synth :fm, note: ModularForms.gauss_sum_triv(5, j_func.tick(:j)) + 51,
+    release: 0.25, attack: 0.012,
+    pan: ModularForms.analytic_conductor(j_func.look % 15, 4) * rrand_i(-1, 1) if
+  spread(ModularForms.gamma1_index(p), 11).tick(:g)
+
+  sample 90, release: 0.125 if pol.tick == 'x'
+  sample :ambi_drone, beat_stretch: 0.7 if pol.look == '*'
+  sample :bass_drop_c, release: 0.125 if pol.look == '6'
+
+  sleep ModularForms.padic_norm(p, p)
+end
+```
+
+## 🧩 Implemented Modular Forms, Elliptic Curves, and Related Definitions
 
 ### Eisenstein Series
 
@@ -87,7 +168,7 @@ end
 7. `ModularForms.eta_product(eta1, eta2, prec = nil)`
 8.  `ModularForms.eta_quotient(num_eta, den_eta, prec)`
 
-### Theta functions
+### Theta Functions
 
 9. `ModularForms.jacobi_theta_function(jacobi_index = 3, square_coefs = false)`
 10.  `ModularForms.jacobi_theta_function_pow(jacobi_index, power, prec)`
@@ -147,6 +228,32 @@ end
 43. `ModularForms.cardinality_fp(curve)`
 44. `ModularForms.quadratic_twist_fp(curve)`
 
-### Newforms Invariants
+### Newform Invariants
+
 45. `ModularForms.analytic_conductor(level_n, weight_k)`
 
+### L-functions
+
+46. `ModularForms.a_p(p, cardinality)`
+
+### p-adic Fields
+
+47. `ModularForms.padic_valuation(num_b10, p)`
+48. `ModularForms.padic_norm(num_b10, p)`
+49. `ModularForms.padic_expansion(num_b10, p, prec = 11, reverse_trim = false)`
+50. `ModularForms.def_pol_2deg(p = 2, c = 0, num = 1)`
+
+## 🧪 Testing
+
+Install dependencies first:
+
+```bash
+bundle install
+```
+
+This project features a comprehensive Minitest suite covering core functionality across modules. While not all edge cases are tested, the main mathematical functions are well validated to ensure correctness and stability.
+
+Run tests with:
+```bash
+rake test
+```

data/lib/modular_forms/core/dirichlet_characters.rb

@@ -1,7 +1,6 @@
 # frozen_string_literal: true
 
 require_relative '../core/numeric_helpers/numeric_helpers'
-require 'bigdecimal'
 
 module ModularForms
   module Core

data/lib/modular_forms/core/l_functions.rb

@@ -0,0 +1,15 @@
+# frozen_string_literal: true
+
+module ModularForms
+  module Core
+    # ModularForms::Core::LFunctions
+    #
+    # The LFunctions module provides methods related to L-functions,
+    # particularly those associated with modular forms and elliptic curves.
+    module LFunctions
+      def self.a_p(p, cardinality)
+        p + 1 - cardinality
+      end
+    end
+  end
+end

data/lib/modular_forms/core/padic_fields.rb

@@ -0,0 +1,68 @@
+# frozen_string_literal: true
+
+require 'prime'
+require_relative './padic_pol/padic_pol'
+
+module ModularForms
+  module Core
+    # ModularForms::Core::PAdicFields
+    #
+    # This module provides concepts and functionality related to p-adic fields.
+    module PAdicFields
+      def self.padic_valuation(num_b10, p)
+        raise "#{p} is not a prime number" unless Prime.prime?(p)
+
+        index = 0
+        return index if num_b10 < 1
+
+        index += 1 while num_b10 % (p**index) == 0 # rubocop:disable Style/NumericPredicate
+        index - 1
+      end
+
+      def self.padic_norm(num_b10, p)
+        return 0 if num_b10 == 0 # rubocop:disable Style/NumericPredicate
+
+        power = padic_valuation(num_b10, p)
+        1.0 / p**power
+      end
+
+      def self.padic_expansion(num_b10, p, prec = 11, reverse_trim = false) # rubocop:disable Metrics/AbcSize,Metrics/MethodLength,Style/OptionalBooleanParameter
+        inverse_limit_arr = []
+        (1...prec).each do |index|
+          power = p**index
+          inverse_limits = ((num_b10 % power) + power) % power
+          inverse_limit_arr << inverse_limits
+        end
+        p_adic_expansion_array = [inverse_limit_arr[0]]
+
+        (0...(prec - 2)).each do |index|
+          a1 = inverse_limit_arr[index + 1]
+          a0 = inverse_limit_arr[index]
+          power = p**(index + 1)
+          digits = (a1 - a0) / power
+          p_adic_expansion_array << digits
+        end
+
+        if reverse_trim
+          reversed = p_adic_expansion_array.reverse
+          first_nonzero_index = reversed.index { |digit| digit != 0 } || 0
+          reversed[first_nonzero_index..-1] # rubocop:disable Style/SlicingWithRange
+        else
+          p_adic_expansion_array
+        end
+      end
+
+      def self.def_pol_2deg(p = 2, c = 0, num = 1)
+        raise ArgumentError, "Residue field characteristic #{p} must be a prime number" unless Prime.prime?(p)
+
+        unless (0..3).include?(c)
+          raise ArgumentError,
+                "The specified discriminant exponent (#{c}) is invalid. Must be in range (0..3)."
+        end
+        raise ArgumentError, "Invalid value: #{num}. Must be in range (1..4)." unless (1..4).include?(num)
+
+        PAdicPol.get_defining_polynomial(PAdicPol::DEF_POL_HASH_2DEG, p, c, num)
+      end
+    end
+  end
+end

data/lib/modular_forms/core/padic_pol/data_pol_2deg.rb

@@ -0,0 +1,152 @@
+# frozen_string_literal: true
+
+module ModularForms
+  module Core
+    module PAdicFields
+      module PAdicPol # rubocop:disable Metrics/ModuleLength
+        DATA_POL_2DEG = { '2.2.0.1' => [1, 1, 1],
+                          '2.2.2.1' => [2, 2, 1],
+                          '2.2.2.2' => [6, 2, 1],
+                          '2.2.3.1' => [2, 4, 1],
+                          '2.2.3.2' => [10, 4, 1],
+                          '2.2.3.3' => [2, 0, 1],
+                          '2.2.3.4' => [10, 0, 1],
+                          '3.2.0.1' => [2, 2, 1],
+                          '3.2.1.1' => [6, 0, 1],
+                          '3.2.1.2' => [3, 0, 1],
+                          '5.2.0.1' => [2, 4, 1],
+                          '5.2.1.1' => [5, 0, 1],
+                          '5.2.1.2' => [10, 0, 1],
+                          '7.2.0.1' => [3, 6, 1],
+                          '7.2.1.1' => [21, 0, 1],
+                          '7.2.1.2' => [7, 0, 1],
+                          '11.2.0.1' => [2, 7, 1],
+                          '11.2.1.1' => [22, 0, 1],
+                          '11.2.1.2' => [11, 0, 1],
+                          '13.2.0.1' => [2, 12, 1],
+                          '13.2.1.1' => [13, 0, 1],
+                          '13.2.1.2' => [26, 0, 1],
+                          '17.2.0.1' => [3, 16, 1],
+                          '17.2.1.1' => [17, 0, 1],
+                          '17.2.1.2' => [51, 0, 1],
+                          '19.2.0.1' => [2, 18, 1],
+                          '19.2.1.1' => [38, 0, 1],
+                          '19.2.1.2' => [19, 0, 1],
+                          '23.2.0.1' => [5, 21, 1],
+                          '23.2.1.1' => [115, 0, 1],
+                          '23.2.1.2' => [23, 0, 1],
+                          '29.2.0.1' => [2, 24, 1],
+                          '29.2.1.1' => [29, 0, 1],
+                          '29.2.1.2' => [58, 0, 1],
+                          '31.2.0.1' => [3, 29, 1],
+                          '31.2.1.1' => [93, 0, 1],
+                          '31.2.1.2' => [31, 0, 1],
+                          '37.2.0.1' => [2, 33, 1],
+                          '37.2.1.1' => [37, 0, 1],
+                          '37.2.1.2' => [74, 0, 1],
+                          '41.2.0.1' => [6, 38, 1],
+                          '41.2.1.1' => [41, 0, 1],
+                          '41.2.1.2' => [123, 0, 1],
+                          '43.2.0.1' => [3, 42, 1],
+                          '43.2.1.1' => [86, 0, 1],
+                          '43.2.1.2' => [43, 0, 1],
+                          '47.2.0.1' => [5, 45, 1],
+                          '47.2.1.1' => [235, 0, 1],
+                          '47.2.1.2' => [47, 0, 1],
+                          '53.2.0.1' => [2, 49, 1],
+                          '53.2.1.1' => [53, 0, 1],
+                          '53.2.1.2' => [106, 0, 1],
+                          '59.2.0.1' => [2, 58, 1],
+                          '59.2.1.1' => [118, 0, 1],
+                          '59.2.1.2' => [59, 0, 1],
+                          '61.2.0.1' => [2, 60, 1],
+                          '61.2.1.1' => [61, 0, 1],
+                          '61.2.1.2' => [122, 0, 1],
+                          '67.2.0.1' => [2, 63, 1],
+                          '67.2.1.1' => [134, 0, 1],
+                          '67.2.1.2' => [67, 0, 1],
+                          '71.2.0.1' => [7, 69, 1],
+                          '71.2.1.1' => [497, 0, 1],
+                          '71.2.1.2' => [71, 0, 1],
+                          '73.2.0.1' => [5, 70, 1],
+                          '73.2.1.1' => [73, 0, 1],
+                          '73.2.1.2' => [365, 0, 1],
+                          '79.2.0.1' => [3, 78, 1],
+                          '79.2.1.1' => [237, 0, 1],
+                          '79.2.1.2' => [79, 0, 1],
+                          '83.2.0.1' => [2, 82, 1],
+                          '83.2.1.1' => [166, 0, 1],
+                          '83.2.1.2' => [83, 0, 1],
+                          '89.2.0.1' => [3, 82, 1],
+                          '89.2.1.1' => [89, 0, 1],
+                          '89.2.1.2' => [267, 0, 1],
+                          '97.2.0.1' => [5, 96, 1],
+                          '97.2.1.1' => [97, 0, 1],
+                          '97.2.1.2' => [485, 0, 1],
+                          '101.2.0.1' => [2, 97, 1],
+                          '101.2.1.1' => [101, 0, 1],
+                          '101.2.1.2' => [202, 0, 1],
+                          '103.2.0.1' => [5, 102, 1],
+                          '103.2.1.1' => [309, 0, 1],
+                          '103.2.1.2' => [103, 0, 1],
+                          '107.2.0.1' => [2, 103, 1],
+                          '107.2.1.1' => [214, 0, 1],
+                          '107.2.1.2' => [107, 0, 1],
+                          '109.2.0.1' => [6, 108, 1],
+                          '109.2.1.1' => [109, 0, 1],
+                          '109.2.1.2' => [218, 0, 1],
+                          '113.2.0.1' => [3, 101, 1],
+                          '113.2.1.1' => [113, 0, 1],
+                          '113.2.1.2' => [339, 0, 1],
+                          '127.2.0.1' => [3, 126, 1],
+                          '127.2.1.1' => [381, 0, 1],
+                          '127.2.1.2' => [127, 0, 1],
+                          '131.2.0.1' => [2, 127, 1],
+                          '131.2.1.1' => [262, 0, 1],
+                          '131.2.1.2' => [131, 0, 1],
+                          '137.2.0.1' => [3, 131, 1],
+                          '137.2.1.1' => [137, 0, 1],
+                          '137.2.1.2' => [411, 0, 1],
+                          '139.2.0.1' => [2, 138, 1],
+                          '139.2.1.1' => [278, 0, 1],
+                          '139.2.1.2' => [139, 0, 1],
+                          '149.2.0.1' => [2, 145, 1],
+                          '149.2.1.1' => [149, 0, 1],
+                          '149.2.1.2' => [298, 0, 1],
+                          '151.2.0.1' => [6, 149, 1],
+                          '151.2.1.1' => [453, 0, 1],
+                          '151.2.1.2' => [151, 0, 1],
+                          '157.2.0.1' => [5, 152, 1],
+                          '157.2.1.1' => [157, 0, 1],
+                          '157.2.1.2' => [314, 0, 1],
+                          '163.2.0.1' => [2, 159, 1],
+                          '163.2.1.1' => [326, 0, 1],
+                          '163.2.1.2' => [163, 0, 1],
+                          '167.2.0.1' => [5, 166, 1],
+                          '167.2.1.1' => [835, 0, 1],
+                          '167.2.1.2' => [167, 0, 1],
+                          '173.2.0.1' => [2, 169, 1],
+                          '173.2.1.1' => [173, 0, 1],
+                          '173.2.1.2' => [346, 0, 1],
+                          '179.2.0.1' => [2, 172, 1],
+                          '179.2.1.1' => [358, 0, 1],
+                          '179.2.1.2' => [179, 0, 1],
+                          '181.2.0.1' => [2, 177, 1],
+                          '181.2.1.1' => [181, 0, 1],
+                          '181.2.1.2' => [362, 0, 1],
+                          '191.2.0.1' => [19, 190, 1],
+                          '191.2.1.1' => [1337, 0, 1],
+                          '191.2.1.2' => [191, 0, 1],
+                          '193.2.0.1' => [5, 192, 1],
+                          '193.2.1.1' => [193, 0, 1],
+                          '193.2.1.2' => [965, 0, 1],
+                          '197.2.0.1' => [2, 192, 1],
+                          '197.2.1.1' => [197, 0, 1],
+                          '197.2.1.2' => [394, 0, 1],
+                          '199.2.0.1' => [3, 193, 1],
+                          '199.2.1.1' => [597, 0, 1],
+                          '199.2.1.2' => [199, 0, 1] }.freeze
+      end
+    end
+  end
+end

data/lib/modular_forms/core/padic_pol/padic_pol.rb

@@ -0,0 +1,48 @@
+# frozen_string_literal: true
+
+require_relative './data_pol_2deg'
+
+module ModularForms
+  module Core
+    module PAdicFields
+      # ModularForms::Core::PAdicPol
+      #
+      # Provides a constructor for a hash of defining polynomials of finite extensions of p-adic fields.
+      module PAdicPol
+        def self.build_padic_label_hash(data_padic)
+          padic_coefs_by_label = {}
+          data_padic.each do |label, coefs|
+            p, _y, z, w = label.split('.').map(&:to_i)
+
+            padic_coefs_by_label[p] ||= {}
+            padic_coefs_by_label[p][z] ||= {}
+            padic_coefs_by_label[p][z][w] = coefs
+          end
+          padic_coefs_by_label
+        end
+
+        def self.gen_quadratic_pol_string(coefs)
+          c, b, = coefs
+          terms = []
+          terms << 'x**2'
+
+          if (b != 0) && (b == 1)
+            terms << '+ x'
+          elsif b > 1
+            terms << "+ #{b}*x"
+          end
+
+          terms << "+ #{c}"
+
+          terms.join(' ')
+        end
+
+        def self.get_defining_polynomial(hash, p, z, w)
+          gen_quadratic_pol_string(hash.dig(p, z, w))
+        end
+
+        DEF_POL_HASH_2DEG = build_padic_label_hash(DATA_POL_2DEG)
+      end
+    end
+  end
+end

data/lib/modular_forms/core.rb

@@ -11,3 +11,5 @@ require_relative 'core/sl2z_groups'
 require_relative 'core/elliptic_curves_fp'
 require_relative 'core/dirichlet_characters'
 require_relative 'core/newform_invariants'
+require_relative 'core/l_functions'
+require_relative 'core/padic_fields'

data/lib/modular_forms/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module ModularForms
-  VERSION = '0.0.1'
+  VERSION = '0.0.3'
 end

data/lib/modular_forms.rb

@@ -192,4 +192,24 @@ module ModularForms # rubocop:disable Metrics/ModuleLength
   def analytic_conductor(level_n, weight_k)
     Core::NewFormInvariants.analytic_conductor(level_n, weight_k)
   end
+
+  def a_p(p, cardinality)
+    Core::LFunctions.a_p(p, cardinality)
+  end
+
+  def padic_valuation(num_b10, p)
+    Core::PAdicFields.padic_valuation(num_b10, p)
+  end
+
+  def padic_norm(num_b10, p)
+    Core::PAdicFields.padic_norm(num_b10, p)
+  end
+
+  def padic_expansion(num_b10, p, prec = 11, reverse_trim = false) # rubocop:disable Style/OptionalBooleanParameter
+    Core::PAdicFields.padic_expansion(num_b10, p, prec, reverse_trim)
+  end
+
+  def def_pol_2deg(p = 2, c = 0, num = 1)
+    Core::PAdicFields.def_pol_2deg(p, c, num)
+  end
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: modular_forms
 version: !ruby/object:Gem::Version
-  version: 0.0.1
+  version: 0.0.3
 platform: ruby
 authors:
 - Edgar Armando Delgado Vega
@@ -17,8 +17,10 @@ email:
 executables: []
 extensions: []
 extra_rdoc_files:
+- LICENSE
 - README.md
 files:
+- LICENSE
 - README.md
 - lib/modular_forms.rb
 - lib/modular_forms/core.rb
@@ -29,8 +31,12 @@ files:
 - lib/modular_forms/core/elliptic_curves_q.rb
 - lib/modular_forms/core/hecke_operators.rb
 - lib/modular_forms/core/klein_j_invariant.rb
+- lib/modular_forms/core/l_functions.rb
 - lib/modular_forms/core/newform_invariants.rb
 - lib/modular_forms/core/numeric_helpers/numeric_helpers.rb
+- lib/modular_forms/core/padic_fields.rb
+- lib/modular_forms/core/padic_pol/data_pol_2deg.rb
+- lib/modular_forms/core/padic_pol/padic_pol.rb
 - lib/modular_forms/core/ramanujan_tau_function.rb
 - lib/modular_forms/core/sl2z_groups.rb
 - lib/modular_forms/core/theta_functions.rb