add constant summand for newton pool

Author: cl1ckname

src/components/pool/Pool.tsx

@@ -13,6 +13,8 @@ export interface PoolParams {
     max_its: number,
     rea: number,
     ima: number,
+    rec: number,
+    imc: number,
 }
 export default function Pool() {
     const defaultParams: PoolParams = {
@@ -26,6 +28,8 @@ export default function Pool() {
         roots: [0,4*Math.PI / 3, 2*Math.PI /3,0,1,0,1,1,1,1,1,1],
         rea: 1,
         ima: 0,
+        rec: 0,
+        imc: 0,
     }
     const [params, setParams] = useState(defaultParams)
 

src/components/pool/PoolCanvas.tsx

@@ -159,6 +159,13 @@ function putSettingsUniforms(gl: WebGLRenderingContext, program: WebGLProgram, s
     }
     gl.uniform2f(aId, settings.rea, settings.ima)
 
+    const cId = gl.getUniformLocation(program, "point_c")
+    if (!cId) {
+        console.error("c not found")
+        return;
+    }
+    gl.uniform2f(cId, settings.rec, settings.imc)
+
     const resolutionId = gl.getUniformLocation(program, "resol")
     if (!resolutionId) {
         console.error("resolution not found")

src/components/pool/PoolSettings.tsx

@@ -39,5 +39,16 @@ export default function PoolSettings(props: PoolSettingsProps) {
             settingsCopy.ima = n
             props.onChange(settingsCopy)
         }}/>
+        <label>Re(c)</label>
+        <Slider min={0} max={5} value={settingsCopy.rec} step={0.05} onChange={n => {
+            settingsCopy.rec = n
+            props.onChange(settingsCopy)
+        }}/>
+
+        <label>Im(c)</label>
+        <Slider min={0} max={5} value={settingsCopy.imc} step={0.05} onChange={n => {
+            settingsCopy.imc = n
+            props.onChange(settingsCopy)
+        }}/>
     </form>
 }

src/logic/pool/Shader.tsx

@@ -19,6 +19,7 @@ uniform vec2 position;
 uniform vec3 colors[12];
 uniform int max_its;
 uniform vec2 point_a;
+uniform vec2 point_c;
 
 const int MAX_ITS = 300;
 
@@ -58,7 +59,7 @@ vec2 pos = ((gl_FragCoord.xy + position / escale) / resol.xy + vec2(-0.5)) * esc
 vec2 zi = pos;
 for (int i = 0; i < MAX_ITS; i++) {
 	if (i == max_its) {break;}
-	zi = zi - cmul(cdiv(cpow(zi, n) - vec2(100., 0.), float(n) * cpow(zi, n-1)), point_a) + point_a;
+	zi = zi - cmul(cdiv(cpow(zi, n) - vec2(100., 0.), float(n) * cpow(zi, n-1)), point_a) + point_c;
 	// zi = zi - cmul(cdiv(sin(zi) - vec2(1., 0.), cos(zi)), point_a);
 }