我正在使用fipy来解决一个计算域中具有孤立区域的扩散问题。请参见示意图。下图:
在fipy.FaceVariable的帮助下,在@Daniel Wheeler的帮助下,该问题在Fipy下进行了建模,并使用了fipy.FaceVariable定义的可变系数。
然而,计算速度无法满足我的命令,这比使用有限差分法的cython代码要慢得多。如果我想加速fipy计算,我该怎么做?以下是我的演示代码:
from pylab import *
from numpy import *
import fipy
from scipy.spatial import Delaunay
from fipy.vari
首先,如果我犯了一个愚蠢的错误,很抱歉--这是我的第一个tensorflow程序之一。为了我的辩护--我已经搜查过这件事,没有找到任何合理的解释。
密码
代码是从改编的。它使用一个梯度来训练一个简单的线性模型:
import tensorflow as tf
# Define a linear model. The initial values are wrong (W=0.3, b=0.3),
# and should converge to (W=-1, b=1).
# Variables
W = tf.Variable([.3], tf.float32)
b = tf.Variabl
我正在尝试在Matlab中模拟Kuramoto的隐含。我试着用ode45来解决这个问题。我还看到其他人使用Runge-kutta方法。我知道ode45使用龙格-库塔方法,然而,我从每个方法中获得的值有可疑的不同。
kuramoto= @(x,K,N,Omega)Omega+(K/N)*sum(sin(x*ones(1,N)-(ones(N,1)*x')))'
%Kuramoto is a model of N coupled ocilators (such as multiple radiowaves)
%The solution to the model is the pha
我想知道odeint中的步长是否固定。在中
基本的步进概念。遵循这一步进概念的基本步骤能够执行ODE的解x(t)的一步,从而使用给定的步长dt获得x(t+dt)。
在下面的代码中,
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
/* The type of container used to hold the state vector */
typedef std::vec
我最近了解了的特性。我开始在bsplines和一个显式组件中使用这个特性。我目前的问题是,我的设计变量是从两个不同的组件输入的,类似于下面的xsdm。据我所知,不可能为不同的设计变量设置不同的有限差分步骤。因此,再次查看xsdm,控制点x和z应该具有相同的FD步骤,即
model.approx_totals(step=1)
有效,但是
model.approx_totals(step=np.ones(5))
行不通的。我想,一种补救方法是使用相对步长,但我的一些输入范围从0到xx变化,所以可能相对步长不是最好的。有没有一种方法可以把FD步长作为向量或者类似的东西;
for out in out
我正在创建一个程序,它有一个circle对象,它在画布上做圆周运动。我已经创建了圆和渲染类,现在正在进行移动。到目前为止,这是我的代码,不幸的是它不能工作。我不确定从这里到哪里去,或者我是否在正确的轨道上。
import java.util.*;
public class Rotation
{
public static void main(String[] args)
{
int x;
int y;
int radius = 30;
double step = .01;
Scanner input = new Scanner(System.
在bash中,我尝试将PPID导出为环境变量。但是,在我看来,结果在不同的Linux发行版中是非常不一致的。
在Ubuntu18.04上,下面的大多数命令都是“按预期”工作的。在NixOS上,我得到以下行为:
# normally PPID is not part of the environment
env | grep PPID
# --> no output
# exporting for a single command does not work
PPID=foo some-command
# --> bash: PPID: readonly variable
# e
我试图通过使用一个结构数组来创建MIDI内容,用于顺序播放。每个结构创建一个“度量值”,每个度量都具有location、numberOfBeats和tempo的属性。我想循环每一个节拍内的每一个度量,并播放两个MIDI音符中的一个,根据这个节奏是否是测量中的第一个音符。
下面的代码将正确打印(click1,单击2,click2,click1,单击2,click1,click2,click2),但是MIDI内容(附在wav上)将只播放一次,而不会进入数组中的下一个度量值。我一直试图让这个功能按照顺序使用不同的循环类型来完成一系列的度量,但是不能让它正确地前进。谢谢!
import UIKit
i
我正在尝试通过解一个ODE来用Scilab做一些电路分析。但我需要根据函数的当前值更改ODE。我已经使用RK4方法在Scala中实现了这个解决方案,它工作得很好。现在,我正在尝试使用Scilab中的标准函数来做同样的事情。而且它也不起作用。我试着分别解决这两个常微分方程,这是可以的。
clear
state = 0 // state 0 is charging, 1 is discharging
vb = 300.0; vt = 500.0;
r = 100.0; rd = 10.0;
vcc = 600;
c = 48.0e-6;
function dudx = curfunc(t, u
我在模拟指数衰减。为什么它不工作?下面是代码
from math import *
from random import random
from time import sleep
class atom():
def __init__(self, pos):
self.pos = pos
self.isalive = True
class substance():
def __init__(self, halflife):
self.halflife = float(halflife)
self.rate =