Attachment "symdiff.test" to
ticket [1328920fff]
added by
kennykb
2005-10-18 02:28:30.
# symdiff.test --
# Test cases for the 'symdiff' package
#
# This file contains a collection of tests for one or more of the Tcllib
# procedures. Sourcing this file into Tcl runs the tests and
# generates output for errors. No output means no errors were found.
#
# Copyright (c) 2005 by Kevin B. Kenny
# All rights reserved.
#
# See the file "license.terms" for information on usage and redistribution
# of this file, and for a DISCLAIMER OF ALL WARRANTIES.
#
# RCS: @(#) $Id: symdiff.test not yet committed $
if {[lsearch [namespace children] ::tcltest] == -1} {
package require tcltest 2.1
namespace import ::tcltest::*
} else {
# Ensure that 2.1 or higher present.
if {![package vsatisfies [package present tcltest] 2.1]} {
puts "Aborting tests for math::statistics."
puts "Requiring tcltest 2.1, have [package present tcltest]"
return
}
}
source [file join [pwd] [file dirname [info script]] symdiff.tcl]
package require math::symdiff 1.0
namespace eval ::math::symdiff::test {
namespace import ::tcltest::test
namespace import ::tcltest::cleanupTests
namespace import ::math::symdiff::*
set prec $::tcl_precision
set ::tcl_precision 17
test symdiff-1.1 {derivative of a constant} {
symdiff {1.0} a
} 0.0
test symdiff-2.1 {derivative of a variable} {
symdiff {$a} a
} 1.0
test symdiff-2.2 {derivative of a variable} {
symdiff {$b} a
} 0.0
test symdiff-3.1 {derivative of a sum, easy cases} {
symdiff {1.0 + 1.0} a
} 0.0
test symdiff-3.2 {derivative of a sum, easy cases} {
symdiff {1.0 + $a} a
} 1.0
test symdiff-3.3 {derivative of a sum, easy cases} {
symdiff {$a + 1.0} a
} 1.0
test symdiff-3.4 {derivative of a sum, easy cases} {
symdiff {$a + $a} a
} 2.0
test symdiff-3.5 {derivative of a sum, easy cases} {
symdiff {$a + $b} a
} 1.0
test symdiff-3.6 {derivative of a sum, easy cases} {
symdiff {$a + $a + $a} a
} 3.0
test symdiff-4.1 {derivative of a difference, easy cases} {
symdiff {1.0 - 1.0} a
} 0.0
test symdiff-4.2 {derivative of a difference, easy cases} {
symdiff {1.0 - $a} a
} -1.0
test symdiff-4.3 {derivative of a difference, easy cases} {
symdiff {$a - 1.0} a
} 1.0
test symdiff-4.4 {derivative of a difference, easy cases} {
symdiff {$a - $a} a
} 0.0
test symdiff-4.5 {derivative of a difference, easy cases} {
symdiff {$a + $b} a
} 1.0
test symdiff-4.6 {derivative of a difference, easy cases} {
symdiff {$a + $a - $a} a
} 1.0
test symdiff-5.1 {derivative of a product, easy cases} {
symdiff {1.0 * 1.0} a
} 0.0
test symdiff-5.2 {derivative of a product, easy cases} {
symdiff {3.0 * $a} a
} 3.0
test symdiff-5.3 {derivative of a product, easy cases} {
symdiff {$a * 3.0} a
} 3.0
test symdiff-5.4 {derivative of a product, easy cases} {
symdiff {$a * $a} a
} {($a + $a)}
test symdiff-5.5 {derivative of a product, easy cases} {
symdiff {$a * $b} a
} {$b}
test symdiff-5.6 {derivative of a product, easy cases} {
symdiff {($a + $b) * ($a + $b)} a
} {(($a + $b) + ($a + $b))}
test symdiff-5.7 {derivative of a linear function} {
symdiff {$a*$x + $b} x
} {$a}
test symdiff-6.1 {derivative of a sum} {
symdiff {($a*$x+$b)+($c*$x+$d)} x
} {($a + $c)}
test symdiff-7.1 {derivative of a difference} {
symdiff {($a*$x+$b)-($c*$x+$d)} x
} {($a - $c)}
test symdiff-8.1 {derivative of a product} {
symdiff {($a*$x+$b)*($c*$x+$d)} x
} {(($c * (($a * $x) + $b)) + ($a * (($c * $x) + $d)))}
test symdiff-9.1 {derivative of a quotient} {
symdiff {$x/1.0} x
} 1.0
test symdiff-9.2 {derivative of a quotient} {
symdiff {$x/-1.0} x
} -1.0
test symdiff-9.3 {derivative of a quotient} {
symdiff {1.0/$x} x
} {-(((1.0 / $x) / $x))}
test symdiff-9.4 {derivative of a quotient} {
symdiff {($a*$x+$b)/($c*$x+$d)} x
} {(($a - (($c * (($a * $x) + $b)) / (($c * $x) + $d))) / (($c * $x) + $d))}
test symdiff-10.1 {derivative of an exponent} {
symdiff {pow($a*$x+$b,3.5)} x
} {($a * (3.5 * pow((($a * $x) + $b), 2.5)))}
test symdiff-10.2 {derivative of an exponent, slightly harder case} {
symdiff {pow(10.0,$x)} x
} {(pow(10.0, $x) * 2.3025850929940459)}
test symdiff-10.3 {derivative of an exponent, awkward case} {
symdiff {pow($a*$x+$b,$c*$x+$d)} x
} {(pow((($a * $x) + $b), (($c * $x) + $d)) * ((($a * (($c * $x) + $d)) / (($a * $x) + $b)) + ($c * log((($a * $x) + $b)))))}
test symdiff-11.1 {derivative of a unary negation} {
symdiff {-($a*$x + $b)} x
} {-($a)}
test symdiff-11.2 {derivative of a unary plus} {
symdiff {+($a*$x + $b)} x
} {$a}
test symdiff-12.1 {derivative of acos} {
symdiff {acos($x)} x
} {(-1.0 / sqrt((1.0 - ($x * $x))))}
test symdiff-12.2 {derivative of acos} {
symdiff {acos($a*$x+$b)} x
} {-(($a / sqrt((1.0 - ((($a * $x) + $b) * (($a * $x) + $b))))))}
test symdiff-13.1 {derivative of acos} {
symdiff {asin($x)} x
} {(1.0 / sqrt((1.0 - ($x * $x))))}
test symdiff-13.2 {derivative of asin} {
symdiff {asin($a*$x+$b)} x
} {($a / sqrt((1.0 - ((($a * $x) + $b) * (($a * $x) + $b)))))}
test symdiff-14.1 {derivative of atan} {
symdiff {atan($x)} x
} {(1.0 / (1.0 + ($x * $x)))}
test symdiff-14.2 {derivative of atan} {
symdiff {atan($a*$x+$b)} x
} {($a / (1.0 + ((($a * $x) + $b) * (($a * $x) + $b))))}
test symdiff-15.1 {derivative of atan2} {
symdiff {atan2($x,1.0)} x
} {(1.0 / (($x * $x) + 1.0))}
test symdiff-15.2 {derivative of atan2} {
symdiff {atan2(1.0,$x)} x
} {(-1.0 / (1.0 + ($x * $x)))}
test symdiff-15.3 {derivative of atan2} {
symdiff {atan2($x,$y)} x
} {($y / (($x * $x) + ($y * $y)))}
test symdiff-15.4 {derivative of atan2} {
symdiff {atan2($y,$x)} x
} {-(($y / (($y * $y) + ($x * $x))))}
test symdiff-15.5 {derivative of atan2} {
symdiff {atan2($a*$x+$b,$c*$x+$d)} x
} {((($a * (($c * $x) + $d)) - ((($a * $x) + $b) * $c)) / (((($a * $x) + $b) * (($a * $x) + $b)) + ((($c * $x) + $d) * (($c * $x) + $d))))}
test symdiff-16.1 {derivative of cos} {
symdiff {cos($x)} x
} {-(sin($x))}
test symdiff-16.2 {derivative of cos} {
symdiff {cos($a*$x + $b)} x
} {-(($a * sin((($a * $x) + $b))))}
test symdiff-17.1 {derivative of cosh} {
symdiff {cosh($x)} x
} {sinh($x)}
test symdiff-17.2 {derivative of cosh} {
symdiff {cosh($a*$x + $b)} x
} {($a * sinh((($a * $x) + $b)))}
test symdiff-18.1 {derivative of exp} {
symdiff {exp($x)} x
} {exp($x)}
test symdiff-18.2 {derivative of exp} {
symdiff {exp($a*$x+$b)} x
} {($a * exp((($a * $x) + $b)))}
test symdiff-19.1 {derivative of hypot} {
symdiff {hypot(0.0,$a)} a
} {($a / hypot(0.0, $a))}
test symdiff-19.2 {derivative of hypot} {
symdiff {hypot($b,$a)} a
} {($a / hypot($b, $a))}
test symdiff-19.3 {derivative of hypot} {
symdiff {hypot($a*$x+$b,$c*$x+$d)} x
} {((($a * (($a * $x) + $b)) + ($c * (($c * $x) + $d))) / hypot((($a * $x) + $b), (($c * $x) + $d)))}
test symdiff-20.1 {derivative of log} {
symdiff {log($x)} x
} {(1.0 / $x)}
test symdiff-20.2 {derivative of log} {
symdiff {log($a*$x+$b)} x
} {($a / (($a * $x) + $b))}
test symdiff-21.1 {derivative of log10} {
symdiff {log10($x)} x
} {(1.0 / (2.3025850929940459 * $x))}
test symdiff-21.2 {derivative of log10} {
symdiff {log10($a * $x + $b)} x
} {($a / (2.3025850929940459 * (($a * $x) + $b)))}
test symdiff-22.1 {derivative of sin} {
symdiff {sin($x)} x
} {cos($x)}
test symdiff-22.2 {derivative of sin} {
symdiff {sin($a*$x+$b)} x
} {($a * cos((($a * $x) + $b)))}
test symdiff-22.1 {derivative of sinh} {
symdiff {sinh($x)} x
} {cosh($x)}
test symdiff-22.2 {derivative of sinh} {
symdiff {sinh($a*$x+$b)} x
} {($a * cosh((($a * $x) + $b)))}
test symdiff-23.1 {derivative of sqrt} {
symdiff {sqrt($x)} x
} {(1.0 / (2.0 * sqrt($x)))}
test symdiff-23.2 {derivative of sqrt} {
symdiff {sqrt($a*$x+$b)} x
} {($a / (2.0 * sqrt((($a * $x) + $b))))}
test symdiff-24.1 {derivative of tan} {
symdiff {tan($x)} x
} {(1.0 / (cos($x) * cos($x)))}
test symdiff-24.2 {derivative of tan} {
symdiff {tan($a*$x+$b)} x
} {($a / (cos((($a * $x) + $b)) * cos((($a * $x) + $b))))}
test symdiff-24.1 {derivative of tanh} {
symdiff {tanh($x)} x
} {(1.0 / (cosh($x) * cosh($x)))}
test symdiff-24.2 {derivative of tanh} {
symdiff {tanh($a*$x+$b)} x
} {($a / (cosh((($a * $x) + $b)) * cosh((($a * $x) + $b))))}
test symdiff-25.1 {error handling} {
list [catch {symdiff {[foo $x]} x} result] $result
} {1 {symdiff doesn't handle command expressions}}
test symdiff-25.2 {error handling} {
list [catch {symdiff {$x(1)} x} result] $result
} {1 {symdiff doesn't handle arrays}}
test symdiff-25.3 {error handling} {
list [catch {symdiff {$a & $b} a} result] $result
} {1 {symdiff can't differentiate the "&" operator}}
test symdiff-25.4 {error handling} {
list [catch {symdiff {int($a)} a} result] $result
} {1 {symdiff can't differentiate the "int" function}}
test symdiff-25.5 {error handling} {
list [catch {symdiff {$a ? $b : $c} a} result] $result
} {1 {symdiff can't differentiate the ternary ?: operator}}
test symdiff-26.1 {unary minus optimization} {
symdiff {$a * $x + -$b * $x} x
} {($a - $b)}
test symdiff-26.2 {unary minus optimization} {
symdiff {-$a * $x - $b * $x} x
} {-(($a + $b))}
test symdiff-26.3 {unary minus optimization} {
symdiff {$a * $x - -$b * $x} x
} {($a + $b)}
test symdiff-26.4 {unary minus optimization} {
symdiff {-$a * $x * $b} x
} {-(($a * $b))}
test symdiff-26.5 {unary minus optimization} {
symdiff {$a * $x * -$b} x
} {-(($a * $b))}
test symdiff-26.6 {unary minus optimization} {
symdiff {---($a*$x+$b)} x
} {-($a)}
test symdiff-26.7 {unary minus optimization} {
symdiff {-$x * $x} x
} {-(($x + $x))}
test symdiff-27.1 {power optimizations} {
symdiff {pow($x,1)} x
} 1.0
test symdiff-27.2 {power optimizations} {
symdiff {pow($x,2.0)} x
} {(2.0 * $x)}
test symdiff-28.1 {quotient optimization} {
symdiff {($x * $x) / 1.0} x
} {($x + $x)}
test symdiff-28.2 {quotient optimization} {
symdiff {($x * $x) / -1.0} x
} {-(($x + $x))}
test symdiff-28.3 {quotient optimization - error case} {
list [catch {symdiff {($x * $x) / 0.0} x} result] $result
} {1 {requested expression will result in division by zero at run time}}
test symdiff-29.1 {product optimization} {
symdiff {(2. * $x) * 3.0} x
} 6.0
test symdiff-29.2 {product optimization} {
symdiff {($a * $x) * -1.0} x
} {-($a)}
test symdiff-30.0 {illustration of Newton's method - find a root of sin(x)-0.5 near 0.5} {
proc root {expr var guess} {
upvar 1 $var v
set deriv [symdiff $expr $var]
set v $guess
set updateExpr [list expr "\$$var - ($expr) / ($deriv)"]
for { set i 0 } { $i < 4 } { incr i } {
set v [uplevel 1 $updateExpr]
}
return $v
}
set r [root {sin($x)-0.5} x 0.5]
expr {sin($r)}
} 0.5
# End of test cases
set ::tcl_precision $prec
cleanupTests
}
namespace delete ::math::symdiff::test
# Local Variables:
# mode: tcl
# End:
