示例#1
0
 def __init__ (self, r_inter, r_intra, kF, mlist):
    Qm_inter = rpakernel.kernel_m_inter(r_inter.r, mlist, r_inter.label)
    self.Qm_inter = Kernel(r_inter.r, Qm_inter)
    Qm_intra = np.zeros((len(r_intra.r), len(r_intra.r)))
    if (kF * r_intra.r.max() > 0.01):
        Qm_intra = rpam.kernel_m_intra(r_intra.r, mlist, kF, '3' + r_intra.label) 
        #rpakernel.kernel_m_intra(r, mlist, kF)
    self.Qm_intra = Kernel(r_intra.r, Qm_intra)
示例#2
0
 def __init__(self, r_inter, r_intra, kF, mlist):
     Qm_inter = rpakernel.kernel_m_inter(r_inter.r, mlist, r_inter.label)
     self.Qm_inter = Kernel(r_inter.r, Qm_inter)
     Qm_intra = np.zeros((len(r_intra.r), len(r_intra.r)))
     if (kF * r_intra.r.max() > 0.01):
         Qm_intra = rpam.kernel_m_intra(r_intra.r, mlist, kF,
                                        '3' + r_intra.label)
         #rpakernel.kernel_m_intra(r, mlist, kF)
     self.Qm_intra = Kernel(r_intra.r, Qm_intra)
示例#3
0
import rpam
import rpakernel

import numpy as np
r = np.linspace(0.01, 50.0, 500)
m = np.arange(0.0, 31.0, 1.0)
print r
print m
kF = 3.0

Q = rpam.kernel_m_intra(r, m, kF, '3')

#Q1 = rpakernel.kernel_m_intra(r, m, kF)
Q1 = rpam.kernel_m_intra(r, m, kF, '4')

import pylab as pl

for r_i in [0.1, 1.0, 2.0, 5.0, 7.0]:
    i = np.abs(r - r_i).argmin()
    pl.figure()
    pl.plot(r, Q[i, :], label='Q0: r = %g' % r[i])
    pl.plot(r, Q1[i, :], label='Q1: r = %g' % r[i])
    pl.plot(r, Q1[i, :] - Q[i, :], label='diff: r = %g' % r[i])
    pl.legend()
pl.show()