Denis St-Onge
Denis St-Onge
The Barcelona team mentioned that the redistribution routines uses for the mirror term might be better performed using the MPI_alltoall library functions, rather than the routines in redistribute.f90. I don't...
As it stands, dissipation.f90 is clocking in at 6182 lines of code, which means it's probably too big. I've already removed the sources/sinks from the file, and perhaps I can...
This is a similar issue to #59, but seems to be much stronger and present even when zed_upwind = 0. When doing the standard RH test case with adiabatic electrons...
Currently, `stella` uses an LU decomposition to perform the linear solve of the response matrix for the parallel streaming term. This isn't the only option, however; GS2 uses a full...
Currently, `stella` uses a sequential approach to operator splitting, where an operator partially advances $g_s$ and passes off to the next operator in the sequence. This isn't the only way...
Currently stella's scaleability is limited by its parallelisation of velocity space alone, relying on a redistribution between a space-local grid and a velocity-local grid. I think the fastest and most...
Recently, adiabatic_option has been moved from dist_fn_knobs to zgrid_parameters. This makes old input files obsolete, so we may want to either throw a deprecation notice/error if dist_fn_knobs is found, or...
Integrating the exact quasineutrality solve with the response matrix approach is a work in progress that will soon be completely. This however, doesn't include any radial variation from terms in...
Currently when the flux-tube train is used in Stella, the k_y = 0 modes of the train are uncoupled, and so these modes are treated as periodic within a single...
The mirror term can be treated fully implicitly, even with radial profile variation. The idea is to inverse Fourier transform to x space then multiply by the appropriate factor give...