相册

25个线程要做完1000个任务

import threading,time,random,os,threadpool
from shutil import copyfile

def dealing(cata):
    print threading.current_thread().name
    print 'entering ',cata
    os.popen('sh DORFTN '+cata).read()
#    os.chdir('..')

if __name__ == "__main__":
    cata=[c for c in os.listdir('.') if os.path.isdir(c)]
   pool=threadpool.ThreadPool(25)
   requests=threadpool.makeRequests(dealing,cata)
    [pool.putRequest(req) for req in requests]
    pool.wait()


如果函数内有多个参数，可如下处理：
import threadpool, os, glob

out=[]

def processing(sample,trace):
   tmp=sample.split('/')
   tracefile=trace+'/'+tmp[-1]
   ret=os.popen('./doublet_search %s %s' % (sample,tracefile)).read()
    print(ret)

if __name__=='__main__':
   traces=glob.glob('Trace/2*')
   samples=glob.glob('Template/2*')
    arglist=[]
    for sample in samples:
       Zs=glob.glob(sample+'/*.?HZ')
       for Z in Zs:
          sam=Z[0:-1]
           for trace in traces:
             arglist.append(([sam,trace],None))
#    print(arglist)
   pool_t=threadpool.ThreadPool(8)
   requests=threadpool.makeRequests(processing,arglist)
    [pool_t.putRequest(req) for req in requests]
    pool_t.wait()

****注：一般在函数中不要用到os.chdir，因为chdir是全局的，如果要用，则要用绝对路径。

a[:,2:4]=10

a[a>5]=10

np.where(a>5,a,10) # if a>5, values equal to a themselves, else equal to 10

a.clip(3,10) # if a<3, =3, if a<10, =10

sudo vim /etc/ssh/ssh_config

ServiceAliveInterval 20

ServiceAliveCountMax 3

import matplotlib as mpl
from matplotlib.backends.backend_pdf import PdfPages

mpl.rcParams['pdf.fonttype']=42
mpl.rcParams['font.family']='WenQuanYi Zen Hei'

with PdfPages('test.pdf') as pdf:
xxxx # plot some figure
pdf.savefig()

#questions/45870858/runtimeerror-failed-to-open-truetype-font-in-matplotlib-backends-pdf-in-python stackoverflow

import pandas as pd
import numpy as np
 
dates=pd.date_range('20180101','20181231',freq='W-WED')

# W-MON：从指定星期几开始算起，每周# 星期几缩写：MON/TUE/WED/THU/FRI/SAT/SUN

#转换日期
dates_array=dates.to_pydatetime()
dates_list=np.vectorize(lambda x: x.strftime('%Y-%m-%d'))(dates_array).

fig,ax=plt.subplots(figsize=(8,8))
cmap=plt.get_cmap('seismic')
cmap.set_bad(color='gray')
cb=ax.matshow(matdiff,cmap=cmap)
fig.colorbar(cb,ax=ax)
cb.set_clim(vmin=-1.0,vmax=1.0)
plt.show()

方法一：

import numpy as np
a = np.array([0,1,2,3,4,5,5,6,7,8,9])
print(a)
print(a.sum())
# [0 1 2 3 4 5 5 6 7 8 9]
# 50
a = np.ma.array(a, mask=False)
a.mask[3] = True
print(a)
print(a.sum())
# [0 1 2 -- 4 5 5 6 7 8 9]
# 47
方法二
a_new = np.delete(a, 3, 0)
这个3可以是个索引列表

想把一列台站名输出成一行，用空格隔开，写了个C。

fgets(buff,1024,fp);
strncpy(temp,buff,strlen(buff)-1);
printf('%s ',temp);

发现当台站名由长变短时，会“继承”上一个台站的后面几个字符，比如前一个台站名为L5317，而后一个为Y02时，前一个正确，而后一个会变成Y0217.怀疑是temp保存了上一次的’\0’的位置，所以在每次printf后强行temp[0]=’\0’，发现没用，考虑是上次的字符串终止符还在。
我试了以下几种方法，都行：

1 每次printf后，temp[strlen(a)-1]=’0’；

2 用malloc给temp分配内存，读到了新的，就realloc一次，记得最后free掉；   temp=(char *)realloc(temp,sizeof(char)*(strlen(buff)-1))

3 直接把buff数组的最后一个值换成’\0’，然后输出buff。 buff[strlen(buff)-1]=’\0′

某列中的数据要能被完整parse为时间值，可以用dateutil.parse中的parse来测试一下，不足的位数要补齐。比如一个文件temp.txt包含两列数，一列时间，一列观测数据200101011001 -3.9

import pandas as pd

import matplotlib.pyplot as plt

dat=pd.read_csv(‘temp.txt’,head=0,sep=’\s+’,names=[‘time’,’values’],parse_dates=[‘time’])
dat.plot(x=’time’,y=’values’)

plt.show()

https://support.apple.com/zh-cn/HT201341

Mac 不再随附 X11，但 XQuartz 项目会提供 X11 服务器和客户端库。

Apple 创建了 XQuartz 项目，共同致力于在 Mac 上进一步开发和支持 X11。XQuartz 项目最初基于 Mac OS X v10.5 中随附的 X11 版本。目前有多个版本的 XQuartz，其中包含修复程序、针对新功能的支持以及 X11 体验的更多改进功能。Apple 是 XQuartz 项目的贡献者，并致力于确保 X11 与 macOS 和最新版本的 XQuartz 正常配合使用。 

XQuartz 项目提供适用于 MacOS 的 X11 服务器和客户端库，网址是 www.xquartz.org。下载可用的最新版本。

////////////////////////////////////////////////////////////////////////////////////////////

专业软件用到了x11，但编译的时候报错说是没有x11，按照以上通知，安装了XQuartz后，程序运行正常。

brew install xquartz