-# Copyright (C) 2008-2010 Alberto Gomez-Casado
+# Copyright (C) 2008-2011 Alberto Gomez-Casado
# Fabrizio Benedetti
# Massimo Sandal <devicerandom@gmail.com>
# W. Trevor King <wking@drexel.edu>
"""
import copy
+import os.path
+import re
import numpy
+import yaml
from ..command import Command, Argument, Failure
from ..command_stack import CommandStack
pass # no need to place duplicate calls on the stack.
else:
curve.command_stack.append(CommandMessage(
- self.name, params))
+ self.name, dict(params)))
class BlockCommand (CurveCommand):
def __init__(self):
super(CurvePlugin, self).__init__(name='curve')
self._commands = [
- GetCommand(self), InfoCommand(self), DeltaCommand(self),
- ExportCommand(self), DifferenceCommand(self),
+ GetCommand(self), InfoCommand(self), BlockInfoCommand(self),
+ DeltaCommand(self), ExportCommand(self), DifferenceCommand(self),
DerivativeCommand(self), PowerSpectrumCommand(self),
ClearStackCommand(self)]
Argument(name='all', type='bool', default=False, count=1,
help='Get all curve information.'),
]
- self.fields = ['name', 'path', 'experiment', 'driver', 'filetype',
- 'note', 'command stack', 'blocks', 'block sizes']
+ self.fields = ['name', 'path', 'driver', 'note', 'command stack',
+ 'blocks', 'block names', 'block sizes']
for field in self.fields:
args.append(Argument(
name=field, type='bool', default=False, count=1,
def _get_path(self, curve):
return curve.path
- def _get_experiment(self, curve):
- return curve.info.get('experiment', None)
-
def _get_driver(self, curve):
return curve.driver
- def _get_filetype(self, curve):
- return curve.info.get('filetype', None)
-
def _get_note(self, curve):
return curve.info.get('note', None)
def _get_blocks(self, curve):
return len(curve.data)
+ def _get_block_names(self, curve):
+ return [block.info['name'] for block in curve.data]
+
def _get_block_sizes(self, curve):
return [block.shape for block in curve.data]
+class BlockInfoCommand (BlockCommand):
+ """Get selected information about a :class:`hooke.curve.Curve` data block.
+ """
+ def __init__(self, plugin):
+ super(BlockInfoCommand, self).__init__(
+ name='block info', arguments=[
+ Argument(
+ name='key', count=-1, optional=False,
+ help='Dot-separted (.) key selection regexp.'),
+ Argument(
+ name='output',
+ help="""
+File name for the output (appended).
+""".strip()),
+ ],
+ help=self.__doc__, plugin=plugin)
+
+ def _run(self, hooke, inqueue, outqueue, params):
+ block = self._block(hooke, params)
+ values = {'index': self._block_index(hooke, params)}
+ for key in params['key']:
+ keys = [(0, key.split('.'), block.info)]
+ while len(keys) > 0:
+ index,key_stack,info = keys.pop(0)
+ regexp = re.compile(key_stack[index])
+ matched = False
+ for k,v in info.items():
+ if regexp.match(k):
+ matched = True
+ new_stack = copy.copy(key_stack)
+ new_stack[index] = k
+ if index+1 == len(key_stack):
+ vals = values
+ for k in new_stack[:-1]:
+ if k not in vals:
+ vals[k] = {}
+ vals = vals[k]
+ vals[new_stack[-1]] = v
+ else:
+ keys.append((index+1, new_stack, v))
+ if matched == False:
+ raise ValueError(
+ 'no match found for %s (%s) in %s'
+ % (key_stack[index], key, sorted(info.keys())))
+ if params['output'] != None:
+ curve = self._curve(hooke, params)
+ with open(params['output'], 'a') as f:
+ yaml.dump({curve.name:{
+ 'path': curve.path,
+ block.info['name']: values
+ }}, f)
+ outqueue.put(values)
+
+
class DeltaCommand (BlockCommand):
"""Get distance information between two points.
def _run(self, hooke, inqueue, outqueue, params):
data = self._block(hooke, params)
- with open(params['output'], 'w') as f:
+ with open(os.path.expanduser(params['output']), 'w') as f:
if params['header'] == True:
f.write('# %s \n' % ('\t'.join(data.info['columns'])))
numpy.savetxt(f, data, delimiter='\t')
def _run(self, hooke, inqueue, outqueue, params):
self._add_to_command_stack(params)
- params = self.__setup_params(hooke=hooke, params=params)
+ params = self._setup_params(hooke=hooke, params=params)
data_A = self._get_column(hooke=hooke, params=params,
block_name='block A',
column_name='column A')
column_name='output column',
values=out)
- def __setup_params(self, hooke, params):
+ def _setup_params(self, hooke, params):
curve = self._curve(hooke, params)
if params['block A'] == None:
params['block A'] = curve.data[0].info['name']
def _run(self, hooke, inqueue, outqueue, params):
self._add_to_command_stack(params)
- params = self.__setup_params(hooke=hooke, params=params)
+ params = self._setup_params(hooke=hooke, params=params)
x_data = self._get_column(hooke=hooke, params=params,
column_name='x column')
f_data = self._get_column(hooke=hooke, params=params,
column_name='output column',
values=d)
- def __setup_params(self, hooke, params):
+ def _setup_params(self, hooke, params):
curve = self._curve(hooke, params)
x_name,x_unit = split_data_label(params['x column'])
f_name,f_unit = split_data_label(params['f column'])
def _run(self, hooke, inqueue, outqueue, params):
self._add_to_command_stack(params)
- params = self.__setup_params(hooke=hooke, params=params)
+ params = self._setup_params(hooke=hooke, params=params)
data = self._get_column(hooke=hooke, params=params)
bounds = params['bounds']
if bounds != None:
values=power)
outqueue.put(b)
- def __setup_params(self, hooke, params):
+ def _setup_params(self, hooke, params):
if params['output block'] in self._block_names(hooke, params):
raise Failure('output block %s already exists in %s.'
% (params['output block'],