Complex numbers for Dahlia

complex.fuse

record complex {
    real: float;
    imag: float
}

def complex_mul(a: complex, b: complex): complex = {
    let cr: float = (a.real * b.real) - (a.imag * b.imag);
    let ci: float = (a.real * b.imag) + (a.imag * b.real);
    let c: complex = { real=cr; imag=ci };
    return c;
}

def complex_add(a: complex, b: complex): complex = {
    let cr: float = a.real + b.real;
    let ci: float = a.imag + b.imag;
    let c: complex = { real=cr; imag=ci };
    return c;
}

def complex_fma(a: complex, b: complex, c: complex): complex = {
    let cr: float = c.real + (a.real * b.real) - (a.imag * b.imag);
    let ci: float = c.imag + (a.real * b.imag) + (a.imag * b.real);
    let rc: complex = { real=cr; imag=ci };
    return rc;
}

def complex_scale(s: float, c: complex): complex = {
  let sr: float = c.real * s;
  let si: float = c.imag * s;
  let sc: complex = { real=sr; imag=si };
  return sc;
}