<commit_before>#-*- coding:Utf-8 -*-
from pylayers.simul.link import *
import numpy as np
import matplotlib.pyplot as plt 
import time
print "======================="
print " start test_reciprocity.py "
print "======================="

DL=DLink(L='defstr.ini')

DL.a=np.array([759,1114,1.0])
DL.b=np.array([767,1114,1.5])
DL.fGHz=np.arange(2,11,0.1)
DL.wav = wvf.Waveform(fcGHz=5,bandGHz=3)
# DL.eval(diffraction=True,force=True,ra_vectorized=True,alg=20152,si_reverb=2)
DL.eval(diffraction=True,force=True,ra_vectorized=False,alg=5,cutoff=2,applywav=True)

DL2=DLink(L='defstr.ini')

DL2.a=np.array([759,1114,1.0])
DL2.b=np.array([767,1114,1.5])
DL2.fGHz=np.arange(2,11,0.1)
DL2.wav = wvf.Waveform(fcGHz=5,bandGHz=3)
# ##########################################################################
# #!!never use the following part except for this particular test file !!!!
# ##########################################################################
DL2.load(DL2.Si,DL.dexist['sig']['grpname'])
r = DL2.Si.raysv(DL2.a,DL2.b)
rr=r.reciprocal()
import ipdb
ipdb.set_trace()<commit_msg>fix reciprocity test<commit_after>#-*- coding:Utf-8 -*-
from pylayers.simul.link import *
import numpy as np
import matplotlib.pyplot as plt 
import time
print "======================="
print " start test_reciprocity.py "
print "======================="

DL=DLink(L='defstr.ini')

DL.a=np.array([759,1114,1.0])
DL.b=np.array([767,1114,1.5])
DL.fGHz=np.arange(2,11,0.1)
DL.wav = wvf.Waveform(fcGHz=5,bandGHz=3)
# DL.eval(diffraction=True,force=True,ra_vectorized=True,alg=20152,si_reverb=2)
ra_vectorized = False

DL.eval(diffraction=True,force=True,ra_vectorized=ra_vectorized,alg=5,cutoff=2,applywav=True)
DL2=DLink(L='defstr.ini')

DL2.a=np.array([759,1114,1.0])
DL2.b=np.array([767,1114,1.5])
DL2.fGHz=np.arange(2,11,0.1)
DL2.wav = wvf.Waveform(fcGHz=5,bandGHz=3)
# ##########################################################################
# #!!never use the following part except for this particular test file !!!!
# ##########################################################################
DL2.load(DL2.Si,DL.dexist['sig']['grpname'])
if ra_vectorized:
    r = DL2.Si.raysv(DL2.a,DL2.b)
else:
    r = DL2.Si.rays(DL2.a,DL2.b)
rr=r.reciprocal()
{'fraction_non_alphanumeric': 0.3343792021515016, 'fraction_numerical': 0.0632003585835948, 'mean_word_length': 7.857142857142857, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '5515908', 'n_tokens_mistral': 966, 'n_tokens_neox': 881, 'n_words': 109}