Comments (3)
Interestingly, the test is passing when I double the number of cells in each dimension of the tally mesh from 20 to 40. Here is the original input deck
import numpy as np
import mcdc
# =============================================================================
# Set model
# =============================================================================
# A shielding problem based on Problem 2 of [Coper NSE 2001]
# https://ans.tandfonline.com/action/showCitFormats?doi=10.13182/NSE00-34
# Set materials
SigmaT = 5.0
c = 0.8
m_barrier = mcdc.material(capture=np.array([SigmaT]), scatter=np.array([[SigmaT * c]]))
SigmaT = 1.0
m_room = mcdc.material(capture=np.array([SigmaT]), scatter=np.array([[SigmaT * c]]))
# Set surfaces
sx1 = mcdc.surface("plane-x", x=0.0, bc="reflective")
sx2 = mcdc.surface("plane-x", x=2.0)
sx3 = mcdc.surface("plane-x", x=2.4)
sx4 = mcdc.surface("plane-x", x=4.0, bc="vacuum")
sy1 = mcdc.surface("plane-y", y=0.0, bc="reflective")
sy2 = mcdc.surface("plane-y", y=2.0)
sy3 = mcdc.surface("plane-y", y=4.0, bc="vacuum")
# Set cells
mcdc.cell([+sx1, -sx2, +sy1, -sy2], m_room)
mcdc.cell([+sx1, -sx4, +sy2, -sy3], m_room)
mcdc.cell([+sx3, -sx4, +sy1, -sy2], m_room)
mcdc.cell([+sx2, -sx3, +sy1, -sy2], m_barrier)
# =============================================================================
# iQMC Parameters
# =============================================================================
N = 2e2
Nx = Ny = 20
maxit = 2
tol = 1e-3
x = np.linspace(0, 4, num=Nx + 1)
y = np.linspace(0, 4, num=Ny + 1)
generator = "halton"
# fixed source in lower left corner
fixed_source = np.zeros((Nx, Ny))
fixed_source[0 : int(0.25 * Nx), 0 : int(0.25 * Nx)] = 1
phi0 = np.ones((Nx, Ny))
mcdc.iQMC(
x=x,
y=y,
fixed_source=fixed_source,
phi0=phi0,
maxitt=maxit,
tol=tol,
generator=generator,
score=["tilt-x", "tilt-y", "tilt-xy"],
)
# =============================================================================
# Set tally, setting, and run mcdc
# =============================================================================
# Setting
mcdc.setting(N_particle=N)
# Run
mcdc.run()
from mcdc.
Could we change the test to do that? I know that's an unsatisfying result but this would allow us to move on
from mcdc.
Sorry forgot to follow up on this, but that's what I did. The test is passing now.
from mcdc.
Related Issues (20)
- mesh_get_index returns a value of -1 when a particle lies on mesh minimum boundary HOT 1
- Function Timers HOT 1
- Documentation website initilization HOT 1
- MCDC (iQMC) fails with large processor count
- Updating Regression Tests and Examples HOT 6
- "Particle Lost" Monte Carlo Eigenvalue Bug HOT 5
- SEED_SPLIT_UQ OverflowError on PC HOT 5
- Visualization merging HOT 4
- Issue on Quartz: HYDU_sock_connect unable to connect HOT 2
- Performance testing via `chron` jobs HOT 3
- Visualizer crashes on 3D C5G7 and SMRG7
- adding `cache`ing functionality
- adding CSG operators
- Non-uniform structors of arrays in numba HOT 2
- adding functinality for cell detectors
- `kernel.split_particle` returns by reference in Numba
- `cache` issues HOT 3
- merging `harmonize-integration` with `main`
- Outside GPU libraries for GPU functions HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mcdc.