• പൈത്തൺ
  • റസ്റ്റ്
  • ഹാസ്കെൽ
  • സി
  • പ്രകടനം
  • പൈത്തൺ പരിഹരിക്കുന്നു
  • ഹോം
  • പോസ്റ്റുകൾ
  • കുറിപ്പുകൾ
  • പുസ്തകശാല
  • രചയിതാവ്
🇺🇸 en 🇫🇷 fr 🇨🇳 zh

Nathaniel Thomas

ബാസൽ പ്രശ്നം (ഹലോ, ലോകം!)

2023, ജൂലൈ 28

ഹലോ, ലോകം! ഇതാണ് എന്റെ ആദ്യത്തെ പോസ്റ്റ്, ഈ വെബ്സൈറ്റിന്റെ പ്രവർത്തനം പരിശോധിക്കാൻ മാത്രമാണ് ഇത് ഉപയോഗിക്കുന്നത്.

ബാസൽ പ്രശ്നം കണക്കാക്കുന്ന വിവിധ ഭാഷകളിലെ ചില കോഡ് സ്നിപ്പറ്റുകൾ ഇതാ:

ആദ്യം, LATE​X

n=1∑∞​n21​=6π2​=1.6449340668482264

പൈത്തൺ

def pi_squared_over_6(N: int) -> float:
    return sum(x**(-2) for x in range(1,N))

റസ്റ്റ്

fn pi_squared_over_6(N: u64) -> f64 {
    (1..N).map(|x| 1.0 / ((x*x) as f64)).sum()
}

ഹാസ്കെൽ

piSquaredOver6 :: Integer -> Double
-- no capital N in Haskell :(
piSquaredOver6 n = sum $ map (\x -> 1 / fromIntegral (x * x)) [1..n]

സി

double pi_squared_over_6(unsigned int N) {
    double sum = 0.0;
    for (int i = 1; i < N; i++) {
        sum += 1.0 / (i*i);
    }
    return sum;
}

നിങ്ങളുടെ പ്രിയപ്പെട്ട പരിഹാരം ഏതാണ്?

പ്രകടനം

ഒരു എം1 പ്രോയിൽ N=109 എന്നതിന് അവ എങ്ങനെ പ്രകടനം കാണിക്കുന്നുവെന്ന് നോക്കാം.

ഭാഷ സമയം (ms, μ±σ)
റസ്റ്റ് (സമാന്തരം) 112.6±3.5
റസ്റ്റ് (–release) 937.9±0.4
സി (-O3) 995.3±0.8
ഹാസ്കെൽ (-O3) 13454±205
പൈത്തൺ (3.10) 67720±0

പൈത്തൺ പരിഹരിക്കുന്നു

പൈത്തൺ കോഡ് രാവണികളായ സമയം എടുത്തു, അതിനാൽ വെക്റ്ററൈസ്ഡ് സി കോഡിലേക്ക് വിളിക്കുന്ന numpy ഉപയോഗിച്ച് ഇത് പരിഹരിക്കാം.

import numpy as np

def pi_squared_over_6(N: int) -> float:
    x = np.ones(N)
    r = np.arange(1,N)
    sq = np.square(r)
    div = np.divide(x, sq)
    return float(np.sum(div))

കുറച്ച് നല്ലത്, പക്ഷേ btm പരിശോധിക്കുമ്പോൾ, അമിതമായ മെമ്മറി ഉപയോഗം സൂചിപ്പിക്കുന്നത് ബില്യണുകളായ ഫ്ലോട്ടുകൾ ചലിപ്പിക്കുന്നതാണ് പ്രധാന പ്രവർത്തനം, യഥാർത്ഥ ഗണിതമല്ല എന്നാണ്. ഇത് ചങ്കുകളായി വിഭജിക്കാൻ ശ്രമിക്കാം:

def pi_squared_over_6(N: int) -> float:
    CHUNKS = 25000
    SIZE = N // CHUNKS
    s = 0.0
    x = np.ones(N // CHUNKS - 1)
    for i in range(CHUNKS):
        N_tmp = i * SIZE
        r = np.arange(N_tmp + 1, N_tmp + SIZE)
        sq = np.square(r)
        div = np.divide(x, sq)
        s += np.sum(div)
        # മെമ്മറി ഡീലോക്കേറ്റ് ചെയ്യുക
        del sq
        del div
        del r
        
    return s

വളരെ നല്ലത്! ഇപ്പോൾ ഇത് 2 സെക്കൻഡിനുള്ളിൽ പ്രവർത്തിക്കുന്നു!


പൈത്തണിലെ അത്രയും സാധാരണമല്ലാത്ത പ്രകടന ഒപ്റ്റിമൈസേഷൻ
→

back to top