PyMC3-自定义theano Op进行数值积分
PyMC3是一个用于概率编程的Python库，它提供了一种灵活且高效的方式来建立和推断概率模型。PyMC3基于Theano，一个用于定义、优化和评估数学表达式的Python库。
数值积分是在数学和物理领域中常见的一种计算方法，用于计算函数在给定区间上的面积或体积。PyMC3允许用户通过自定义Theano Op来进行数值积分。
自定义Theano Op是一种用户可以定义和实现的Theano操作符。用户可以使用Python编写自定义Theano Op来实现特定的数值积分算法。通过将自定义Theano Op与PyMC3的概率模型结合使用，用户可以在概率模型中进行数值积分。
使用PyMC3进行数值积分的优势包括：
- 灵活性：PyMC3允许用户根据具体需求自定义数值积分算法，从而满足不同的数学模型和问题。
- 高效性：PyMC3基于Theano，利用了Theano的优化和并行计算功能，可以提高数值积分的计算效率。
- 可扩展性：PyMC3支持并行计算和分布式计算，可以在多个计算节点上进行数值积分，从而提高计算速度和处理大规模数据的能力。
数值积分在很多领域都有广泛的应用场景，例如：
- 物理学：数值积分可以用于计算物体的质心、惯性矩阵等物理量。
- 金融学：数值积分可以用于计算金融衍生品的价格和风险。
- 机器学习：数值积分可以用于计算概率模型的边缘概率和条件概率。
腾讯云提供了一系列与云计算相关的产品，其中包括计算、存储、数据库、人工智能等多个领域。对于PyMC3中自定义Theano Op进行数值积分的应用，可以使用腾讯云的云服务器（CVM）来搭建计算环境，使用云数据库（TencentDB）来存储数据，使用人工智能平台（AI Lab）来进行模型训练和推理。
更多关于腾讯云产品的介绍和详细信息，请参考腾讯云官方网站：https://cloud.tencent.com/
相关·内容
1回答
我试图实现一个自定义的Op，但我得到了"Bad input argument to theano function“错误。下面是代码。据我所知，问题是:如何将PyMC3变量转换为theano可理解的类型？
import numpy as np
import theano
import theano.tensor as t
from theano import config
config.compute_test_value = 'off'
#true_Data = [1,2]
#values=[]
class trial_Op(theano.Op):
__p
浏览 0提问于2017-02-22得票数 0
1回答
我正在尝试使用pymc3定义一个复杂的自定义似然函数。似然函数涉及大量迭代，因此我尝试使用theano的扫描方法直接在theano中定义迭代。这里有一个非常简单的例子，它说明了我所面临的挑战。我试图定义的(伪)似然函数是两个pymc3随机变量p和θ的和。当然，我可以简单地返回p+theta，但我尝试编写的实际似然函数更复杂，而且我认为我需要使用theano.scan，因为它涉及很多迭代。
import pymc3 as pm
from pymc3 import Model, Uniform, DensityDist
import theano.tensor as T
import theano
浏览 10提问于2016-09-03得票数 0
我试图使用在Theano中内置的函数。但是，当我从theano.function中定义ExpmGrad时，我会看到一个错误，即输出必须由theano变量来完成。
我不确定使用这个ExpmGrad函数的正确方法是什么，因为我没有在网上找到任何使用它的例子。
这就是我试过的：
import theano
from theano.tensor import T
J1 = T.dscalar('J1')
H = np.arange(16).reshape(4, 4) * J1
gJ = theano.tensor.slinalg.ExpmGrad(H)
f = theano.funct
浏览 4提问于2016-12-01得票数 0
回答已采纳
3回答
我正在尝试编写一个自定义Theano Op，它以数字方式集成两个值之间的函数。Op是PyMC3的一种自定义似然，它涉及到一些积分的数值计算。我不能简单地使用@as_op装饰器，因为我需要使用HMC来完成HMC步骤。任何帮助都将不胜感激，因为这个问题似乎已经出现了几次，但从未得到解决(例如，)。
显然，一种解决方案是在Theano中编写一个数值积分器，但当已经有非常好的积分器可用时，例如通过scipy.integrate，这似乎是浪费精力。
为了保持这是一个最小的例子，让我们试着在一个Op中集成一个介于0和1之间的函数。下面在Op之外集成一个Theano函数，并在我的测试过程中产生正确的结果。
浏览 4提问于2017-03-09得票数 6
回答已采纳
给定几个要获取的符号变量，我需要知道哪些占位符是依赖的。
在Theano，我们有：
import theano as th
import theano.tensor as T
x, y, z = T.scalars('xyz')
u, v = x*y, y*z
w = u + v
th.gof.graph.inputs([w]) # gives [x, y, z]
th.gof.graph.inputs([u]) # gives [x, y]
th.gof.graph.inputs([v]) # gives [y, z]
th.gof.graph.inputs([u,
浏览 0提问于2017-08-08得票数 3
我创建了一个theano.Op，它返回两个输入集合的每对之间的距离，转换scipy cdist
class Cdist(theano.Op):
__props__ = ()
def __init__(self):
#self.fn = scipy_cdist2
super(Cdist, self).__init__()
def make_node(self, x, w):
#print('make_node')
return gof.Apply(self, [x, w], [x.typ
浏览 3提问于2017-06-29得票数 1
回答已采纳
1回答
我正在尝试定义一个带有梯度的自定义theano操作，以便与pymc3一起使用它，但我不明白如何定义grad方法。
下面的代码是我卡住的地方。函数phi()是一个模拟函数(实际上，它是一个外部程序)；对于标量输入x，它返回一个向量(phi_0(x), phi_1(x), ...)。函数phi_diff() (也是一个模拟函数)返回向量(dphi_0/dx, dphi_1/dx, ...)。
我将phi()和phi_diff()封装在theano.Op对象中，但grad函数的实现不起作用。theano的文档包含了一些简单的例子，我不明白在这种情况下如何调整它们。任何帮助都将不胜感激。
import
浏览 0提问于2018-10-11得票数 1
回答已采纳
1回答
import os
import theano, numpy
from theano import tensor as T
from collections import OrderedDict
class RNNSLU(object):
""" Elman neural net"""
def __init__(self, nh, nc, ne, de, cs):
"""
Hyperparameters used for initialization
nh : dimension of the hidden lay
浏览 2提问于2015-07-29得票数 0
回答已采纳
1回答
Matlab中数值积分的精度
、、、、
我正在尝试在矩形区域上集成一个解析函数( sqrt和trig函数的组合)。它在该区域没有奇异点，似乎是使用dblquad的完美候选者。我的问题是如何评估Matlab提供给我的数值的准确性。在不知道积分的确切值的情况下，我们如何证明有效数字的合理性？当需要给出一个精度为特定位数的值时，您应该能够对齐。考虑到这个值是用Matlab计算的，有可能实现吗？
浏览 4提问于2013-05-16得票数 4
b=T.fvector()
theano.gradient.hessian( T.gammaln(b).sum(), [b] )
给出一个未命名的NotImplementedError。由于等价的代码在T.log而不是T.gammaln上运行良好，我猜theano只是不知道gammaln的第二个派生
浏览 13提问于2017-02-02得票数 0
回答已采纳
3回答
我是提亚诺的新手，我试图实现一个反应扩散系统的数值积分器--这个版本的：
目前，我的表达方式是：
import theano as th
import theano.tensor as T
u = T.dmatrix('u')
v = T.dmatrix('v')
e = T.dscalar('e')
a0 = T.dscalar('a0')
a1 = T.dscalar('a1')
dudt = u - u**3 -v
dvdt = e*(u - a1*v - a0)
所以我还没有实现有限
浏览 7提问于2016-02-21得票数 8
回答已采纳
我试图计算python中高斯的积分，如下所示：
from math import exp
from scipy import stats, integrate
import scipy.interpolate as interpolate
from numpy import cumsum, random, histogram, linspace, zeros, inf, pi,sqrt
import matplotlib.pyplot as plt
A = 1
mu = 0
sigma = 1
p = lambda x: A * exp(-(((x-mu)**2))/(2*(sigma**2
浏览 1提问于2015-08-13得票数 0
回答已采纳
1回答
我最近在Sympy中遇到了一些符号矩阵导数的性能瓶颈(具体地说，通过使用lambdas进行替换来计算符号矩阵的单行代码占用了大约90%的程序运行时间)，因此我决定尝试Theano。
它以前的应用是计算高斯过程超参数上的偏导数，其中使用Sympy symbols (MatrixSymbol)的(1，k)维矩阵在迭代该列表和区分每个项目上的矩阵方面工作得很好。
然而，这并没有延续到Theano中，文档似乎也没有详细说明如何做到这一点。在Theano中索引符号向量将返回Subtensor类型，该类型对于计算上的梯度无效。
下面是我正在尝试做的事情的一个简单版本(但在算法上完全不正确--精简到我试图获
浏览 9提问于2016-07-27得票数 0
当我想通过RCaller运行一些代码时，Java会抛出以下异常：
Exception in thread "JavaFX Application Thread" com.github.rcaller.exception.ExecutionException: Can not send the source code to R file due to: java.io.IOException: The pipe is being closed Maximum number of retries exceeded.
这是我的密码：
protected void initialize
浏览 2提问于2017-06-05得票数 2
回答已采纳
我未能将属于类实例的方法作为确定性函数与PyMc3相匹配。你能告诉我怎么做吗？
为了简单起见，下面用一个简单的例子对我的情况进行了总结。在现实中，我的约束是，一切都是通过图形用户界面完成的，像‘find_MAP’这样的操作应该在链接到pyqt按钮的方法中。
我想在数据点上拟合“FunctionIWantToFit”函数。问题，请使用以下代码：
import numpy as np
import pymc3 as pm3
from scipy.interpolate import interp1d
import theano.tensor as tt
import theano.compile
浏览 5提问于2017-02-13得票数 8
回答已采纳
1回答
3d滑动窗口操作在Theano？
、、、、
TL.DR。是否有的三维友好实现
我想使用神经网络对体积进行体素分类( nxnxn )，该神经网络接收nxnxn图像，其中N>n。要对体积中的每个体素进行分类，我必须迭代每个体素。对于每次迭代，我获得邻域体素并将其作为输入传递给神经网络。这只是一个简单的滑动窗口操作，该操作是神经网络。
虽然我的神经网络是用Theano实现的，但滑动窗口实现是用python/numpy实现的。由于这不是纯粹的Theano操作，因此分类需要永远(> 3小时)来将所有体素分类到一个体积中。对于二维滑动窗口操作，Theano有一个辅助方法，theano.tensor.nnet.neighbours.ima
浏览 39提问于2016-02-29得票数 21
下面的代码是从在线教程中借来的。我看到下面的行写在代码的主要方法中
c = broadcasted_add(a, b)
增加维数的张量a (2,1,2,2)和维数的张量b(2，2，2)。即使我们在make_tensor方法中将可广播性声明为“false”，它如何能够正确地添加？难道我们不应该宣布广播为真，这样它就可以添加不同的尺寸？它不应该抛出一个错误，说尺寸不匹配？我对广播的理解是错误的吗？
import numpy as np
from theano import function
import theano.tensor as T
def make_tensor(dim):
&
浏览 1提问于2015-10-13得票数 3
回答已采纳
1回答
如何在PyMC3中定义自定义似然？在PyMC2中，我可以使用@pymc.potential。我尝试在PyMC3中使用pymc.Potential，但是，似乎布尔操作不能应用于参数(当我这样做时，我得到一个类似的错误)。例如，以下代码不起作用：
from pymc import *
with Model() as model:
x = Normal('x', 1, 1)
def z(u):
if u > 0: #comparisons like this are not supported
# if theano.tens
浏览 4提问于2014-12-18得票数 7
2回答
计算密度估计曲线下面积，即概率
、、、、
我对我的数据曲线(见下图)有一个密度估计(使用density函数)，我需要找到概率Pr(learningTime > c)，即密度曲线下的面积从给定的数字c (红色垂直线)到曲线末端。有什么想法吗？
浏览 2提问于2016-11-29得票数 8
回答已采纳
2回答
我试图计算特定区域上二元正态分布的概率，这是java中的一个特定多边形。
数学描述是将特定复杂区域上二元正态分布的概率密度函数(、pdf、)结合起来。
我的第一种方法是在NormalDistribution库的帮助下使用两个apache-commons-math对象。给定维度1的dataset x和维度2的dataset y，我计算了每个NormalDistribution的均值和标准差。
使用public double probability(double x0, double x1) from org.apache.commons.math3.distribution.Norma
浏览 13提问于2015-04-29得票数 1
设D是以一系列点{x_i，y_i} (1<=i<=N)为界的区域。(该区域不一定是凸的，且这些点沿边界曲线移动。)
设f是定义在D上的函数，但我们只知道它在给定的点集(有限和离散)上的值，例如{x'_i，y'_i，f(x'_i，y'_i)} (1<=i<=N‘)(给定的数据集在D中不需要“稠密”)。
我怎样才能对D上的f进行数值积分呢？
以下是我的想法：
( 1)首先，我们应该用一系列点之间的分段来逼近D的边界。
2)然后对给定的数据集进行插值。然而，，在二维插值并不总是可能的，然后我被卡住了.
你能帮忙吗？谢谢。
浏览 1提问于2011-05-17得票数 1
我试图在第23章的贝叶斯数据分析方法的基础上，使用PyMC3对一个序数预测变量进行建模。我希望使用find_MAP确定一个良好的起始值，但是正在收到一个优化错误。
模式：
import pymc3 as pm
import numpy as np
import theano
import theano.tensor as tt
# Some helper functions
def cdf(x, location=0, scale=1):
epsilon = np.array(1e-32, dtype=theano.config.floatX)
location = tt.
浏览 1提问于2016-11-04得票数 0
1回答
我在使用pymc2复制pymc3代码时遇到了问题。
我相信这是因为pymc3使用的是theano类型的变量，这些变量与我使用的numpy操作不兼容。所以我使用的是@theano.decorator：
我的职能是：
with pymc3.Model() as model:
z_stars = pymc3.Uniform('z_star', self.z_min_ssp_limit, self.z_max_ssp_limit)
Av_stars = pymc3.Uniform('Av_star', 0.0,
浏览 4提问于2017-09-29得票数 3
回答已采纳
6回答
集成的C数学库？
、、
是否有与科学/数学相关的C/C++编程库，其中包含数学集成函数？这可能有助于查找累积分布函数的值，等等。符号操作是理想的，但不是必需的。谢谢!
附注:我看过GSL，但它不包含特定于集成的功能。
浏览 2提问于2009-10-14得票数 4
回答已采纳
1回答
对MATLAB集成的误区
、、、
所有人！
你们能不能解释一下为什么精确集成的答案
quad( @(x)x.^3-2.*x,0,2)
是
-2.220446049250313e-016
而不是零，就像它应该是的那样？
浏览 1提问于2011-11-02得票数 0
回答已采纳
我试图重构我的代码，以便更容易地更改架构。目前，我正在构建一个递归神经网络，如下所示。
# input (where first dimension is time)
x = T.matrix()
# target (where first dimension is time)
t = T.matrix()
# recurrent weights as a shared variable
W_hh = theano.shared(numpy.random.uniform(size=(n, n), low=-.01, high=.01))
# input to hidden layer wei
浏览 2提问于2014-12-11得票数 0
回答已采纳
2回答
C#中的定积分计算
、、、
如何更好地计算定积分？我使用一个函数进行集成，另一个函数递归地找到阶乘。
我想要更好的算法，效率，甚至准确性。
public static double testStatistic(double meanTreatmentSumOfSquares, double meanErrorSumOfSquares)
{
return (meanTreatmentSumOfSquares / meanErrorSumOfSquares);
}
public static double pValue(double fStatistic, int degre
浏览 0提问于2012-07-25得票数 3
1回答
我遵循theano教程给出了简单随机梯度下降的。但是，在这里，我无法理解p_y_given_x和y_pred的值是如何根据W和b的值自动更新的，因为以后运行test_logistic()时，我们只更新W和b的值？谢谢
class LogisticRegression(object):
def __init__(self, input, n_in, n_out):
self.W = theano.shared(
value=numpy.zeros(
(n_in, n_out),
dt
浏览 1提问于2015-01-18得票数 1
回答已采纳
1回答
我目前正在进行一个项目，使用opencv和python来测量通常是弯曲的对象，例如，下面所示的箭头尽可能准确。
我认为一种策略可能是使用枕叶Voronoi函数来获得箭头中心脊柱的点，但现在有麻烦了。这是我的代码：
img = cv2.imread('example_rubystreak_2.PNG')
img.shape
gray = cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
ret, thresh = cv2.threshold(gray,190,255,cv2.THRESH_BINARY)
countimage, contours,
浏览 10提问于2015-08-07得票数 2
回答已采纳
我是西亚诺的新手。我试图实现简单的线性回归，但我的程序抛出了以下错误：
"/home/akhan/Theano-Project/uog/theano_application/linear_regression.py:36“：(‘TypeError函数的输入参数错误，名称为0(基于0)’，”期望有一个类似数组的对象，但是找到了一个变量:也许您试图在一个(可能是共享的)变量上调用一个函数，而不是一个数字数组？“)
这是我的代码：
import theano
from theano import tensor as T
import numpy as np
import matplo
浏览 3提问于2015-05-27得票数 3
回答已采纳
1回答
我有一个密度函数f_N，其定义如下(K_nu(z)是修改后的贝塞尔函数)：
我想为N的每个值计算下面的积分
以下是上述方法在matlab中的实现。
for N=1:100
syms z
f =@(z) (1/(gamma(N)*sqrt(pi))*(z/2).^(N-0.5).*besselk(0.5-N,z));
g = @(z) f(z).*log(f(z));
val=integral(g,0,Inf);
但当我运行上面的代码时，它总是为N的各种值返回NaN，并显示以下警告：
Warning: Infinite or Not-a-Num
浏览 42提问于2017-02-05得票数 0
我正在研究一个非线性联立方程组。这两个变量是u>0和b>0。如何在Matlab、Python或Fortran中解决此问题？谢谢。
浏览 2提问于2012-05-02得票数 0
回答已采纳
3回答
我试着用数值计算一个双积分。问题是，(我想)我需要符号积分和数值积分的混合。积分看起来是这样的：
我不能使用numpy.integrate，因为它不只是一个双积分，因为它的幂(1/a)在中间。我无法得到最里面的积分(然后提高到幂)的数字，因为它最终是一个依赖于x的函数，然后我需要对它进行积分。
我尝试使用符号演算，使用嵌套的sym.integrate，如下所示
sym.integrate((sym.integrate(sym.exp(-(w**2)/(2*sigmaw)-alpha*((x-w)**2)/(2*sigma)),(w,-sym.oo, sym.oo)))**(1/alpha),
浏览 6提问于2022-01-18得票数 0
我有一个关于加速度a( t )的函数，如果我对t积分两次，就可以找到位置x(t)。我正试图找到t= 10秒的位置。
因为a(t)不是多变量，所以我在使用Scipy dblquad函数来计算我需要的二重积分时遇到了问题。请看我到目前为止所拥有的：
def a(t):
return (2.5 / (1 + math.exp((t-8)/0.8)))
def upperbound():
return 10
def lowerbound():
return 0
x = dblquad(a,0,10,lowerbound,upperbound)
浏览 3提问于2017-10-26得票数 3
1回答
西亚诺有操作和功能。
有什么区别？
函数似乎很好，很容易定义，如：
x = T.dmatrix('x')
linmax = function([x], T.maximum(x,0))
行动的定义似乎很复杂。所有抽象类和类似的东西(如theano.tesnor.tanh和theano.tensor.nnet.sigmoid )都被定义为Ops。我不确定有什么区别。
我将如何将上面的linmax函数写成一个操作程序？
浏览 1提问于2013-11-27得票数 2
回答已采纳
我想从一个我不知道归一化常数的分布中画出随机偏差。分布是未知形状和scale=1的gamma似然的共轭先验。
这个pdf的积分没有封闭的形式。在这里使用scipy.stats.rv_continuous是一个选项吗？如果不是，任何人还能推荐什么来绘制任意pdf的偏差？
浏览 1提问于2016-04-29得票数 0
我想在matlab中计算int[0,10]int[0,y] f(x,y) dxdy。当然，我可以使用integral2或quad2d：integral2(f,0,10,0,@(x)x)。但是，我也可以处理g(y)=nt[0,y] f(x,y) dx并做积分int[0,10] g(y) dy。我不想使用integral2，而是只想使用integral2进行计算。我写的代码是：
g=@(y)integral(@(x)f(x,y),0,y)
integral(g,0,10)
但是，它给出的错误消息是
Error using integral (line 85) A and B must be flo
浏览 18提问于2016-07-26得票数 0
1回答
我正在使用R来做一些多变量分析。对于这项工作，我需要对三变量PDF.Since进行积分，我想在最大似然估计中使用它，一个积分向量。有没有办法让Integrate带来一个向量而不是一个值。
下面是一个简单的例子：
f1=function(x, y, z) {dmvnorm(x=as.matrix(cbind(x,y,z)), mean=c(0,0,0), sigma=sigma)}
f1(x=c(1,1,1), y=c(1,1,1), z=c(1,1,1))
integrate(Vectorize(function(x) {f1(x=c(1,1,1), y=c(1,1,1), z=c(1,1,1
浏览 0提问于2013-04-18得票数 0
回答已采纳
1回答
我需要用MATLAB集成一个表达式。这些功能之间有什么区别？
QUAD，QUAD8，QUADL，DBLQUAD
浏览 2提问于2011-02-23得票数 0
我需要通过扫描多次执行theano函数，以便对成本函数求和并在梯度计算中使用它。我熟悉做这件事的深度学习教程，但是我的数据切片和其他一些复杂性意味着我需要做一些不同的事情。下面是我正在尝试做的事情的一个简化版本。
tn = testnet()
cost = tn.single_cost( )
x = theano.shared(numpy.asarray([7.1,2.2,3.4], dtype='float32'))
index = T.lscalar('index')
test_fn = theano.function(inputs=[index], out
浏览 1提问于2015-04-01得票数 5
1回答
我想知道R中是否有一个函数可以求解下面的方程式中的x
P(X < x，Y< a) = b，
其中(X，Y)遵循二元正态分布或(斜)学生t分布，a和b是已知的。
换句话说，我试图计算二重积分的上界x。
浏览 4提问于2016-01-05得票数 1
2回答
我正在阅读Bergstra等人的 2012，并且不理解在这种情况下“符号”的含义(例如，指向目录的符号链接)：
在Python循环中使用Theano的symbolically-defined实现可以通过迭代过程阻止符号的区分，并阻止应用某些图优化。将循环完全展开为符号链通常会导致不可管理的大图，并且不允许具有可变迭代次数的“while”-style循环。
浏览 2提问于2015-09-24得票数 2
回答已采纳
我试着用conv3d2d制作3d CNN。下面出现错误(我的代码位于)：
回溯(最近一次调用)：
文件"/home/shome/workspace_temp/conv3d_test/conv3d_test.py"，第1行，从convnet3d导入ConvLayer
文件"/home/shome/softwares/theano/Theano-3D-ConvNet-master/convnet3d/convnet3d.py"，第12行，从conv3d2d导入conv3d
文件"/home/shome/softwares/theano/Theano-3D-
浏览 2提问于2016-01-19得票数 0
回答已采纳
我已经写了一些python或matlab的神经网络代码，但没有使用任何框架或自动微分，但正如我们所知，Theano和TensorFlow使用自动微分，你建立一个计算图，他们为你做计算(反向传播)，但有时写的程序可以运行，但肯定不像我希望的那样运行，所以我想知道有一些方法来确保我的程序是正确的？打印构造的计算图？但是当NN层的数量很大时似乎很复杂，比如Imagenet的获胜者采用了152层，或者使用简单的matlab或python代码编写另一个程序，然后将此输出与使用框架的程序进行比较？
浏览 0提问于2016-02-19得票数 1
1回答
我想仔细看看全连接神经网络中每一层的雅可比矩阵，即∂y/∂x，其中x是该层的输入向量(激活上一层)，y是该层的输出向量(激活这一层)。
在在线学习方案中，这可以很容易地完成如下操作：
import theano
import theano.tensor as T
import numpy as np
x = T.vector('x')
w = theano.shared(np.random.randn(10, 5))
y = T.tanh(T.dot(w, x))
# computation of Jacobian
j = T.jacobian(y, x)
当学习批处理时，
浏览 0提问于2017-06-02得票数 1
4回答
我想解决这个问题：
我不想用"int"，我想用"quad“家族(quad，dblquad，triplequad)，但我不能。你能帮我吗？
浏览 0提问于2010-07-27得票数 0
如何在Mathematica中数值计算二重积分？
Integrate[Exp[-0.099308 s]
* Integrate[Exp[0.041657423 u] Exp[-3.1413 s + 3.12 u]
* ((u/(s - u))^(1/2) BesselI[1,2 (u (s - u))^(1/2)]
+ 0.293 BesselI[0,2 (u (s - u))^(1/2)]),{u,0,s}],{s,0,10}]
浏览 0提问于2012-03-06得票数 3
回答已采纳
不言而喻--根据的说法，通用的集成方法是“四边形”。命名约定继续适用于双重集成(dblquad)、三重集成(tplquad)等。
据我所知，积分和术语"quad“之间没有数学上的联系。那么，为什么scipy会这样命名这些方法呢？
浏览 1提问于2020-04-09得票数 0
在使用joblib并行化一些涉及Theano函数的模型拟合代码的过程中，我偶然发现了一些对我来说很奇怪的行为。
考虑下面这个非常简单的例子：
from joblib import Parallel, delayed
import theano
from theano import tensor as te
import numpy as np
class TheanoModel(object):
def __init__(self):
X = te.dvector('X')
Y = (X ** te.log(X ** 2)).sum()
浏览 0提问于2014-08-18得票数 7
回答已采纳
1回答
我正试图在之后的两个GPU中并行处理NN。我有所有的依赖项，但是cuda运行时初始化失败了，下面的消息如下。
ERROR (theano.sandbox.cuda): ERROR: Not using GPU. Initialisation of device 0 failed:
cublasCreate() returned this error 'the CUDA Runtime initialization failed'
Error when trying to find the memory information on the GPU: invalid device
浏览 1提问于2015-09-24得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
