1 SUBROUTINE set_dual(data1, hi, yi, yi2, hp, yp, yp2, wten, alsq,
4 USE vmec_input,
ONLY: tensi, tensi2, tensp, fpolyi
5 USE vparams,
ONLY: zero
10 INTEGER :: niota, npres, nots, info
11 REAL(rprec),
DIMENSION(nots) :: data1
12 REAL(rprec),
DIMENSION(niota) :: hi, yi, yi2
13 REAL(rprec),
DIMENSION(npres) :: hp, yp, yp2
14 REAL(rprec),
DIMENSION(nots) :: wten
15 REAL(rprec),
DIMENSION(nots,nots) :: alsq
26 CALL add_tension (alsq, wten, hi, tensi, tensi2, fpolyi, niota, nb
28 CALL add_tension (alsq, wten, hp, tensp, zero, zero, npres, nb,
35 CALL solver (alsq, data1, nots, 1, info)
36 yi(:niota) = data1(:niota)
37 yp(:npres) = data1(1+niota:npres+niota)
42 CALL gety2 (yi, yi2, hi, niota, nb)
43 CALL gety2 (yp, yp2, hp, npres, nb)
45 END SUBROUTINE set_dual