Skip to content

Instantly share code, notes, and snippets.

View liaocs2008's full-sized avatar
🎯
Focusing

Leo liaocs2008

🎯
Focusing
View GitHub Profile
"""
newA
[[1 2 1]
[1 1 1]
[2 4 2]]
L, U = LU([(3-1,2-1,0),(3-1,2-1,3)], newA)
L= [[ 1. 0. 0.]
[-0. 1. -0.]
[ 0. 3. 1.]]
"""
real part
[[ 2.98023224e-08 0.00000000e+00 2.98023224e-08]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
img part
[[ 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 2.98023224e-08 0.00000000e+00]]
real part
@liaocs2008
liaocs2008 / gist:ac4091d5f1b28bf517623edede76fba9
Last active December 19, 2016 23:00
Understanding in convolutions
import numpy as np
# http://www.songho.ca/dsp/convolution/convolution.html#separable_convolution
# https://www.mathworks.com/help/signal/ug/linear-and-circular-convolution.html
# make circular convolution equivalent to linear convolution
# M = 13
# N = 7