File: test-convex.tcl

package info (click to toggle)
tklib 0.6%2B20190108-1
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 15,008 kB
  • sloc: tcl: 75,757; sh: 5,789; ansic: 792; pascal: 359; makefile: 70; sed: 53; exp: 21
file content (84 lines) | stat: -rwxr-xr-x 2,003 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
#!/usr/bin/env tclsh
## -*- tcl -*-

package require Plotchart

proc convex {u v} {

    #expr {0.1*($u-$v)**2 - ($u + $v)/ sqrt(2.0) + 2.5}
    #expr {3.0 - 0.5 * ($u**2 + $v**2 - 2.0 * $u * $v) - ($u+$v)/sqrt(2.0)}

    expr {10.0 - $u**2/1.1 - $v**2/1.2}
}

pack [canvas .c -width 500 -height 500]
set p [::Plotchart::createXYPlot .c {-5 5 2.5} {-5 5 2.5}]

#::Plotchart::colorMap jet

#set xv {-2 -1.75 -1.5 -1.25 -1.0 -0.75 -0.5 -0.25 0.0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 2.0 2.5 3.0 3.5 4.0 4.5}
#set yv {-2 -1.75 -1.5 -1.25 -1.0 -0.75 -0.5 -0.25 0.0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 2.0 2.5 3.0 3.5 4.0 4.5}

for {set i 0} { $i < 20 } { incr i } {
    lappend xv [expr {($i-10)*0.5}]
    lappend yv [expr {($i-10)*0.5}]
}

catch { console show }
set zv {}
foreach y $yv {
    set row {}
    foreach x $xv {
        lappend row [convex $x $y]
    }
    lappend zv $row
}

foreach y {-1.0 0.0 1.0} {
    foreach x {-1.0 0.0 1.0} {
        $p labeldot $x $y [format %.2f [convex $x $y]]
    }
}
$p labeldot  4 0 [format %.2f [convex 4 0]] s
$p labeldot  0 4 [format %.2f [convex 0 4]]

$p plot axisx -5  0
$p plot axisx  5  0
$p plot axisy  0 -5
$p plot axisy  0  5

#$p contourlinesfunctionvalues $xv $yv $zv {-2.0 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 3.0 4.0}
#$p contourlinesfunctionvalues $xv $yv $zv {-2.0 0.0 2.0}

$p contourlinesfunctionvalues $xv $yv $zv {-2.0 0.0 2.0}
$p legendisolines $zv {-2.0 0.0 2.0}


toplevel .t
pack [canvas .t.c -width 500 -height 500]
set p2 [::Plotchart::createXYPlot .t.c {-5 5 2.5} {-5 5 2.5}]

set zv {}
set xcrd {}
set ycrd {}
foreach y $yv {
    lappend xcrd $xv
    set row {}
    set yc  {}
    foreach x $xv {
        lappend yc  $y
        lappend row [convex $x $y]
    }
    lappend zv   $row
    lappend ycrd $yc
}


$p2 contourfill  $xcrd $ycrd $zv {-2.0 0.0 2.0}
$p2 legendshades $zv {-2.0 0.0 2.0}

#
# contourfill requires full matrices for x and y
#$p contourfill $xv $yv $zv {-2.0 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 3.0 4.0}
#$p contourfill $xv $yv $zv 0.0