Python symdiff示例

示例#1

文件： parsetest.py 项目： cuihantao/symdiff

def maptest(arg):
  print
  print arg;
  try:
    out = symdiff(arg)
  except SymdiffError as x:
    out = x
  print out

示例#2

文件： parsetest.py 项目： QuantumOfMoose/symdiff

def maptest(arg):
    print
    print arg
    try:
        out = symdiff(arg)
    except SymdiffError as x:
        out = x
    print out

示例#3

文件： parsetest.py 项目： devsim/symdiff

def maptest(arg):
  print("")
  print(arg);
  try:
    out = symdiff(arg)
  except SymdiffError as x:
    out = x
  print(out)

示例#4

def maptest(arg):
    print("")
    print(arg)
    try:
        out = symdiff(arg)
    except SymdiffError as x:
        out = x
    print(out)

示例#5

文件： ordered.py 项目： devsim/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(b, a)')
symdiff('define_model(d, b * c)')
mylist = ordered_list('d');
for i in mylist:
  print(('%s' % i))

示例#6

文件： subexpression1.py 项目： devsim/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(x, y * z)')
symdiff('define_model(z, y * z)')
subexpression()
l = model_list()
for i in l:
  print(("%s, %s" % (i, symdiff('model_value(%s)' % i))))

示例#7

文件： utf8.py 项目： cuihantao/symdiff

# -*- coding: utf-8 -*-
from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

alpha = 'α'
print alpha
x = symdiff('%s' % alpha);
x = symdiff('diff(%s^2, %s)' % (alpha, alpha));
print x

示例#8

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(x, y * z)')
symdiff('define_model(z, y * z)')
subexpression()
l = model_list()
for i in l:
    print(("%s, %s" % (i, symdiff('model_value(%s)' % i))))

示例#9

文件： deftest.py 项目： weatherhead99/symdiff

    "simplify(diff(sin(y),x))",
    "simplify(diff(sin(x^2),x))",
    "simplify(diff(sin(y^2),y))",
    "clear",
    "define(f(x,y),diff(x*y^2+y^3,x),diff(x*y^2+y^3,y))",
    "simplify(diff(f(x,y),x))",
    "simplify(diff(f(x,y),y))",
    "simplify(diff(f(y,x),y))",
    "simplify(diff(f(y,x),x))",
    "clear",
    "declare(g(x,y,z))",
    "define(f(x,y,z),g(y,z,x),1,0)",
    "diff(f(x,y,z),x)",
    "diff(f(a,b,c),x)",
    "diff(f(a,b,c),a)",
    "diff(f(a,b,c),b)",
    "diff(f(a,b,c),c)",
    "clear",
    "a+b",
    "$_*$_",
    "expand($_)",
    "simplify($_)",
    "$a = $_ +1",
    "$b = $_ - $a",
    "$c = simplify($_)",
    "$c = simplify(expand($_))",
]

for i in commands:
    print(symdiff(i))

示例#10

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

print(symdiff('declare_model(x)'))
print(symdiff('define_model(y, a*x)'))
print(symdiff('diff(y, x)'))
print(symdiff('clear_model(x)'))
print(ordered_list('y'))
print(symdiff('declare_model(x)'))
print(ordered_list('y'))
print(symdiff('define_model(x,y)'))
# catch cycles
try:
    print(ordered_list('y'))
except SymdiffError as x:
    print(x)

for i in model_list():
    print(i + ' ' + symdiff('model_value(%s)' % i))
    symdiff('clear_model(%s)' % i)

print("")
print(symdiff('declare_model(x)'))
print(symdiff('define_model(y,x^2)'))
print(symdiff('diff(y,x)'))
print(symdiff('clear_model(y)'))
print(symdiff('model_value(y__x)'))
print(symdiff('diff(model_value(y__x),z)'))

for i in model_list():

示例#11

文件： table.py 项目： cuihantao/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('declare_model(x)')
for i, v in enumerate(symdiff_table('a*x + b*c')):
  print '%s %s' % (i, v)

示例#12

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('declare_model(y)')
symdiff('define_model(x, 3 * y + z)')
symdiff('diff(x, z)')
print('%s' % symdiff('model_value(x)'))
print('%s' % symdiff('model_value(x__z)'))
print('%s' % symdiff('model_value(y__z)'))

示例#13

文件： utf8.py 项目： QuantumOfMoose/symdiff

# -*- coding: utf-8 -*-
from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

alpha = 'α'
print alpha
x = symdiff('%s' % alpha)
x = symdiff('diff(%s^2, %s)' % (alpha, alpha))
print x

示例#14

    "simplify(diff(sin(y),x))",
    "simplify(diff(sin(x^2),x))",
    "simplify(diff(sin(y^2),y))",
    "clear",
    "define(f(x,y),diff(x*y^2+y^3,x),diff(x*y^2+y^3,y))",
    "simplify(diff(f(x,y),x))",
    "simplify(diff(f(x,y),y))",
    "simplify(diff(f(y,x),y))",
    "simplify(diff(f(y,x),x))",
    "clear",
    "declare(g(x,y,z))",
    "define(f(x,y,z),g(y,z,x),1,0)",
    "diff(f(x,y,z),x)",
    "diff(f(a,b,c),x)",
    "diff(f(a,b,c),a)",
    "diff(f(a,b,c),b)",
    "diff(f(a,b,c),c)",
    "clear",
    "a+b",
    "$_*$_",
    "expand($_)",
    "simplify($_)",
    "$a = $_ +1",
    "$b = $_ - $a",
    "$c = simplify($_)",
    "$c = simplify(expand($_))",
]

