modular_forms 0.0.2 → 0.0.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 749de4a6eccae53fa56728d33efae75d5f9c0d4376d2031760202dddb0815666
-  data.tar.gz: 0c4d325ea237d6c6478c2fae75f862e0b1f8e116d1b822b5a35eac5481ea88c4
+  metadata.gz: f9b49161da093d436c9afd089b0d2935a43fafb71d1b1eece1af75cc3ef08833
+  data.tar.gz: 54672d7ef27bc439d12747a919a024e9c9fa58d76a3d4b3ae1fa3f807226be6b
 SHA512:
-  metadata.gz: 7f9b7939543427e63283e38d4c3c4bd796823962809c65348ef0023c6a66a94cde414bddf5cfffe2f2ed0c88e342f90da8aa60774de78cf932d09d36a89f8d18
-  data.tar.gz: cc6ea4af26b7e9502d607e6128f522a383e03f4da9c1cf2b328158fb94f6b85a226cb3c2f9837fc00fbc2a1f0a1370d58b88df7c7a930f1e3968d03235829949
+  metadata.gz: d8ad361a6edc83f140ac815bb4b04a6a627100204f523e2192edfaa4cdc8e7b1022cd1c4e6651526582d287818231193d546dd2fb99b859d39401cb5002b3861
+  data.tar.gz: 2a7223689108395ebc2f389d1e1ae9b6e61a5165e9da9885624b7758500131ca01f983c433dc045928c433608d773b6ea22410d2196b43b2261cb88c40933b27

data/README.md

@@ -14,37 +14,41 @@ A creative toolkit for exploring modular forms and elliptic curves through [Soni
 ## ⚙️ 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 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
 
 - **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
 
-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:
+### Overview
 
-- [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)
-- [L-functions](#l-functions)
+Given the vastness of the field, this tool intentionally focuses on a limited subset of definitions, without covering all aspects of each.
+
+- 🧩 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
@@ -60,31 +64,20 @@ If you are **use Ruby**, then import via
 require 'modular_forms'
 ```
 
-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:
+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:
 
-```rb
-require "<PATH>/modular_forms.rb"
+```bash
+gem which modular_forms
 ```
 
-## 🧰 How to use?
-
-### 🎶 Play an Eisenstein Series as Music for Warm-Up
-
-Here is a simple example of how to use **modular_forms** to generate a basic musical pattern:
+Then require the file like this:
 
 ```rb
 require "<PATH>/modular_forms.rb"
-
-# Generate an Eisenstein series of weight k = 4
-eisenstein_melody = ModularForms.eisenstein_serie(4)
-
-# Play the melody in a loop with a mathematical transformation
-120.times do
-  play eisenstein_melody.next % 12 * 7
-  sleep 0.5
-end
 ```
 
+## 🧰 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**.
@@ -122,6 +115,43 @@ live_loop :modularity_music do
 end
 ```
 
+### 🎶 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
@@ -138,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)`
@@ -198,10 +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/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

@@ -12,3 +12,4 @@ 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.2'
+  VERSION = '0.0.3'
 end

data/lib/modular_forms.rb

@@ -196,4 +196,20 @@ module ModularForms # rubocop:disable Metrics/ModuleLength
   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.2
+  version: 0.0.3
 platform: ruby
 authors:
 - Edgar Armando Delgado Vega
@@ -34,6 +34,9 @@ files:
 - 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