- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我想在 sympy 中评估类似形式(更复杂)的函数。
y = a * b / np.sum( a*( b + c) )
其中所有变量都是长度为n
的向量。评估将在优化例程的每个时间步进行。因此,我想有效地实现它。最有可能的是,最好编译这些函数,但 autowrap 模块给我带来了奇怪的错误。
什么有效:
import numpy as np
import sympy as sp
from __future__ import division
a = sp.IndexedBase('a')
b = sp.IndexedBase('b')
c = sp.IndexedBase('c')
n = 4
expr_fun = lambda x: a[x] * (b[x] + c[x])
expr = [ a[i]*b[i] / np.sum([expr_fun(i) for i in range(n)]) for i in range(n)]
我可以直接在 sympy 中计算这个表达式:
r = np.random.random(n)
subs_dict = {}
[ subs_dict.update({a[i]:r[i],b[i]:r[i] }) for i in range(n) ]
[expr[i].subs(subs_dict) for i in range(n)]
给我(如预期):
0.0786923966864026*c[0] + 0.403977159637609*c[1] + 0.598011208186539*c[2] + 0.0896229978341944*c[3] + 0.535039725662632
但是我编译这个表达式失败了。我从几个小时开始就在阅读博客和手册,但要么我太累了,要么没有找到正确的信息。非常感谢任何帮助。
编辑:回应 Eric:我现在不知道如何在 theano 或 autowrap 中实现向量的总和。我使用 lambda 函数尝试了不同的版本,但遇到了各种错误。也许最可重复的与输入的维度有关:
from sympy.printing.theanocode import theano_function
from sympy.printing.theanocode import sympy as sp
from sympy.printing.theanocode import dim_handling
import numpy as np
symbols = ['a', 'b', 'c', 'd']
a, b, c, d = map(sp.Symbol, symbols )
expr = a + b*(c+d)/np.sum(b + c*d)
n = 1
dim = {} # collections.OrderedDict()
[dim.update( {i: n} ) for i in [a, b, c, d] ]
dt = {} # collections.OrderedDict()
[dt.update( {i: 'float64'} ) for i in [a, b, c, d] ]
f = theano_function( [a, b, c, d], [expr], dims = dim, dtypes=dt )
in_var = np.array([ [1,2,3,4] ])
f(in_var.T)
TypeError: ('Bad input argument to theano function at index 0(0-based)', 'Wrong number of dimensions: expected 1, got 2 with shape (4, 1).')
如果我尝试使用自动换行编译一个简单的表达式:
import sympy as sp
from sympy.utilities.autowrap import autowrap
m, n = sp.symbols('m n', integer=True)
a, b, c,d = map(sp.IndexedBase, ['a', 'b', 'c', 'd'])
i = sp.Idx('i',m)
j = sp.Idx('j',n)
instruction = sp.Eq(a[i], b[i]*(c[i] + d[i]) )
f = autowrap(instruction)
我得到:
UnicodeDecodeError: 'ascii' codec can't decode byte 0xe2 in position 3749: ordinal not in range(128)
最佳答案
您的方法存在多个问题。 1)您不能将 numpy 函数与 SymPy 类型混合,然后期望代码生成器能够处理混合类型。 SymPy 中的代码生成器仅适用于 SymPy 类型。您可以传入映射到 SymPy 函数的外部函数。例如,您可以使用 sympy.Sum(),然后编写从 SymPy.Sum() 到 numpy.sum() 的映射以供代码生成器使用。 2)据我所知,代码生成器不支持 Eq() 。 3) 索引类型在传递给 autowrap 时会执行非常具体的操作。您需要仔细阅读有关它们的文档。
关于python - 如何使用 theano 或 autowrap 评估和编译这样的函数?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/29723803/
是否有任何库或框架旨在促进从另一种成熟的编程语言中构建项目? 在 C++、java 等编程语言中指定逻辑、集合和复杂规则非常容易,但在 Makefile 中完成这些事情似乎是一场艰苦的战斗。我还没有深
我有这段代码可以用 clang 编译得很好(即使使用 -Weverything),但是 gcc 会发出错误。 #include #include #include using namespace
我有以下 block 头文件 BKE_mesh.h: /* Connectivity data */ typedef struct IndexNode { struct IndexNode *
我在我的一个项目中遇到了一个奇怪的问题。我的代码库依赖于一个外部库,其中包含一个名为 Dataset 的类. Dataset类私有(private)继承自 std::vector (其中 Sample
当使用 gcc、g++ 或 make 在终端中编译一个小型 C 或 C++ 项目时,我收到以下错误: /tmp/ccG1caGi.o: In function `main': main.c:(.tex
我正在尝试从 CVS 为 Windows 上的 Emacs 23.1.50 编译 CEDET,但在“第 6 步:打开 EDE...”时出现错误:“defvar:作为变量的符号值是无效的:cedet-m
我正在(重新)学习编程,我从 C 开始。我的 IDE(如果我可以这么说)是 Windows7 上的 cygwin(32 位)和 Visual-Studio 2010。我总是编译我用 gcc (cygw
我喜欢在模板类中使用本地类来执行类似“static if”的构造。但是我遇到了 gcc 4.8 不想编译我的代码的问题。但是 4.7 可以。 这个例子: #include #include #in
我有一个项目,必须仅使用 java 1.4 进行编译。但我计划使用mockito 编写一些单元测试。我想要一种在 pom 中指定的方法,以便 src/main/java 使用 jdk 1.4 编译,但
我想了解 PHP 编译过程是如何工作的。 假设我有一个名为funcs.php 的文件并且这个文件有三个函数,如果我include 或require 它,所有的在文件加载期间编译三个函数?或者源代码会被
编译工具链 我们写程序的时候用的都是集成开发环境 (IDE: Integrated Development Environment),集成开发环境可以极大地方便我们程序员编写程序,但是配置起来
当我编写一些 Scala 代码时,在尝试编译代码时收到一条奇怪的错误消息。我将代码分解为一个更简单的代码(从语义的角度来看这完全没有意义,但仍然显示了错误)。 scala> :paste // Ent
我正在编译一个 SCSS 文件,它似乎删除了我的评论。我可以使用什么命令来保留所有评论? >SASS input.scss output.css 我在 SCSS 中看到两种类型的注释。 // Comm
这是我的代码: #include typedef struct { const char *description; float value; int age; } swag
当您编译 grails war 时,我知道 .groovy 代码被编译为字节码类文件,但我不明白容器(例如 tomcat)如何在请求 GSP 时知道如何编译它们。容器了解 GSP 吗?安装在服务器上的
我正在努力将多个文件编译成一个通用程序。我收到一个错误: undefined reference to 'pi' 这是我的代码和 Makefile 的框架。我做错了什么?谢谢! 文件:calcPi.c
我尝试使用 LD_PRELOAD 来 Hook sprintf function ,所以我将打印到缓冲区的结果: #define _GNU_SOURCE #include #include int
我正在寻找最简单的方法来自动将 CoffeeScript 重新编译为 JS。 阅读documentation但仍然很难得到我想要的东西。 我需要它来监视文件夹 src/ 中的任何 *.coffee 文
我想使用定制waveformjs 。我发现this on SO但是,我不知道如何编译/安装波形来开始。我从 GitHub 克隆它并进行了更改,但是我不知道如何将其转换为 .js 文件。 最佳答案 为了
很难说出这里问的是什么。这个问题是含糊的、模糊的、不完整的、过于宽泛的或修辞性的,无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开它,visit the help center 。 已关
我是一名优秀的程序员,十分优秀!