for i in commands:
    print symdiff(i)

示例#15

文件： deftest.py 项目： cuihantao/symdiff

"clear",
"define(f(x,y),diff(x*y^2+y^3,x),diff(x*y^2+y^3,y))",
"simplify(diff(f(x,y),x))",
"simplify(diff(f(x,y),y))",
"simplify(diff(f(y,x),y))",
"simplify(diff(f(y,x),x))",

"clear",
"declare(g(x,y,z))",
"define(f(x,y,z),g(y,z,x),1,0)",
"diff(f(x,y,z),x)",
"diff(f(a,b,c),x)",
"diff(f(a,b,c),a)",
"diff(f(a,b,c),b)",
"diff(f(a,b,c),c)",

"clear",
"a+b",
"$_*$_",
"expand($_)",
"simplify($_)",
"$a = $_ +1",
"$b = $_ - $a",
"$c = simplify($_)",
"$c = simplify(expand($_))",
]

for i in commands:
  print symdiff(i)

示例#16

文件： remove1.py 项目： weatherhead99/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(x, 0)')
symdiff('define_model(y, x + z)')
remove_zeros()
print(('%s' % symdiff('model_value(y)')))

示例#17

文件： ordered.py 项目： QuantumOfMoose/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(b, a)')
symdiff('define_model(d, b * c)')
mylist = ordered_list('d')
for i in mylist:
    print('%s' % i)

示例#18

文件： arrhenius.py 项目： cuihantao/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(arrhenius, simplify(Cf * pow(T, Etar) * exp(-Ea/T)))')
symdiff('define_model(KC1      , simplify(A3 * R * T * exp((A1 * (-pow(T,(-1)) + pow(A2,(-1))) * pow(R,(-1))))))')
symdiff('define_model(kb_model , simplify(arrhenius / KC1))')
symdiff('declare_model(fprod)')

symdiff('define_model(rprod    , simplify(H * CN * AR))')
symdiff('define_model(omega_dot_H     , simplify((vr_H - vf_H) * (arrhenius * fprod - kb_model * rprod)))')
symdiff('diff(omega_dot_H, T)')
symdiff('diff(omega_dot_H, H)')
symdiff('diff(omega_dot_H, HCN)')
symdiff('diff(omega_dot_H, AR)')
symdiff('diff(omega_dot_H, CN)')

subexpression()
remove_zeros()
mylist = ordered_list('omega_dot_H', 'omega_dot_H__T', 'omega_dot_H__H', 'omega_dot_H__HCN', 'omega_dot_H__AR', 'omega_dot_H__CN')
for i in mylist:
  mv = symdiff("model_value(%s);" % i)
  print "%s, %s" % (i, mv)

示例#19

文件： models1.py 项目： devsim/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

print(symdiff('declare_model(x)'))
print(symdiff('define_model(y, a*x)'))
print(symdiff('diff(y, x)'))
print(symdiff('clear_model(x)'))
print(ordered_list('y'))
print(symdiff('declare_model(x)'))
print(ordered_list('y'))
print(symdiff('define_model(x,y)'))
# catch cycles
try:
  print(ordered_list('y'))
except SymdiffError as x:
  print(x)

for i in model_list():
  print(i + ' ' + symdiff('model_value(%s)' % i))
  symdiff('clear_model(%s)' % i)

print("")
print(symdiff('declare_model(x)'))
print(symdiff('define_model(y,x^2)'))
print(symdiff('diff(y,x)'))
print(symdiff('clear_model(y)'))
print(symdiff('model_value(y__x)'))
print(symdiff('diff(model_value(y__x),z)'))

for i in model_list():

示例#20

文件： arrhenius.py 项目： weatherhead99/symdiff

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('define_model(arrhenius, simplify(Cf * pow(T, Etar) * exp(-Ea/T)))')
symdiff(
    'define_model(KC1      , simplify(A3 * R * T * exp((A1 * (-pow(T,(-1)) + pow(A2,(-1))) * pow(R,(-1))))))'
)
symdiff('define_model(kb_model , simplify(arrhenius / KC1))')
symdiff('declare_model(fprod)')

symdiff('define_model(rprod    , simplify(H * CN * AR))')
symdiff(
    'define_model(omega_dot_H     , simplify((vr_H - vf_H) * (arrhenius * fprod - kb_model * rprod)))'
)
symdiff('diff(omega_dot_H, T)')
symdiff('diff(omega_dot_H, H)')
symdiff('diff(omega_dot_H, HCN)')
symdiff('diff(omega_dot_H, AR)')
symdiff('diff(omega_dot_H, CN)')

subexpression()
remove_zeros()
mylist = ordered_list('omega_dot_H', 'omega_dot_H__T', 'omega_dot_H__H',
                      'omega_dot_H__HCN', 'omega_dot_H__AR', 'omega_dot_H__CN')
for i in mylist:
    mv = symdiff("model_value(%s);" % i)
    print("%s, %s" % (i, mv))

示例#21

from symdiff import *
# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff('declare_model(x)')
for i, v in enumerate(symdiff_table('a*x + b*c')):
    print '%s %s' % (i, v)

示例#22

文件： subexpression1.py 项目： weatherhead99/symdiff

from symdiff import *

# Any copyright is dedicated to the Public Domain.
# http://creativecommons.org/publicdomain/zero/1.0/

symdiff("define_model(x, y * z)")
symdiff("define_model(z, y * z)")
subexpression()
l = model_list()
for i in l:
    print("%s, %s" % (i, symdiff("model_value(%s)" % i)))