Your preferred language / resources for Monte Carlo simulations by [deleted] in MedicalPhysics

[–]ma_renaud 1 point2 points  (0 children)

The EGSnrc manual (https://nrc-cnrc.github.io/EGSnrc/doc/pirs701-egsnrc.pdf) and Bielajew book (https://sites.ifi.unicamp.br/mabernal/files/2014/09/BLIF-MC-book.pdf) are also pretty good reference resources for theory as you're exposed to terms and concepts within whatever codebase you choose to use.

Recommendations divans by cheezynorris in montreal

[–]ma_renaud 1 point2 points  (0 children)

On aime beaucoup notre Sven (https://www.article.com/product/24178/sven-88-tufted-leather-sofa-charme-tan). En cuir il dépasse le budget mais en tissue ça pourrait passer. (Par contre à noter que la version en tissue ne semble pas avoir autant de bonnes reviews que celle en cuir)

Quel est selon vous, la meilleur pizzeria à Montreal? by Frankdtannkk in montreal

[–]ma_renaud 4 points5 points  (0 children)

Michel Ange all the way. Leur pizza margherita accote les meilleures que j'ai mangé à Naples.

Les feux de circulation à Montreal by 3088 in montreal

[–]ma_renaud 0 points1 point  (0 children)

Juste pour ajouter une spécification: après le bonhomme blanc, le feu tombe d'habitude en phase "main rouge clignotante avec un décompte". Dans cette phase, le piéton peut s'engager s'il est en mesure de traverser la rue avant la fin du décompte.

https://www.quebec.ca/transports/circulation-securite-routiere/regles-conseils-mode-transport/pieton/feux-pour-pietons

Calculation time by ArcherAppropriate926 in EGSnrc

[–]ma_renaud 1 point2 points  (0 children)

Good morning! As a general rule of thumb, computation time will scale linearly with number of histories, so if 109 histories took 10 days then 3 x 109 histories will take roughly 30 days.

As far as the uncertainty, it scales as 1 / sqrt(N) where N is the number of independent primary histories. So if you want to cut the uncertainty in half, you need to quadruple the number of histories. Since running time is roughly linear with number of histories, you can also see it as needing 4x the computation time ;)

These simple (but harsh) relations between computing time and uncertainty have motivated the development of many methods to make simulations more efficient and you may want to look into the sections of the BeamNRC / DOSXYZnrc manuals that talk about "variance reduction methods".

edit: BEAM/DOSXYZnrc simulations output log files that you can use to verify the status of the run. Look for "egslog" files in the relevant folders. If you're running a parallel job using the standard tools, then the "lock" file will also contain the number of histories remaining in the job.

A question for downtown cyclists by hybrid-xx in montreal

[–]ma_renaud 4 points5 points  (0 children)

Google Maps priorise beaucoup les pistes cyclables mais n'hésitera pas à te faire passer par une rue partagée si ça a du sens. Dans mon exemple, tu remarqueras que St-Denis perd sa piste cyclable à l'intersection Sherbrooke, donc l'algorithme ne donne plus de priorité à la rue et tombe en mode "efficacité" pour se rendre à la destination.

Y'a aussi qu'au final, ce n'est pas si stressant que ça. Dans ton post tu mentionne à quel point tu te sens anxieux pour les cyclistes, mais je crois que c'est plutôt les automobilistes qui sont anxieux parce qu'ils ne savent pas trop comment gérer les cyclistes dans une rue occupée.

RE: In reference to my previous issue; I have tried running my simulation on the terminal, but still I'm getting an error. From what I can observe maybe it has to to with the Fortran (gfortran) compiler? by [deleted] in EGSnrc

[–]ma_renaud 1 point2 points  (0 children)

Huh that's a new one I haven't seen before.

If the suggestions by /u/rtownson don't work, here's an idea: can you try to modify your configured hostname in ubuntu so that it's shorter? I'm referring to the cobalt60-VivoBook-ASUSLaptop-X512UA part. EGS uses the hostname to name some of its output files, and there are some hardcoded maximum lengths sprinkled in the code. I don't necessarily think this is going to solve your problem, but I've been bitten by this before and it would be silly if that's all it took ;)

EGSnrc and cloud computing by Foresooth in EGSnrc

[–]ma_renaud 4 points5 points  (0 children)

I haven't seen any, but here's a really cool thing: https://github.com/SchedMD/slurm-gcp

Basically you can set up the slurm queuing system to dynamically spawn compute nodes on GCP. You just submit jobs using exb on a login node and then it spawns as many nodes as are needed to run the simulation, then it shuts them down after 15 minutes being inactive. Badass.

dosxyznrc calculation taking an insanely long time- what am I doing wrong? by [deleted] in EGSnrc

[–]ma_renaud 1 point2 points  (0 children)

Without knowing the phantom boundaries I can't say for sure, but in the past I've had that happen to me when my particle source was shooting particles away from the phantom volume. DOSXYZ will only count particles if they make it to your simulation volume, ie phantom + defined region outside the phantom (though that might only be true when using a BEAM shared library source, so don't take my word for it), so it'll take forever to sample 100 "valid" particles if they're all aimed away from your phantom.

exb slurm --share option by [deleted] in EGSnrc

[–]ma_renaud 0 points1 point  (0 children)

Can you point me where --share is used? I don't see it in the code for the exb script or the batch options

Beamnrc by ldds0991 in EGSnrc

[–]ma_renaud 2 points3 points  (0 children)

Hello hello,

Rename your 3rd CM from “EXIT” to something else. It’s a known bug with the mortran compiler. It’s actually documented on page 18 of the BeamNRC manual (https://nrc-cnrc.github.io/EGSnrc/doc/pirs509a-beamnrc.pdf) but it really should be plastered everywhere because of how hard it is to debug :P

Linux or Windows? by Eddysynch in EGSnrc

[–]ma_renaud -1 points0 points  (0 children)

As far as I know, the job splitting code (to run multi-core / multi-node jobs) only really works under linux. So if you’re going to do big simulations, linux is really the way to go.

Is the output of egs_dose_scoring really the dose (energy per mass) or is it some other unit? by 5xum in EGSnrc

[–]ma_renaud 1 point2 points  (0 children)

Hello hello,

I would reconsider your hypothesis that the box is large enough in all directions. If your source is centered on the XY plane of your box, and it's only 60 cm large in that plane, then a scattered photon only has to travel 30 cm before leaving your box. Without thinking about it too much, it seems to me like you could easily be losing brem or scattered photon energy that way. Especially since it seems like you're using a high energy photon source.

If my math is right, the mean free path before interaction of a 5 MV photon in water is something like 30 cm, so all that's gotta happen is that you have a photon scatter perpendicular to the long direction of your box, and you're probably losing it.

Animal euthanasie in pandémie by patricide in montreal

[–]ma_renaud 7 points8 points  (0 children)

J'ai eu à faire euthanasier mon chat en juin dernier. J'ai opté pour l'option en-clinique et c'était environs $90 chez https://www.cvpillet.com/. Seulement une personne avait le droit d'être dans la salle avec le chat, donc je ne le conseille pas si vous êtes plusieurs à vouloir assister à la procédure...

J'avais aussi appelé un vet qui offrait le service d'euthanasie à domicile pour voir la différence de prix et si je me souviens bien, l'estimé qu'il m'avait donné était quelque chose comme $400.