modular_forms 0.0.3 → 0.0.5

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f9b49161da093d436c9afd089b0d2935a43fafb71d1b1eece1af75cc3ef08833
-  data.tar.gz: 54672d7ef27bc439d12747a919a024e9c9fa58d76a3d4b3ae1fa3f807226be6b
+  metadata.gz: 3db834376767beecfe9ac15dba1a8ac97e3e9d845db19729668bbd32d511278d
+  data.tar.gz: d12f3d74498c06a7a16f264ab2d659b372e5629d2688470a4f39d2c3ad71883d
 SHA512:
-  metadata.gz: d8ad361a6edc83f140ac815bb4b04a6a627100204f523e2192edfaa4cdc8e7b1022cd1c4e6651526582d287818231193d546dd2fb99b859d39401cb5002b3861
-  data.tar.gz: 2a7223689108395ebc2f389d1e1ae9b6e61a5165e9da9885624b7758500131ca01f983c433dc045928c433608d773b6ea22410d2196b43b2261cb88c40933b27
+  metadata.gz: a599f05653f0abea5792617200ffba19726ee0abe421de1e04b13a4beb92e3f1911e61688ee05b55aae7aa8d8114fa63378e188fceb91a568160e20307c77232
+  data.tar.gz: bddeae428139325379f8172ebac90d324da39c633470823b40fbebbd0211358c7805dcb4db0c464e1e5156619e59a912d988a28006a9fbc9e08af29cbc4ee856

data/README.md

