# -*- coding: utf-8 -*-
"""Modify the PhysicalQuantities class for more convenient interactive use.

Also redefine some math functions to operate on PhysQties with no need for
special method syntax. This just means moving them out to the global
namespace.

This module should always be loaded *after* math or Numeric, so it can
overwrite math functions with the versions that handle units.

Authors
-------
- Fernando Perez <Fernando.Perez@berkeley.edu>
"""

#*****************************************************************************
#       Copyright (C) 2008-2009 The IPython Development Team
#       Copyright (C) 2001-2007 Fernando Perez <fperez@colorado.edu>
#
#  Distributed under the terms of the BSD License.  The full license is in
#  the file COPYING, distributed as part of this software.
#*****************************************************************************

from Scientific.Physics.PhysicalQuantities import PhysicalQuantity

# This code can be set up to work with Numeric or with math for providing the
# mathematical functions. Uncomment the one you prefer to use below.

# If you use math, sin(x) won't work for x an array, only float or PhysQty
import math

# If you use Numeric, sin(x) works for x a float, PhysQty an array.
#import Numeric as math

class PhysicalQuantityFunction:
    """Generic function wrapper for PhysicalQuantity instances.

    Calls functions from either the math library or the instance's methods as
    required.  Allows using sin(theta) or sqrt(v**2) syntax irrespective of
    whether theta is a pure number or a PhysicalQuantity.

    This is *slow*. It's meant for convenient interactive use, not for
    speed."""

    def __init__(self,name):
        self.name = name
        
    def __call__(self,x):
        if isinstance(x,PhysicalQuantity):
            return PhysicalQuantity.__dict__[self.name](x)
        else:
            return math.__dict__[self.name](x)

class PhysicalQuantityInteractive(PhysicalQuantity):
    """Physical quantity with units - modified for Interactive use.

    Basically, the __str__ and __repr__ methods have been swapped for more
    convenient interactive use. Powers are shown as ^ instead of ** and only 4
    significant figures are shown.

    Also adds the following aliases for commonly used methods:
      b = PhysicalQuantity.inBaseUnits
      u = PhysicalQuantity.inUnitsOf
      
    These are useful when doing a lot of interactive calculations.
    """
    
    # shorthands for the most useful unit conversions
    b = PhysicalQuantity.inBaseUnits  # so you can just type x.b to get base units
    u = PhysicalQuantity.inUnitsOf

    # This can be done, but it can get dangerous when coupled with IPython's
    # auto-calling. Everything ends up shown in baseunits and things like x*2
    # get automatically converted to k(*2), which doesn't work.
    # Probably not a good idea in general...
    #__call__ = b
    
    def __str__(self):
        return PhysicalQuantity.__repr__(self)
 
    def __repr__(self):
        value = '%.4G' % self.value
        units = self.unit.name().replace('**','^')
        return value + ' ' + units

# implement the methods defined in PhysicalQuantity as PhysicalQuantityFunctions
sin = PhysicalQuantityFunction('sin')
cos = PhysicalQuantityFunction('cos')
tan = PhysicalQuantityFunction('tan')
sqrt = PhysicalQuantityFunction('sqrt')