blob: ec3eebd46624f6a6ef0a9051a0ed135b16c7ee2e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
#! /usr/bin/python
from Numeric import *
import sys
import re
recpat = re.compile(" *([\d\.E+-]+) +([\d\.E+-]+) +([\d\.E+-]+)")
f=open(sys.argv[1])
card = f.readline()
m=recpat.match(card)
v=array(map(float,(m.group(1),m.group(2), m.group(3))))
while 1:
card = f.readline()
if len(card)<3: break
m=recpat.match(card)
w=array(map(float,(m.group(1),m.group(2), m.group(3))))
print (w-v)[0], (w-v)[1], (w-v)[2]
v=w
|