| Server IP : 194.181.228.65 / Your IP : 216.73.216.135 Web Server : LiteSpeed System : Linux wn13.webd.pl 5.14.0-687.46.1.el9_8.x86_64 #1 SMP PREEMPT_DYNAMIC Fri Sep 11 09:03:19 EDT 2026 x86_64 User : prusit ( 133039) PHP Version : 5.6.40 Disable Function : exec, system, shell_exec, passthru, proc_open, proc_close, popen MySQL : ON | cURL : ON | WGET : ON | Perl : ON | Python : OFF | Sudo : ON | Pkexec : ON Directory : /lib64/python3.9/site-packages/scipy/signal/ |
Upload File : |
# Author: Eric Larson
# 2014
import numpy as np
#pythran export _max_len_seq_inner(intp[], int8[], int, int, int8[])
#pythran export _max_len_seq_inner(int[], int8[], int, int, int8[])
# Fast inner loop of max_len_seq.
def _max_len_seq_inner(taps, state, nbits, length, seq):
# Here we compute MLS using a shift register, indexed using a ring buffer
# technique (faster than using something like np.roll to shift)
n_taps = taps.shape[0]
idx = 0
for i in range(length):
feedback = state[idx]
seq[i] = feedback
for ti in range(n_taps):
feedback ^= state[(taps[ti] + idx) % nbits]
state[idx] = feedback
idx = (idx + 1) % nbits
# state must be rolled s.t. next run, when idx==0, it's in the right place
return np.roll(state, -idx, axis=0)