DMelt:Math/1 Functions

From HandWiki
Member



Mathematical functions

Function is a relation between inputs values and a set of outputs values such that each input is related to exactly one output.

DataMelt can be used to evaluate, plot, simplify functions. There are several approaches to work with parametric or non-parametric functions using DataMelt Java API:

  • User 2D/3D function plotter from DataMelt IDE (see the "Tools" menu)
  • Simplify and plot functions using ether Octave or Matlab syntax. Plotting functions in such approach is a bit too simplistic and there is no good integration with Java API and data structures
  • Use the full-blown DataMelt Java libraries to create parametric, non-parametric, special functions as they are provided by Java API. One can also integrate, differentiate and simplify such functions. In this case one can use the Java language, or any Java scripting languages (such as Jython or BeanShell). There is complete integration with data structures API and functions can be plotted in 2D/3D using several powerful interactive canvases.

In this section we will discuss the third option as it is the most flexible for integration with Java libraries. We will start with simple analytical functions with one variable ("x") and then will consider more complicated cases with several variables ("y", "z"). Finally, we will show how to make function of arbitrary complexity with any number of variables

Functions with one variable

Functions can be defined using either strings or programically. First we will discuss how to show functions using strings and then we will show other, more flexible methods.

We will start this tutorial with a simple example of how to create and display a function with one and only one variable ("x"). First, let us create an object representing a function:

[math]\displaystyle{ f(x)=2*exp(-x*x/50)+sin(pi*x)/x }[/math]

Such functions are represented by the jhplot.F1D jhplot.F1D Java class (1D means one-dimensional function).

from jhplot  import  *
f1 = F1D("2*exp(-x*x/50)+sin(pi*x)/x")

Here the input function is given by java.lang.String java.lang.String. One can evaluate this function at a certain point, say at 10, as:

print f1.eval(10)
0.270670566473

By default, all trigonometric functions are implemented using the fast calculations with e-14 accuracy. Note, you can also use exact calculations using </javadoc>java.lang.Math</javadoc> using additional parameter:

from jhplot  import  *
HParam.setExactMath()
f1 = F1D("2*exp(-x*x/50)+sin(pi*x)/x")

The method "HParam.setExactMath()" tells that full precision of Java Math package is not required. In this case, calculations of trigonometric functions is a factor 2-3 slower. Look at the class jhplot.HParam jhplot.HParam.

The code shown below benchmark calculations of F1D function containing trigonometric functions with and without the fast math enabled.

from jhplot import *
import time


HParam.setFastMath() # apply fast calculations for functions 
                     # comment out this code for full precision (java Math)
print "Fast calculations are set to ", HParam.isFastMath(), " working .."

a=F1D("acos(x)+asin(x)+atan(2*x)+cos(x)")
s=0;
start = time.time()
for i in range (10000000):
       s=s+a.eval(0.99-0.00000001*i)
end = time.time()
print "fast calculation result=",s, " time used (sec)",end - start



HParam.setFastMath(False)
print "Fast calculations are set to ", HParam.isFastMath(),  " working .."

s=0;
start = time.time()
for i in range (10000000):
       s=s+a.eval(0.99-0.00000001*i)
end = time.time()
print "fast calculation result=",s, " time used (sec)",end - start

