from NeuroCore.Models.Conventions import Domains, Steps
from NeuroCore.Approximations.Space.FEM.GalerkinApproximation import GalerkinApproximation
from NeuroCore.Approximations.Time.BackwardEuler import BackwardEuler
from NeuroCore.Models.Conditions.BCs import BoundaryConditions
from NeuroCore.Models.Conditions.Dirichlets import KilledEnd
from Analytics.Solutions.Validations.Transient import ValidationWithF
from Plotting.Simulation import Simulation
from Plotting.IDataPlot import IDataPlot
from numpy import arange
########################################
### Simulation Variable ###
########################################
simDomain = Domains(space=[0, 1], time=[0, 10])
simSteps = Steps(space=0.1, time=0.1)
diffusionValue = 1
reactionValue = 1
kValue = 1
boundaryConditions = BoundaryConditions(KilledEnd(), KilledEnd())
sElements = arange(0, simDomain.space[-1] + simSteps.space,\
simSteps.space)
########################################
### Setting the Simulation ###
########################################
import numpy
import math
import time
########################################
### Simulation Variable ###
########################################
length = 1
T = 10
delta_x = 0.1
delta_t = 0.1
simDomain = Domains(time=[0,T],space=[0,length])
simSteps = Steps(time=delta_t,space=delta_x)
spatial_elements = numpy.arange(0, simDomain.space[-1] + simSteps.space,\
simSteps.space)
time_elements = numpy.arange(0, simDomain.time[-1] + simSteps.time, simSteps.time)
V_0 = []
def u0(x):
return math.sin(2*math.pi*x)
for x in spatial_elements:
V_0.append(u0(x))
from NeuroCore.Approximations.Space.FEM.GalerkinApproximation import GalerkinApproximation
from NeuroCore.Models.Conditions.BCs import BoundaryConditions
from NeuroCore.Models.Conditions.Neumanns import SealedEnd
from NeuroCore.Models.Conventions import Domains, Steps
from NeuroCore.Models.GeneralModel import GeneralModel
from NeuroCore.Neuron.Segment.Base import ISegment
from Plotting.IDataPlot import IDataPlot
from Plotting.IDataPlots import IDataPlots
from Plotting.Simulation import Simulation
########################################
### Simulation Variable ###
########################################
simDomain = Domains(space=[0, 12])
simSteps = Steps(space=0.01)
diffusionValue = 100
reactionValue = 10
boundaryConditions = BoundaryConditions(SealedEnd(), SealedEnd())
sElements = numpy.arange(0, simDomain.space[-1] + simSteps.space, \
simSteps.space)
########################################
### Setting the Simulation ###
########################################
analytical = ValidationWithF(domain=simDomain, steps=simSteps, \
BCs=boundaryConditions, \
from NeuroCore.Approximations.Space.FEM.GalerkinApproximation import GalerkinApproximation
#from Approximations.Time.BackwardEuler import BackwardEuler
from NeuroCore.Models.Conditions.BCs import BoundaryConditions
from NeuroCore.Models.Conditions.Dirichlets import KilledEnd
from math import pi, pow, sin
import numpy
import math
import time
########################################
### Simulation Variable ###
########################################
simDomain = Domains(time=[0, 10], space=[0, 1])
simSteps = Steps(time=0.1, space=0.1)
spatial_elements = numpy.arange(0, simDomain.space[-1] + simSteps.space,\
simSteps.space)
time_elements = numpy.arange(0, simDomain.time[-1] + simSteps.time,
simSteps.time)
V_0 = []
def u0(x):
return math.sin(2 * math.pi * x)
for x in spatial_elements:
from NeuroCore.Models.Conditions.BCs import BoundaryConditions
from NeuroCore.Models.Conditions.Neumanns import Neumann
from Analytics.Solutions.Validations.Stationary import ValidationZeroF as analyCableModel
from NeuroCore.Approximations.Space.FEM.Integration.Trapezoidal import Trapezoidal
from Plotting.Simulation import Simulation
from Plotting.IDataPlot import IDataPlot
import numpy
import time
########################################
### Simulation Variable ###
########################################
simDomain = Domains(space=[0, 1])
simSteps = Steps(space=0.1)
########################################
### Setting the Simulation ###
########################################
def f(x):
return 0
diffusionValue = -0.01
reactionValue = 1
boundaryConditions = BoundaryConditions(Neumann(bcValue=1), Neumann(bcValue=1))