@@ -3,32 +3,31 @@
 A creative toolkit for exploring modular forms and elliptic curves through [Sonic Pi](https://sonic-pi.net/).
 
 <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">
+  <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" />
 </p>
 
-[![Modular Forms - Image](https://raw.githubusercontent.com/edelveart/modular_forms/main/modular_forms.png)](https://rubygems.org/gems/modular_forms)
+<p align="center">
+  <img src="https://raw.githubusercontent.com/edelveart/modular_forms/main/modular_forms.png" width="400" alt="Modular Forms Image" />
+</p>
 
-## ⚙️ Project Status
+## 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.
+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. 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.
+Future updates might include a **DSL**, depending on how the library is used and the interest from the community. A key challenge is turning mathematical ideas into music that feels deliberate, expressive, and naturally musical.
 
-## 🧊 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
-
-### Overview
+##  Scope
 
 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:
+- List of implemented modules:
   - [Eisenstein Series](#eisenstein-series)
   - [Eta Functions and Eta Quotients](#eta-functions-and-eta-quotients)
   - [Theta Functions](#theta-functions)
@@ -42,16 +41,14 @@ Given the vastness of the field, this tool intentionally focuses on a limited su
   - [Newform Invariants](#newform-invariants)
   - [L-functions](#l-functions)
   - [p-adic Fields](#p-adic-fields)
+  - [Number Fields](#number-fields)
 
 ### Not Optimized for Computational Efficiency
 
-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**.
+> 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.
 
@@ -76,9 +73,9 @@ Then require the file like this:
 require "<PATH>/modular_forms.rb"
 ```
 
-## 🧰 How to use?
+## How to use?
 
-### 🎶 Explore Fermat’s Last Theorem in Sonic Pi
+### 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**.
 
@@ -115,24 +112,36 @@ live_loop :modularity_music do
 end
 ```
 
-### 🎶 Merging Machine and Organism - Warp Up
+### 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)
+eisenstein_melody = ModularForms.eisenstein_series(8)
 e8 = eisenstein_melody.take(58)
 hecke_op = ModularForms.hecke_operator_prime_non_cusp(e8, p, 8, 20)
+
+ellc = ModularForms.elliptic_curve_q([2, 5])
+disc = ModularForms.discriminant_q(ellc)
+
 j_func = ModularForms.j_function(40)
 newform_ac = ModularForms.analytic_conductor(15, 2)
 pol = ModularForms.def_pol_2deg(41)
 
+x = 7
+dirichlet_char_group = []
+(1..x).each do |i|
+  dirichlet_char_group << ModularForms.conrey_p_pminus1(x, i)
+end
+
 (disc * -1).times do
+  synth :zawa,
+    note: ModularForms.zeta_coeffs_deg2(dirichlet_char_group, 30)
+  .tick(:zeta) + 72, amp: rrand(0.3, 0.6)
+
   synth :subpulse, note: ModularForms.padic_valuation(eisenstein_melody.next, p) % 7 + 70,
     release: 0.25 if (spread (disc % 6), 7).tick(:d)
 
@@ -152,7 +161,7 @@ pol = ModularForms.def_pol_2deg(41)
 end
 ```
 
-## 🧩 Implemented Modular Forms, Elliptic Curves, and Related Definitions
+## Implemented Modular Forms, Elliptic Curves, and Related Definitions
 
 ### Eisenstein Series
 
@@ -209,41 +218,46 @@ end
 27. `ModularForms.discriminant_q(curve)`
 28. `ModularForms.j_invariant_q(curve)`
 29. `ModularForms.point_on_curve_q?(curve, point)`
-30. `ModularForms.point_addition_q(curve, p, q)`
-31. `ModularForms.scalar_mul_point_q(curve, n, point)`
-32. `ModularForms.isogeny_2deg_q(curve, point_2tor)`
-33. `ModularForms.isogeny_ndeg_q(curve, point_ntor, order)`
-34. `ModularForms.weil_height(x_point)`
-35. `ModularForms.canonical_height(curve, point, prec = 64)`
+30. `ModularForms.negate_p(point)`
+31. `ModularForms.point_addition_q(curve, p, q)`
+32. `ModularForms.scalar_mul_point_q(curve, n, point)`
+33. `ModularForms.isogeny_2deg_q(curve, point_2tor)`
+34. `ModularForms.isogeny_ndeg_q(curve, point_ntor, order)`
+35. `ModularForms.weil_height(x_point)`
+36. `ModularForms.canonical_height(curve, point, prec = 64)`
 
 ### Elliptic Curves over Finite Fields
 
-36. `ModularForms.elliptic_curve_fp(p, coefs)`
-37. `ModularForms.point_on_curve_modp?(curve, point)`
-38. `ModularForms.discriminant_modp(curve)`
-39. `ModularForms.j_invariant_modp(curve)`
-40. `ModularForms.point_addition_modp(curve, p_point, q_point)`
-41. `ModularForms.scalar_mul_point_modp(curve, n, point)`
-42. `ModularForms.points_fp(curve, point_at_infinity = false)`
-43. `ModularForms.cardinality_fp(curve)`
-44. `ModularForms.quadratic_twist_fp(curve)`
+37. `ModularForms.elliptic_curve_fp(p, coefs)`
+38. `ModularForms.point_on_curve_modp?(curve, point)`
+39. `ModularForms.discriminant_modp(curve)`
+40. `ModularForms.j_invariant_modp(curve)`
+41. `ModularForms.point_addition_modp(curve, p_point, q_point)`
+42. `ModularForms.scalar_mul_point_modp(curve, n, point)`
+43. `ModularForms.points_fp(curve, point_at_infinity = false)`
+44. `ModularForms.cardinality_fp(curve)`
+45. `ModularForms.quadratic_twist_fp(curve)`
 
 ### Newform Invariants
 
-45. `ModularForms.analytic_conductor(level_n, weight_k)`
+46. `ModularForms.analytic_conductor(level_n, weight_k)`
 
 ### L-functions
 
-46. `ModularForms.a_p(p, cardinality)`
+47. `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)`
+48. `ModularForms.padic_valuation(num_b10, p)`
+49. `ModularForms.padic_norm(num_b10, p)`
+50. `ModularForms.padic_expansion(num_b10, p, prec = 11, reverse_trim = false)`
+51. `ModularForms.def_pol_2deg(p = 2, c = 0, num = 1)`
+
+### Number Fields
+
+52. `ModularForms.zeta_coeffs_deg2(dirichlet_char_group, n)`
 
-## 🧪 Testing
+## Testing
 
 Install dependencies first:
 

data/lib/modular_forms/core/elliptic_curves_fp.rb

@@ -13,7 +13,7 @@ module ModularForms
         ((n % modp) + modp) % modp
       end
 
-      def self.elliptic_curve_fp(p, coefs)
+      def self.elliptic_curve_fp(p, coefs) # rubocop:disable Metrics/MethodLength
         a, b = coefs
         raise "#{p} is not a prime number" if NumericHelpers.prime_number?(p) == false
 
@@ -21,9 +21,13 @@ module ModularForms
 
         a_modp = reduction_modp(a, p)
         b_modp = reduction_modp(b, p)
-        raise "y^2=x^3 #{a_modp}x #{b_modp} defines a singular curve" if d == 0 # rubocop:disable Style/NumericPredicate
 
-        puts "y^2 = x^3 #{a_modp}x #{b_modp} over Finite Field #{p}"
+        if d.zero?
+          raise "#{EllipticCurvesQ.show_curve_equation(a_modp,
+                                                       b_modp)} defines a singular curve over F_#{p}"
+        end
+
+        puts "#{EllipticCurvesQ.show_curve_equation(a_modp, b_modp)} over Finite Field #{p}"
         { a: a, b: b, p: p }
       end
 

data/lib/modular_forms/core/elliptic_curves_q.rb

@@ -7,16 +7,35 @@ module ModularForms
     # ModularForms::Core::EllipticCurves
     #
     # This module provides methods for generating points on Elliptic Curves (short Weierstrass form) over Q
-    module EllipticCurvesQ
+    module EllipticCurvesQ # rubocop:disable Metrics/ModuleLength
       def self.discriminant(a, b)
         -16 * (4 * a**3 + 27 * b**2)
       end
 
+      def self.show_curve_equation(a, b) # rubocop:disable Metrics/MethodLength,Lint/RedundantCopDisableDirective
+        eq = 'y^2 = x^3'
+        eq += a != 0 ? " #{a > 0 ? '+' : '-'} #{a.abs}x" : '' # rubocop:disable Style/NestedTernaryOperator,Style/NumericPredicate
+        eq += b != 0 ? " #{b > 0 ? '+' : '-'} #{b.abs}" : '' # rubocop:disable Style/NestedTernaryOperator,Style/NumericPredicate
+        eq
+      end
+
+      def self.singular?(a, b)
+        delta = discriminant(a, b)
+        return unless delta.zero?
+
+        type = begin
+          x_s = a.zero? ? 0 : Math.sqrt(-a / 3) rescue 0 # rubocop:disable Style/RescueModifier
+          x_s.zero? ? 'cusp' : 'node'
+        end
+
+        # Concatenar la ecuación directamente en el mensaje
+        raise "#{show_curve_equation(a, b)} defines a singular curve (#{type})"
+      end
+
       def self.elliptic_curve_q(coefs)
         a, b = coefs
-        raise "y^2=x^3 #{a}x #{b} defines a singular curve" if discriminant(a, b) == 0 # rubocop:disable Style/NumericPredicate
-
-        puts "y^2 = x^3 #{a}x #{b}"
+        singular?(a, b)
+        puts show_curve_equation(a, b)
         { a: a, b: b }
       end
 
@@ -126,3 +145,6 @@ module ModularForms
     end
   end
 end
+# # [-1, 2],
+# ModularForms::Core::EllipticCurvesQ.elliptic_curve_q([-27, 54])
+# # ModularForms::Core::EllipticCurvesQ.show_curve_equation([-0, 0])

data/lib/modular_forms/core/number_fields.rb

@@ -0,0 +1,31 @@
+# frozen_string_literal: true
+
+module ModularForms
+  module Core
+    # ModularForms::Core::NumberFields
+    #
+    # This module provides methods and objects related to number fields,
+    # such as zeta_coefficients and utility helpers.
+    module NumberFields
+      def self.fill_arr_with_interval(length, num, interval)
+        arr = Array.new(length, 0)
+        (interval - 1).step(length - 1, interval) do |i|
+          arr[i] = num
+        end
+        arr
+      end
+
+      def self.zeta_coeffs_deg2(dirichlet_char_group, num_coeffs)
+        ded_zeta_to_dir_serie = []
+        dirichlet_char_group.each_with_index do |i, j|
+          a = 1
+          while a <= num_coeffs
+            ded_zeta_to_dir_serie << fill_arr_with_interval(num_coeffs, i, j + a)
+            a += dirichlet_char_group.length
+          end
+        end
+        ded_zeta_to_dir_serie.transpose.map { |component| component.sum } # rubocop:disable Style/SymbolProc
+      end
+    end
+  end
+end

data/lib/modular_forms/core.rb

@@ -13,3 +13,4 @@ require_relative 'core/dirichlet_characters'
 require_relative 'core/newform_invariants'
 require_relative 'core/l_functions'
 require_relative 'core/padic_fields'
+require_relative 'core/number_fields'

data/lib/modular_forms/version.rb

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

data/lib/modular_forms.rb

@@ -212,4 +212,8 @@ module ModularForms # rubocop:disable Metrics/ModuleLength
   def def_pol_2deg(p = 2, c = 0, num = 1)
     Core::PAdicFields.def_pol_2deg(p, c, num)
   end
+
+  def zeta_coeffs_deg2(dirichlet_char_group, n)
+    Core::NumberFields.zeta_coeffs_deg2(dirichlet_char_group, n)
+  end
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: modular_forms
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.0.5
 platform: ruby
 authors:
 - Edgar Armando Delgado Vega
@@ -33,6 +33,7 @@ files:
 - 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/number_fields.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