To draw the function, one should specify limits during the function initialization (or later using the methods of the Java object jhplot.F1D jhplot.F1D. Then one should build a canvas and use the method "draw" to plot the function. In the example below, we define the limits for this function (-2,5) and then we plot it:

from jhplot  import *
f1 = F1D("2*exp(-x*x/50)+sin(pi*x)/x", 1.0, 10.0)
c1 = HPlot("Example of  function")
c1.visible()
c1.setAutoRange()
c1.draw(f1)

Here we drop the Python invitation ">>>" since we assume that the this code snippet is executed as a file. Insert these lines to the DataMelt editor and hit the key [F8]. You will see a pop-up window with the output as shown below:

DMelt example: Show a F1D function in the range 1-10

To learn about the method of DataMelt or any Jython or Java objects (like f1 and c1 in the above example):

  • Hit <Ctrl>+<Space> after the dot when using the Jython shell prompt (bottom part of the DataMelt IDE).
  • When using the DataMelt editor, use the key [F4] (after the dot) which brings up a table with all available Java methods.
  • Use Eclipse or NetBean IDE when developing the code using Java. They have their own help system.
  • Read the Java API of these classes.

Information about the syntax

Here is a short description about how to define functions. The following table shows the main mathematical operations:

()     parenthesis           
 +      plus                  
 -      minus                 
 *      times                 
 /      divide                
 ^      raise to a power      
 **     raise to a power

Functions:

String Definition
abs absolute value
log(x) Natural Logarithm
log10(x) Logarithm base 10
sqrt(x) <m>sqrt{x}</m>
exp(x) Exp
ceil nearest upper integer
floor nearest lower integer
cbrt cubic root
sin(x) Sine
cos(x) Cosine
tan(x) Tangent
asin(x) ArcSine
acos(x) ArcCosine
atan(x) ArcTangent
sinh(x) Hyperbolic Sine
cosh(x) Hyperbolic Cosine
tanh(x) Hyperbolic Tangent

One can also use some predefined constants, such as pi or Pi (the pi value, 3.14..). It is possible to use Scientific notation for numbers. The number is split into a significant / mantissa (y) and exponent (x) of the form 'yEx' which is evaluated as 'y * 10^x'. Let us give the examples:

from jhplot  import *
 
f1=F1D("1+x+(x^2)-(x^3)")   # correct answer -1
print f1.eval(2)            
f1=F1D("1+x+(x*x)-(x*x*x)") # correct answer -1
print f1.eval(2)            
f1=F1D("1+x+x^2-x^3")       # correct answer -1
print f1.eval(2)
f1=F1D("1+x+x**2-x**3")     # correct answer -1
print f1.eval(2)

Parameters

Functions can have parameters. If you have parameters, you should avoid "parsing" functions, and apply parsing only after finals substitution. This example shows how to work with such functions:

f1=F1D("(3*alpha)/alpha+10+x",0) # 0 means False
f1.simplify()                    # simplify the function
print f1.getName()            
f1.setPar("alpha",100)           # substitute alpha=100
f1.parse()                       # parse it for evaluation, plotting etc.
print f1.getName()
print f1.eval([1,2,3])

Please also consider the method f1.toString() that prints more information about jhplot.F1D jhplot.F1D.


Function initialization

Functions can be initialized with and without ranges.

When using this syntax:

from jhplot import *
f1=F1D("x^2")

function does not have limits used for plotting and evaluation. For such definition, plot ranges are entirely determined by the plotting canvases (i.e. by the setRange() methods.

When using this syntax:

from jhplot import *
f1=F1D("x^2",-10,10) # limit to [-10,10]

The function includes limits for numeric calculation and plotting. This means you can set auto-range for plotting canvases, since the function includes the X-range in the definition, thus it will be shown in the range [-10,10], as in this example.

Non-parametric functions

The most flexible way to draw functions is to use codding with objects and calling third-party Java libraries directly, instead of using strings with function definitions. This topic will be discussed in Section Non-parametric functions

Symbolic manipulations

Symbolic manipulation with functions will be discussed in more details in Section Introduction. Still, one should note that symbolic manipulations can also be done with the F1D functions using Java (or Jython) coding without using any other specialized language (as Matlab or Octave). One can simplify and expand F1D functions. Functions can also be rewritten in terms of elementary functions (log, sqrt, exp). Finally, one can perform some elementary substitutions before attempting to plot a F1D function.


In this example we simplify a few elementary functions:

from jhplot  import *

s="sqrt(17/12)-(1/6*sqrt(51))"
f1 = F1D(s)
f1.simplify()
print "Simplified "+s+" = ",f1.getName()

s= "(x^2-1)/(x-1)"
f1 = F1D(s)
f1.simplify()
print "Simplified "+s+" = ",f1.getName()

It is often useful to rewrite an expression in terms of elementary functions (log, exp, frac, sqrt, implicit roots), using the "elementary()" method, before simplifying it.

from jhplot  import *
s="cos(x)^2+sin(x)^2"
f1 = F1D(s)
f1.elementary()
print "In ementary functions "+s+" = ",f1.getName()
f1.simplify()
print "Simplified "+s+" = ",f1.getName()

Integration and differentiation

Functions can be numerically integrated. The program supports 5 methods of integration, which vary in evaluation time and accuracy. Below we will integrate the function "cos(x*x)+sin(1/x)*x^2)"

from jhplot import *
import time

f1=F1D('cos(x*x)+sin(1.0/x)*x^2')
methods=['trapezium','gauss4', 'gauss8',  'richardson', 'simpson'] 
for m in methods:
   start = time.clock()
   d=f1.integral(m,5000,10,100)
   t = time.clock()-start 
   print m,' =',d,' time (ms)=',t*1000

This code integrates the function between 10 and 100 using 5000 integration point (the large the number is, the more accurate the results are). The code also performs some benchmarking which gives you ideas about the accuracy of calculations:

trapezium  = 4949.64217622  time (ms)= 31.654937
gauss4     = 4949.64028115  time (ms)= 43.84014
gauss8     = 4949.64028111  time (ms)= 65.27855
richardson = 4949.64028393  time (ms)= 56.030938
simpson    = 4949.64014049  time (ms)= 17.634015

Displaying functions

Here is a more detailed example showing how to plot several functions using different colors:

from java.awt import Color
from jhplot  import *

# first function
f1 = F1D("2*exp(-x*x/50)+sin(pi*x)/x", 1.0, 10.0)
f1.setPenDash(4)

# second function
f2 = F1D("exp(-x*x/50)+pi*x", 1.0, 10.0)
f2.setColor(Color.green)
f2.setPenWidth(1)

# build a canvas with X and Y titles
c1 = HPlot("Canvas")
c1.setGTitle("2 functions", Color.red)
c1.setMarginLeft(90)
c1.setNameX("X")
c1.setNameY("Y")
c1.visible()
c1.setAutoRange()

c1.draw(f1)
c1.draw(f2)


Note that we have imported the Java class java.awt.Color java.awt.Color. The output of this script is shown here

DMelt example: Showing two functions with colors

You can also plot objects on different pads as shown in the Section Plotting canvases.

Exporting functions

As any Java object, F1D functions can be saved into files. Read Input and Output for detailed discussion. In addition, one can convert a function in MathML or Java formats, or just display it as a table with values after evaluation of a function.

>>> from jhplot import F1D
>>> f=F1D("10*sin(x)")
>>> print f.toMathML()  # convert to MathML
>>> print f.toJava()    # convert to Java code
>>> f.toTable()         # show a table with X-Y values (the function should be plotted or evaluated before)

2D Functions

Functions in 2 dimensions can be build analogously using the Java class jhplot.F2D jhplot.F2D. The example below shows how to construct and evaluate a function 2D:

>>> from jhplot  import  *
>>> f1 = F2D("2*exp(-x*y/50)+sin(pi*x)/y")
>>> print f1.eval(10,20)
0.0366312777775

The code prints the answer: 0.270670566473.

Such functions can be displayed as will be shown later. Also, look at the Sections 2D Plots and 3D Plots

3D Functions

Functions in 2 dimensions can be build by analogy using the Java jhplot.F3D jhplot.F3D class. The example below shows how to construct and evaluate a function 3D:

>>> from jhplot  import  *
>>> f1 = F3D("2*exp(-x*y/50)+sin(pi*x)/z")
>>> print f1.eval(10,20,30)
0.0366312777775

The code prints the answer: 0.0366312777775.

Converting functions to histograms

Histograms can be created from 1D and 2D functions. In the example above, we created a 2D function and then used it to generate 2D histogram. Functions can be converted to Histograms with arbitrary number of bins. This often can be used for show selected regions in different color.

Consider the example in which a function is used to create a histogram with fine bins. We use this histogram to highlight a range of our original function.

from java.awt import Color
from jhplot  import *

c1 = HPlot("Canvas",700,400,2, 1)
c1.setAutoRange()
c1.setGTitle("Converting to histogram", Color.red) #put title
c1.setNameX("X")
c1.setNameY("Y")
c1.visible(1)

c1.cd(1,1)
c1.setAutoRange()
f1 = F1D("2*exp(-x*x/50)+sin(pi*x)/x", -2.0, 5.0)
c1.draw(f1)

c1.cd(2,1)
c1.setAutoRange()
c1.setNameX("X")
c1.setNameY("Y")
h=f1.getH1D()
h.setColor(Color.red)
h.setFill(1)
h.setFillColor(Color.red)
h.setErrX(0)
h.setErrY(0)
c1.draw(h)

# export to some image (png,eps,pdf,jpeg...)
# c1.export(Editor.DocMaster

DMelt example: A function (F1D) converted to H1D histogram

Histograms can have arbitrary number of bins, but if the bin number is close to 500 (a typical evaluation step for functions), the difference between function and a histogram will be almost impossible to see.

Plotting in 3D

F2D functions can be shown using 3D canvases. They can be overlayed and shown on a single plot. Here is a small example of showing 2 functions in 3D:

from java.awt import Color
from jhplot import *

c1 = HPlot3D('Canvas',600,600)
c1.setNameX('X')
c1.setNameY('Y')
c1.visible()
f1=F2D('x^2+y^2', -2.0, 2.0, -2.0, 2.0)
f2=F2D('4*x*y', -2.0, 2.0, -2.0, 2.0)
c1.draw(f1,f2)

The output of this code is:

DMelt example: Plot two F2D functions in 3D on a single HPlot3D canvas

One can also overlay several different objects. Let us show an example of how to show functions and histograms in 3D:

DMelt example: H2D histogram and F2D function plotted on same 3D plot (HPlot3D)


The code used to generate this image is given below.

from jhplot  import HPlot3D,H2D,F2D
from java.util import Random

c1 = HPlot3D("Canvas",600,400)
c1.setGTitle("F2D and H2D objects")
c1.setTextBottom("Global X")
c1.setTextLeft("Global Y")


c1.setNameX("X axis")
c1.setNameY("Y axis")

c1.setColorMode(4)
c1.visible(1)

h1 = H2D("My 2D Test 1",30,-3.0, 3.0, 30, -3.0, 3.0)
f1 = F2D("8*(x*x+y*y)", -3.0, 3.0, -3.0, 5.0)
rand = Random()
for i in range(10000):
               h1.fill(0.4*rand.nextGaussian(),rand.nextGaussian())
c1.draw(h1,f1)

# export to some image (png,eps,pdf,jpeg...)
# c1.export(Editor.DocMasterName()+".png")


You can plot an arbitrary function in 3D using jhplot.HJavaView jhplot.HJavaView by defining function programmatically, and rebuilding vertexes for surface plots. Here is how you can do it:

# Show how to visualise a 2D function using Python code
# @author Sergei Chekanov 

from jv.geom import PgElementSet
from jhplot import HJavaView
from java.lang import Math
from math import *

geom = PgElementSet(3)
geom.setName("2D function")
# Set number of lines in x and y directions.
numXLines,numYLines = 50,50

def func(X,Y): # define some function in 2D using Python syntax 
    return 10*sin(sqrt(X*X+Y*Y))/sqrt(X*X+Y*Y);

# ranges in X and Y
xMin,xMax=-10.,10
yMin,yMax=-10.,10
zMin,zMax=0,10

geom.setNumVertices(numXLines*numYLines) # Allocate space for vertices

ind = 0 # Compute the vertices
x=xMin;
for i in xrange(numXLines):
      x = xMin + i*(xMax-xMin)/(numXLines-1)
      for j in xrange(numYLines):
         y = yMin + j*(yMax-yMin)/(numYLines-1)
         geom.setVertex(ind, x, y, func(x,y))
         ind=ind+1
# Compute connectivity of a rectangular mesh.
geom.makeQuadrConn(numXLines, numYLines)
geom.makeElementColorsFromXYZ() # add color and transparancy  
geom.showElementColors(True)
geom.setTransparency(0.5);
geom.showTransparency(True);

# show everythin 
c1= HJavaView()
c1.draw(geom)
view=c1.getView()
disp=view.getDisplay()
disp.showAxes(True)
disp.showBndBox(True);
c1.visible()

The output of this code is:

DMelt example: Show arbirary 2D surface function using JavaView

You can adjust all visual characteristics of this plot using JavaView GUI.

Multidimensional functions

So far we have learned how to build functions in 1D (F1D class), 2D (F2D class) and 3D (F3D class). In addition to these "fixed dimension" classes, there is a jhplot.FND jhplot.FND which can be used to build a function with arbitrary parameters and variables. To build such function is very easy and similar to the above styles, only now one should be more explicit about the names of variables. As example, below we make a function with 4 variables, called xx, yy, bb, kk:

from jhplot import *
f=FND("xx*sqrt(yy)*bb*+kk",'xx,yy,bb,kk')
ans=f.eval('xx=10,yy=1,bb=0,kk=20')
print ans

Such functions can be evaluated at any variable values and plotted.

Integration

Numeric integration of functions is discussed in section Integration. Symbolic integrations will be discussed in dedicated sections.

Let us give a small example showing how to integrate [math]\displaystyle{ cos(x)^3 }[/math] using a trapezium rule. We will integrate this function between 1 and 10 using 10k iterations.

from jhplot import F1D
f1=F1D('cos(x)^3')
print f1.integral(10000,1,10)

Differentiation

For differentiation, use the "diff" method. As argument, you must specify variable used for differentiation.

>>> from jhplot  import *
>>> fn=FND('x^4+x^2+x^4+y+z+h+10+20','x,y,z,h')
>>> fn.diff('x') # differentiate using x
>>> print 'Differentiate=',fn.toString()
Diffirentiate=  4.0*x^3.0+2.0*x+4.0*x^3.0

Simplification

F1D functions can be simplified using the method "simplify". The example below simplifies

>>> from jhplot  import *
>>> fn=FND('x^4+x^2+x^4+y+z+h+10+20','x,y,z,h')
>>> fn.simplify()
>>> print 'Simplify=',fn.toString()

The output is:

Simplify= 30.0+x^4.0+x^2.0+x^4.0+y+z+h

Minimization

Please read the section Minimization

Expression Builder

Multidimensional functions can also be built using the package jhplot.math.exp4j. In the example below we build a function using jhplot.math.exp4j.ExpressionBuilder jhplot.math.exp4j.ExpressionBuilder:

from java.awt import*
from jhplot  import *
from jhplot.math.exp4j import *

e=ExpressionBuilder("3 * sin(x) - 2 / (x - 2)")
e.variables("x")

ff=F1D("Test", e.build() ,10,100)
print ff.eval(10) # check at x=10

c1 = HPlot()
c1.visible(1)
c1.setAutoRange()
c1.draw(ff)

You can build an arbitrary function using Python and plot it, mixing Java and Jython libraries in one code. You can put an arbitrary complex expressions and "if" statements. In the example below we make a custom function using a Jython class and build an "ExpressionBuilder" object. then we evaluate this function for array of "X" values and plot it.

from java.awt import* 
from jhplot  import * 
from jhplot.math.exp4j import *
from jhplot.math.exp4j.function import * 


class cf2(Function):
  def  apply(self,val):
    if (val[0]<100):  
                return val[0]*val[0] 
    if (val[0]>=100): 
                return  val[0] 

e=ExpressionBuilder("cf2(x)")
func=cf2("cf2",1);
e.function(func)
e.variables("x")
fff=e.build()
#fff.setVariable("x", 10);
#print fff.evaluate()

c1 = HPlot()
c1.visible(1)
c1.setAutoRange()
f=F1D(e.build(), 10,1000)
c1.draw(f)

Now we can show how to create a custom function in 3D. Again, we will use a python code to do this.

DMelt example: Showing a custom function using exp4j in 3D

from java.awt import* 
from jhplot  import * 
from jhplot.math.exp4j import *
from jhplot.math.exp4j.function import * 


class cf2d(Function):
  def  apply(self,val):
    if (val[0]<50 or val[1]<50):  
                return val[1]*val[0] 
    if (val[0]>=50 and val[1]>=50): 
                return  val[0]*val[0]+val[1]*val[1] 

e=ExpressionBuilder("cf2d(x,y)")
func=cf2d("cf2d",2);
e.function(func)
e.variables("x")
e.variables("y")
fff=e.build()
#fff.setVariable("x", 10);
#print fff.evaluate()

c1 = HPlot3D()
c1.visible(1)
c1.setAutoRange()
f=F2D(e.build(), 0,200, 0,200)
c1.draw(f)

Fast Math

The implementation of trigonometric functions is a 3-4 times slower performance compare to C/C++. This is especially true for atan2(), acos(), asin() functions. By default, DMelt uses fast calculations for all F1D, F2D etc. functions. The table below shows the performance of FastMath compared to the standard Java implementation blog. The times are in mean nanoseconds.

<html>

<thead> </thead> <tbody> </tbody>
Function Math Mean FastMath Mean Times Faster
acos 58 16 3.6
asin 57 15 3.8
atan 94 15 6.2
atan2 145 23 6.3
cbrt 112 18 6.2
ceil 10 11 0.9
cos 74 13 5.7
cosh 141 28 5.0
exp 73 16 4.6
expm1 113 17 6.6
floor 11 12 0.9
hypot 435 23 18.9
log 31 16 1.9
log10 27 13 2.1
log1p 97 15 6.5
nextAfter 11 10 1.1
nextUp 10 8 1.3
pow 101 38 2.7
round 11 10 1.1
sin 72 16 4.5
sinh 126 23 5.5
sqrt 7 15 0.5
tan 52 14 3.7
tanh 167 26 6.4
toDegrees 7 8 0.9
toRadians 7 7 1.0

</html>

DMelt also include less precise methods for fast evaluation basic functions using jhplot.math.MathUtilsd jhplot.math.MathUtilsd. They are faster than the functions shown in above table, but their accuracy is roughly 1E-05.

Using special functions

You can integrate special functions from third-party Java libraries into your code and show them in 2D and 3D.

Parametric functions

Please read the section Parametric functions.