Rory Conlin
Rory Conlin
- [ ] near axis constraints with omega - [ ] constant offset surface - [ ] Proximal projection - [ ] Add stuff for spectral condensation / optimal fitting...
```python for n in range(3,13): basis = ZernikePolynomial(n,n) grid1 = QuadratureGrid(n,n,0) grid2 = ConcentricGrid(n,n,0) print(f"{n} {grid1.num_nodes/basis.num_modes:.2f} {grid2.num_nodes/basis.num_modes:.2f}") ``` ``` 3 2.80 1.00 4 3.00 1.00 5 3.14 1.00 6 3.25...
Right now we still have cases where the proximal optimizer keeps trying smaller and smaller steps until eventually it quits saying ``Warning: A bad approximation caused failure to predict improvement.``...
Adds functionality for tracing particles in an arbitrary `MagneticField`, using either guiding center or full orbit equations.
There are a bunch of warnings that can get thrown by `InputReader.parse_vmec_input`, but if we're doing `FourierRZToroidalSurface.from_input_file` or similar, most of them don't matter and may seem confusing.
- [x] #901 - [ ] Shape Gradient ([link](https://arxiv.org/pdf/1910.14144.pdf), [another link](https://iopscience.iop.org/article/10.1088/1741-4326/aac1c7), [one more link](https://arxiv.org/pdf/1910.14144.pdf), [final one](https://www.cambridge.org/core/journals/journal-of-plasma-physics/article/computing-the-shape-gradient-of-stellarator-coil-complexity-with-respect-to-the-plasma-boundary/8A0ED88DD11BEA08C92B45B175A50DDE)) - [x] #895 - [x] #896 - [x] #897 - [x] #898 - [...
```python eq = desc.examples.get("W7-X") eq.change_resolution(3,3,3,6,6,6) assert eq.M == 3 assert eq.surface.M == 3 R1 = eq.surface.R_lmn # these are coeffs from truncating boundary representation to M,N=3 R2 = eq.get_surface_at(rho=1.0).R_lmn #...
nstead of doing something like ```python R_modes = np.vstack( ( [0, 0, 0], eq.surface.R_basis.modes[ np.max(np.abs(eq.surface.R_basis.modes), 1) > k, : ], ) ) Z_modes = eq.surface.Z_basis.modes[ np.max(np.abs(eq.surface.Z_basis.modes), 1) > k, :...
# Basics: - each coil/curve has its own parameterization, knows how to compute basic stuff (position along coil etc) - base class takes care of more general derived quantities -...
Nocedal & Wright give a method for log-barrier interior point methods on pp. 576, should be straightforward to adapt to the Coleman-Li interior point method we use.