|
21 | 21 | using PythonPlot: pyplot |
22 | 22 | using LaTeXStrings |
23 | 23 | using Colors |
24 | | - using AugmentedPoissonBoltzmann.ICMPBP: ICMPBP, ICMPBData, AppliedPotentialHalfCell, set_molarity!, calc_cmol, calc_c0mol, calc_χ, get_E, get_φ, get_p, get_c0, |
25 | | - set_κ!, set_φ! |
| 24 | + using AugmentedPoissonBoltzmann.SolverCore: AugmentedPBP, AugmentedPBData, AppliedPotentialHalfCell, set_molarity!, calc_cmol, calc_c0mol, calc_χ, get_E, get_φ, get_p, |
| 25 | + set_κ!, set_φ!, dlcapsweep |
26 | 26 | using DrWatson, PoissonBoltzmannIPAM2025 |
27 | 27 |
|
28 | 28 | end |
@@ -61,7 +61,7 @@ All values are given with respect to SI basic units (m, kg, s, V, A) |
61 | 61 | """ |
62 | 62 |
|
63 | 63 | # ╔═╡ b24b7e23-61ea-41fc-a345-286e904c042b |
64 | | -datavhalf = ICMPBData(χvar = true) |
| 64 | +datavhalf = AugmentedPBData(χvar = true) |
65 | 65 |
|
66 | 66 | # ╔═╡ 1bb47749-edde-4bee-be9f-059a7652b354 |
67 | 67 | begin |
@@ -101,11 +101,11 @@ p1 = let |
101 | 101 | M = molarities[i] |
102 | 102 | set_molarity!(halfcell, M) |
103 | 103 |
|
104 | | - volts, dlcaps = ICMPBP.dlcapsweep(halfcell; φ_max) |
| 104 | + volts, dlcaps = dlcapsweep(halfcell; φ_max) |
105 | 105 | plot(volts, dlcaps / (μF / cm^2), color = colors[i], linestyle = "--") |
106 | 106 |
|
107 | 107 | set_molarity!(halfcelldd, M) |
108 | | - volts, dlcaps = ICMPBP.dlcapsweep(halfcelldd; φ_max, damp_initial = 0.5) |
| 108 | + volts, dlcaps = dlcapsweep(halfcelldd; φ_max, damp_initial = 0.5) |
109 | 109 | plot(volts, dlcaps / (μF / cm^2), label = "M=$(M)", color = colors[i]) |
110 | 110 | end |
111 | 111 | ax.set_xlabel("U/V") |
@@ -140,10 +140,10 @@ p2 = let |
140 | 140 | set_κ!(halfcell, κ) |
141 | 141 | set_κ!(halfcelldd, κ) |
142 | 142 |
|
143 | | - volts, dlcaps = ICMPBP.dlcapsweep(halfcelldd; φ_max, damp_initial = 0.5) |
| 143 | + volts, dlcaps = dlcapsweep(halfcelldd; φ_max, damp_initial = 0.5) |
144 | 144 | plot(volts, dlcaps / (μF / cm^2), label = "κ=$(κ)", color = colors[i]) |
145 | 145 |
|
146 | | - volts, dlcaps = ICMPBP.dlcapsweep(halfcell; φ_max) |
| 146 | + volts, dlcaps = dlcapsweep(halfcell; φ_max) |
147 | 147 | plot(volts, dlcaps / (μF / cm^2), color = colors[i], linestyle = "--") |
148 | 148 |
|
149 | 149 | end |
|
0 commit comments