| line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
|
1
|
|
|
|
|
|
|
=encoding utf8 |
|
2
|
|
|
|
|
|
|
|
|
3
|
|
|
|
|
|
|
=head1 NAME |
|
4
|
|
|
|
|
|
|
|
|
5
|
|
|
|
|
|
|
Math::Symbolic - Symbolic calculations |
|
6
|
|
|
|
|
|
|
|
|
7
|
|
|
|
|
|
|
=head1 SYNOPSIS |
|
8
|
|
|
|
|
|
|
|
|
9
|
|
|
|
|
|
|
use Math::Symbolic; |
|
10
|
|
|
|
|
|
|
|
|
11
|
|
|
|
|
|
|
my $tree = Math::Symbolic->parse_from_string('1/2 * m * v^2'); |
|
12
|
|
|
|
|
|
|
# Now do symbolic calculations with $tree. |
|
13
|
|
|
|
|
|
|
# ... like deriving it... |
|
14
|
|
|
|
|
|
|
|
|
15
|
|
|
|
|
|
|
my ($sub) = Math::Symbolic::Compiler->compile_to_sub($tree); |
|
16
|
|
|
|
|
|
|
|
|
17
|
|
|
|
|
|
|
my $kinetic_energy = $sub->($mass, $velocity); |
|
18
|
|
|
|
|
|
|
|
|
19
|
|
|
|
|
|
|
=head1 DESCRIPTION |
|
20
|
|
|
|
|
|
|
|
|
21
|
|
|
|
|
|
|
Math::Symbolic is intended to offer symbolic calculation capabilities |
|
22
|
|
|
|
|
|
|
to the Perl programmer without using external (and commercial) libraries |
|
23
|
|
|
|
|
|
|
and/or applications. |
|
24
|
|
|
|
|
|
|
|
|
25
|
|
|
|
|
|
|
Unless, however, some interested and knowledgable developers turn up to |
|
26
|
|
|
|
|
|
|
participate in the development, the library will be severely limited by |
|
27
|
|
|
|
|
|
|
my experience in the area. Symbolic calculations are an active field of |
|
28
|
|
|
|
|
|
|
research in CS. |
|
29
|
|
|
|
|
|
|
|
|
30
|
|
|
|
|
|
|
There are several ways to construct Math::Symbolic trees. There are no |
|
31
|
|
|
|
|
|
|
actual Math::Symbolic objects, but rather trees of objects of subclasses of |
|
32
|
|
|
|
|
|
|
Math::Symbolic. The most general but unfortunately also the least intuitive |
|
33
|
|
|
|
|
|
|
way of constructing trees is to use the constructors of |
|
34
|
|
|
|
|
|
|
the Math::Symbolic::Operator, Math::Symbolic::Variable, and |
|
35
|
|
|
|
|
|
|
Math::Symbolic::Constant classes to create (nested) objects of the |
|
36
|
|
|
|
|
|
|
corresponding types. |
|
37
|
|
|
|
|
|
|
|
|
38
|
|
|
|
|
|
|
Furthermore, you may use the overloaded interface to apply the standard |
|
39
|
|
|
|
|
|
|
Perl operators (and functions, see L) to existing |
|
40
|
|
|
|
|
|
|
Math::Symbolic trees and standard Perl expressions. |
|
41
|
|
|
|
|
|
|
|
|
42
|
|
|
|
|
|
|
Possibly the most convenient way of constructing Math::Symbolic trees is |
|
43
|
|
|
|
|
|
|
using the builtin parser to generate trees from expressions such as C<2 * x^5>. |
|
44
|
|
|
|
|
|
|
You may use the Cparse_from_string()> class method for this. |
|
45
|
|
|
|
|
|
|
|
|
46
|
|
|
|
|
|
|
Of course, you may combine the overloaded interface with the parser to |
|
47
|
|
|
|
|
|
|
generate trees with Perl code such as C<$term * 5 * 'sin(omega*t+phi)'> which |
|
48
|
|
|
|
|
|
|
will create a tree of the existing tree $term times 5 times the sine of |
|
49
|
|
|
|
|
|
|
the vars omega times t plus phi. |
|
50
|
|
|
|
|
|
|
|
|
51
|
|
|
|
|
|
|
There are several modules in the distribution that contain subroutines |
|
52
|
|
|
|
|
|
|
related to calculus. These are not loaded by Math::Symbolic by default. |
|
53
|
|
|
|
|
|
|
Furthermore, there are several extensions to Math::Symbolic available |
|
54
|
|
|
|
|
|
|
from CPAN as separate distributions. Please refer to L |
|
55
|
|
|
|
|
|
|
for an incomplete list of these. |
|
56
|
|
|
|
|
|
|
|
|
57
|
|
|
|
|
|
|
For example, L come with C and |
|
58
|
|
|
|
|
|
|
contains routines to compute Taylor Polynomials and the associated errors. |
|
59
|
|
|
|
|
|
|
|
|
60
|
|
|
|
|
|
|
Routines related to vector calculus such as grad, div, rot, and Jacobi- and |
|
61
|
|
|
|
|
|
|
Hesse matrices are available through the L |
|
62
|
|
|
|
|
|
|
module. This module is also able to compute Taylor Polynomials of |
|
63
|
|
|
|
|
|
|
functions of two variables, directional derivatives, total differentials, |
|
64
|
|
|
|
|
|
|
and Wronskian Determinants. |
|
65
|
|
|
|
|
|
|
|
|
66
|
|
|
|
|
|
|
Some basic support for linear algebra can be found in |
|
67
|
|
|
|
|
|
|
L. This includes a routine to compute |
|
68
|
|
|
|
|
|
|
the determinant of a matrix of C trees. |
|
69
|
|
|
|
|
|
|
|
|
70
|
|
|
|
|
|
|
=head2 EXPORT |
|
71
|
|
|
|
|
|
|
|
|
72
|
|
|
|
|
|
|
None by default, but you may choose to have the following constants |
|
73
|
|
|
|
|
|
|
exported to your namespace using the standard Exporter semantics. |
|
74
|
|
|
|
|
|
|
There are two export tags: :all and :constants. :all will export |
|
75
|
|
|
|
|
|
|
all constants and the parse_from_string subroutine. |
|
76
|
|
|
|
|
|
|
|
|
77
|
|
|
|
|
|
|
Constants for transcendetal numbers: |
|
78
|
|
|
|
|
|
|
EULER (2.7182...) |
|
79
|
|
|
|
|
|
|
PI (3.14159...) |
|
80
|
|
|
|
|
|
|
|
|
81
|
|
|
|
|
|
|
Constants representing operator types: (First letter indicates arity) |
|
82
|
|
|
|
|
|
|
(These evaluate to the same numbers that are returned by the type() |
|
83
|
|
|
|
|
|
|
method of Math::Symbolic::Operator objects.) |
|
84
|
|
|
|
|
|
|
B_SUM |
|
85
|
|
|
|
|
|
|
B_DIFFERENCE |
|
86
|
|
|
|
|
|
|
B_PRODUCT |
|
87
|
|
|
|
|
|
|
B_DIVISION |
|
88
|
|
|
|
|
|
|
B_LOG |
|
89
|
|
|
|
|
|
|
B_EXP |
|
90
|
|
|
|
|
|
|
U_MINUS |
|
91
|
|
|
|
|
|
|
U_P_DERIVATIVE (partial derivative) |
|
92
|
|
|
|
|
|
|
U_T_DERIVATIVE (total derivative) |
|
93
|
|
|
|
|
|
|
U_SINE |
|
94
|
|
|
|
|
|
|
U_COSINE |
|
95
|
|
|
|
|
|
|
U_TANGENT |
|
96
|
|
|
|
|
|
|
U_COTANGENT |
|
97
|
|
|
|
|
|
|
U_ARCSINE |
|
98
|
|
|
|
|
|
|
U_ARCCOSINE |
|
99
|
|
|
|
|
|
|
U_ARCTANGENT |
|
100
|
|
|
|
|
|
|
U_ARCCOTANGENT |
|
101
|
|
|
|
|
|
|
U_SINE_H |
|
102
|
|
|
|
|
|
|
U_COSINE_H |
|
103
|
|
|
|
|
|
|
U_AREASINE_H |
|
104
|
|
|
|
|
|
|
U_AREACOSINE_H |
|
105
|
|
|
|
|
|
|
B_ARCTANGENT_TWO |
|
106
|
|
|
|
|
|
|
|
|
107
|
|
|
|
|
|
|
Constants representing Math::Symbolic term types: |
|
108
|
|
|
|
|
|
|
(These evaluate to the same numbers that are returned by the term_type() |
|
109
|
|
|
|
|
|
|
methods.) |
|
110
|
|
|
|
|
|
|
T_OPERATOR |
|
111
|
|
|
|
|
|
|
T_CONSTANT |
|
112
|
|
|
|
|
|
|
T_VARIABLE |
|
113
|
|
|
|
|
|
|
|
|
114
|
|
|
|
|
|
|
Subroutines: |
|
115
|
|
|
|
|
|
|
parse_from_string (returns Math::Symbolic tree) |
|
116
|
|
|
|
|
|
|
|
|
117
|
|
|
|
|
|
|
=cut |
|
118
|
|
|
|
|
|
|
|
|
119
|
|
|
|
|
|
|
package Math::Symbolic; |
|
120
|
|
|
|
|
|
|
|
|
121
|
23
|
|
|
23
|
|
738683
|
use 5.006; |
|
|
23
|
|
|
|
|
93
|
|
|
|
23
|
|
|
|
|
989
|
|
|
122
|
23
|
|
|
23
|
|
164
|
use strict; |
|
|
23
|
|
|
|
|
47
|
|
|
|
23
|
|
|
|
|
955
|
|
|
123
|
23
|
|
|
23
|
|
132
|
use warnings; |
|
|
23
|
|
|
|
|
76
|
|
|
|
23
|
|
|
|
|
844
|
|
|
124
|
|
|
|
|
|
|
|
|
125
|
23
|
|
|
23
|
|
148
|
use Carp; |
|
|
23
|
|
|
|
|
44
|
|
|
|
23
|
|
|
|
|
2596
|
|
|
126
|
|
|
|
|
|
|
|
|
127
|
23
|
|
|
23
|
|
13728
|
use Math::Symbolic::ExportConstants qw/:all/; |
|
|
23
|
|
|
|
|
60
|
|
|
|
23
|
|
|
|
|
6855
|
|
|
128
|
23
|
|
|
23
|
|
30040
|
use Math::Symbolic::AuxFunctions; |
|
|
23
|
|
|
|
|
80
|
|
|
|
23
|
|
|
|
|
884
|
|
|
129
|
|
|
|
|
|
|
|
|
130
|
23
|
|
|
23
|
|
19621
|
use Math::Symbolic::Base; |
|
|
23
|
|
|
|
|
81
|
|
|
|
23
|
|
|
|
|
794
|
|
|
131
|
23
|
|
|
23
|
|
42663
|
use Math::Symbolic::Operator; |
|
|
23
|
|
|
|
|
89
|
|
|
|
23
|
|
|
|
|
908
|
|
|
132
|
23
|
|
|
23
|
|
17838
|
use Math::Symbolic::Variable; |
|
|
23
|
|
|
|
|
66
|
|
|
|
23
|
|
|
|
|
705
|
|
|
133
|
23
|
|
|
23
|
|
13866
|
use Math::Symbolic::Constant; |
|
|
23
|
|
|
|
|
66
|
|
|
|
23
|
|
|
|
|
674
|
|
|
134
|
|
|
|
|
|
|
|
|
135
|
23
|
|
|
23
|
|
143
|
use Math::Symbolic::Derivative; |
|
|
23
|
|
|
|
|
44
|
|
|
|
23
|
|
|
|
|
1000
|
|
|
136
|
|
|
|
|
|
|
|
|
137
|
23
|
|
|
23
|
|
16184
|
use Math::Symbolic::Parser; |
|
|
23
|
|
|
|
|
68
|
|
|
|
23
|
|
|
|
|
765
|
|
|
138
|
23
|
|
|
23
|
|
16450
|
use Math::Symbolic::Compiler; |
|
|
23
|
|
|
|
|
75
|
|
|
|
23
|
|
|
|
|
1041
|
|
|
139
|
|
|
|
|
|
|
|
|
140
|
23
|
|
|
23
|
|
13849
|
use Math::Symbolic::Custom; |
|
|
23
|
|
|
|
|
87
|
|
|
|
23
|
|
|
|
|
7431
|
|
|
141
|
|
|
|
|
|
|
|
|
142
|
|
|
|
|
|
|
require Exporter; |
|
143
|
|
|
|
|
|
|
|
|
144
|
|
|
|
|
|
|
our @ISA = qw(Exporter); |
|
145
|
|
|
|
|
|
|
|
|
146
|
|
|
|
|
|
|
our %EXPORT_TAGS = ( |
|
147
|
|
|
|
|
|
|
all => [ |
|
148
|
|
|
|
|
|
|
@{ $Math::Symbolic::ExportConstants::EXPORT_TAGS{all} }, |
|
149
|
|
|
|
|
|
|
qw{&parse_from_string}, |
|
150
|
|
|
|
|
|
|
], |
|
151
|
|
|
|
|
|
|
constants => [ @{ $Math::Symbolic::ExportConstants::EXPORT_TAGS{all} }, ], |
|
152
|
|
|
|
|
|
|
); |
|
153
|
|
|
|
|
|
|
our @EXPORT_OK = ( @{ $EXPORT_TAGS{'all'} } ); |
|
154
|
|
|
|
|
|
|
our @EXPORT = qw(); |
|
155
|
|
|
|
|
|
|
|
|
156
|
|
|
|
|
|
|
our $VERSION = '0.612'; |
|
157
|
|
|
|
|
|
|
|
|
158
|
|
|
|
|
|
|
=head1 CLASS DATA |
|
159
|
|
|
|
|
|
|
|
|
160
|
|
|
|
|
|
|
The package variable $Parser will contain a Parse::RecDescent |
|
161
|
|
|
|
|
|
|
object that is used to parse strings at runtime. |
|
162
|
|
|
|
|
|
|
|
|
163
|
|
|
|
|
|
|
=cut |
|
164
|
|
|
|
|
|
|
|
|
165
|
|
|
|
|
|
|
our $Parser = Math::Symbolic::Parser->new(); |
|
166
|
|
|
|
|
|
|
|
|
167
|
|
|
|
|
|
|
=head1 SUBROUTINES |
|
168
|
|
|
|
|
|
|
|
|
169
|
|
|
|
|
|
|
=head2 parse_from_string |
|
170
|
|
|
|
|
|
|
|
|
171
|
|
|
|
|
|
|
This subroutine takes a string as argument and parses it using |
|
172
|
|
|
|
|
|
|
a Parse::RecDescent parser taken from the package variable |
|
173
|
|
|
|
|
|
|
$Math::Symbolic::Parser. It generates a Math::Symbolic tree |
|
174
|
|
|
|
|
|
|
from the string and returns that tree. |
|
175
|
|
|
|
|
|
|
|
|
176
|
|
|
|
|
|
|
The string may contain any identifiers matching /[a-zA-Z][a-zA-Z0-9_]*/ which |
|
177
|
|
|
|
|
|
|
will be parsed as variables of the corresponding name. |
|
178
|
|
|
|
|
|
|
|
|
179
|
|
|
|
|
|
|
Please refer to L for more information. |
|
180
|
|
|
|
|
|
|
|
|
181
|
|
|
|
|
|
|
=cut |
|
182
|
|
|
|
|
|
|
|
|
183
|
|
|
|
|
|
|
sub parse_from_string { |
|
184
|
357
|
|
|
357
|
1
|
52802
|
my $string = shift; |
|
185
|
357
|
100
|
|
|
|
1561
|
croak "Missing string argument from parse_from_string() call" |
|
186
|
|
|
|
|
|
|
unless defined $string; |
|
187
|
356
|
100
|
|
|
|
1331
|
if ($string eq 'Math::Symbolic') { |
|
188
|
93
|
100
|
|
|
|
323
|
if (@_) { |
|
189
|
92
|
|
|
|
|
221
|
$string = shift; |
|
190
|
|
|
|
|
|
|
} |
|
191
|
|
|
|
|
|
|
else { |
|
192
|
1
|
|
|
|
|
109
|
croak("Missing string argument from Math::Symbolic->parse_from_string() call"); |
|
193
|
|
|
|
|
|
|
} |
|
194
|
|
|
|
|
|
|
} |
|
195
|
355
|
|
|
|
|
1426
|
$string =~ s/\s+//gso; |
|
196
|
355
|
100
|
|
|
|
1319
|
if ( not defined $Parser ) { |
|
197
|
1
|
|
|
|
|
10
|
$Parser = Math::Symbolic::Parser->new(); |
|
198
|
|
|
|
|
|
|
} |
|
199
|
355
|
|
|
|
|
3385
|
return $Parser->parse($string); |
|
200
|
|
|
|
|
|
|
} |
|
201
|
|
|
|
|
|
|
|
|
202
|
|
|
|
|
|
|
1; |
|
203
|
|
|
|
|
|
|
__END__ |