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)
- Documentation fails to build starting at #dc54a2b HOT 2
- "A particle is lost" recursive loop
- Numba does not reproduce Python results HOT 1
- Setup.py file has mcdc version as 0.1.0
- Update installation script to match pyproject requirements
- `Seg fault` with `*** Process received signal ***` HOT 4
- GPU compilation fails after PR #198 (next VR) HOT 1
- Updated Documentation (next VR) HOT 1
- Updating examples (next VR)
- Resolving issues related to iQMC HOT 4
- Citation box missing information on web-hosted docpages HOT 4
- Minor warning in a function documentation
- Visulizer import deprection
- Regression tests for continuous energy physics HOT 1
- Large files are accidentally Git-tracked? HOT 7
- Conda environment management deprecation
- Early Memory Deallocation/Re-Use in Local Arrays HOT 3
- GPU regression testing and reproducibility HOT 3
- Regression test stdev discrepancy when DD is active HOT 1
- Need to compile Numba-CPU first to be able to compile Numba-GPU per #229 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.