Class: Complex

Inherits:
Numeric show all
Defined in:
opal/opal/corelib/complex.rb

Constant Summary collapse

I = 
new(0, 1)

Instance Attribute Summary collapse

Class Method Summary collapse

Instance Method Summary collapse

Methods inherited from Numeric

#%, #+@, #<=>, #__coerced__, #ceil, #clone, #div, #divmod, #dup, #floor, #i, #integer?, #negative?, #nonzero?, #positive?, #round, #step, #to_c, #to_int, #truncate, #zero?

Methods included from Comparable

#<, #<=, #>, #>=, #between?, #clamp

Constructor Details

#initialize(real, imag = 0) ⇒ Complex

Returns a new instance of Complex.



27
28
29
30
 
# File 'opal/opal/corelib/complex.rb', line 27

def initialize(real, imag = 0)
  @real = real
  @imag = imag
end

Instance Attribute Details

#imag ⇒ Object (readonly) Also known as: imaginary

Returns the value of attribute imag.



25
26
27
 
# File 'opal/opal/corelib/complex.rb', line 25

def imag
  @imag
end

#real ⇒ Object (readonly)

Returns the value of attribute real.



25
26
27
 
# File 'opal/opal/corelib/complex.rb', line 25

def real
  @real
end

Class Method Details

.from_string(str) ⇒ Object 



296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
 
# File 'opal/opal/corelib/complex.rb', line 296

def self.from_string(str)
  %x{
    var re = /[+-]?[\d_]+(\.[\d_]+)?(e\d+)?/,
        match = str.match(re),
        real, imag, denominator;

    function isFloat() {
      return re.test(str);
    }

    function cutFloat() {
      var match = str.match(re);
      var number = match[0];
      str = str.slice(number.length);
      return number.replace(/_/g, '');
    }

    // handles both floats and rationals
    function cutNumber() {
      if (isFloat()) {
        var numerator = parseFloat(cutFloat());

        if (str[0] === '/') {
          // rational real part
          str = str.slice(1);

          if (isFloat()) {
            var denominator = parseFloat(cutFloat());
            return #{Rational(`numerator`, `denominator`)};
          } else {
            // reverting '/'
            str = '/' + str;
            return numerator;
          }
        } else {
          // float real part, no denominator
          return numerator;
        }
      } else {
        return null;
      }
    }

    real = cutNumber();

    if (!real) {
      if (str[0] === 'i') {
        // i => Complex(0, 1)
        return #{Complex(0, 1)};
      }
      if (str[0] === '-' && str[1] === 'i') {
        // -i => Complex(0, -1)
        return #{Complex(0, -1)};
      }
      if (str[0] === '+' && str[1] === 'i') {
        // +i => Complex(0, 1)
        return #{Complex(0, 1)};
      }
      // anything => Complex(0, 0)
      return #{Complex(0, 0)};
    }

    imag = cutNumber();
    if (!imag) {
      if (str[0] === 'i') {
        // 3i => Complex(0, 3)
        return #{Complex(0, `real`)};
      } else {
        // 3 => Complex(3, 0)
        return #{Complex(`real`, 0)};
      }
    } else {
      // 3+2i => Complex(3, 2)
      return #{Complex(`real`, `imag`)};
    }
  }
end

.polar(r, theta = 0) ⇒ Object 



17
18
19
20
21
22
23
 
# File 'opal/opal/corelib/complex.rb', line 17

def self.polar(r, theta = 0)
  unless Numeric === r && r.real? && Numeric === theta && theta.real?
    raise TypeError, 'not a real'
  end

  new(r * Math.cos(theta), r * Math.sin(theta))
end

.rect(real, imag = 0) ⇒ Object 



5
6
7
8
9
10
11
 
# File 'opal/opal/corelib/complex.rb', line 5

def self.rect(real, imag = 0)
  unless Numeric === real && real.real? && Numeric === imag && imag.real?
    raise TypeError, 'not a real'
  end

  new(real, imag)
end

.rectangular ⇒ Object 



14
15
16
17
18
19
20
 
# File 'opal/opal/corelib/complex.rb', line 14

def self.rect(real, imag = 0)
  unless Numeric === real && real.real? && Numeric === imag && imag.real?
    raise TypeError, 'not a real'
  end

  new(real, imag)
end

Instance Method Details

#*(other) ⇒ Object 



76
77
78
79
80
81
82
83
84
85
86
 
# File 'opal/opal/corelib/complex.rb', line 76

def *(other)
  if Complex === other
    Complex(@real * other.real - @imag * other.imag,
      @real * other.imag + @imag * other.real,
    )
  elsif Numeric === other && other.real?
    Complex(@real * other, @imag * other)
  else
    __coerced__ :*, other
  end
end

#**(other) ⇒ Object 



103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
 
# File 'opal/opal/corelib/complex.rb', line 103

def **(other)
  if other == 0
    return Complex.new(1, 0)
  end

  if Complex === other
    r, theta = polar
    ore      = other.real
    oim      = other.imag
    nr       = Math.exp(ore * Math.log(r) - oim * theta)
    ntheta   = theta * ore + oim * Math.log(r)

    Complex.polar(nr, ntheta)
  elsif Integer === other
    if other > 0
      x = self
      z = x
      n = other - 1

      while n != 0
        div, mod = n.divmod(2)
        while mod == 0
          x = Complex(x.real * x.real - x.imag * x.imag, 2 * x.real * x.imag)
          n = div
          div, mod = n.divmod(2)
        end

        z *= x
        n -= 1
      end

      z
    else
      (Rational.new(1, 1) / self)**-other
    end
  elsif Float === other || Rational === other
    r, theta = polar

    Complex.polar(r**other, theta * other)
  else
    __coerced__ :**, other
  end
end

#+(other) ⇒ Object 



56
57
58
59
60
61
62
63
64
 
# File 'opal/opal/corelib/complex.rb', line 56

def +(other)
  if Complex === other
    Complex(@real + other.real, @imag + other.imag)
  elsif Numeric === other && other.real?
    Complex(@real + other, @imag)
  else
    __coerced__ :+, other
  end
end

#-(other) ⇒ Object 



66
67
68
69
70
71
72
73
74
 
# File 'opal/opal/corelib/complex.rb', line 66

def -(other)
  if Complex === other
    Complex(@real - other.real, @imag - other.imag)
  elsif Numeric === other && other.real?
    Complex(@real - other, @imag)
  else
    __coerced__ :-, other
  end
end

#-@ ⇒ Object 



52
53
54
 
# File 'opal/opal/corelib/complex.rb', line 52

def -@
  Complex(-@real, -@imag)
end

#/(other) ⇒ Object Also known as: divide, quo 



88
89
90
91
92
93
94
95
96
97
98
99
100
101
 
# File 'opal/opal/corelib/complex.rb', line 88

def /(other)
  if Complex === other
    if (Number === @real && @real.nan?) || (Number === @imag && @imag.nan?) ||
       (Number === other.real && other.real.nan?) || (Number === other.imag && other.imag.nan?)
      Complex.new(Float::NAN, Float::NAN)
    else
      self * other.conj / other.abs2
    end
  elsif Numeric === other && other.real?
    Complex(@real.quo(other), @imag.quo(other))
  else
    __coerced__ :/, other
  end
end

#==(other) ⇒ Object 



42
43
44
45
46
47
48
49
50
 
# File 'opal/opal/corelib/complex.rb', line 42

def ==(other)
  if Complex === other
    @real == other.real && @imag == other.imag
  elsif Numeric === other && other.real?
    @real == other && @imag == 0
  else
    other == self
  end
end

#abs ⇒ Object Also known as: magnitude 



147
148
149
 
# File 'opal/opal/corelib/complex.rb', line 147

def abs
  Math.hypot(@real, @imag)
end

#abs2 ⇒ Object 



151
152
153
 
# File 'opal/opal/corelib/complex.rb', line 151

def abs2
  @real * @real + @imag * @imag
end

#angle ⇒ Object Also known as: arg 



155
156
157
 
# File 'opal/opal/corelib/complex.rb', line 155

def angle
  Math.atan2(@imag, @real)
end

#coerce(other) ⇒ Object 



32
33
34
35
36
37
38
39
40
 
# File 'opal/opal/corelib/complex.rb', line 32

def coerce(other)
  if Complex === other
    [other, self]
  elsif Numeric === other && other.real?
    [Complex.new(other, 0), self]
  else
    raise TypeError, "#{other.class} can't be coerced into Complex"
  end
end

#conj ⇒ Object Also known as: conjugate 



161
162
163
 
# File 'opal/opal/corelib/complex.rb', line 161

def conj
  Complex(@real, -@imag)
end

#denominator ⇒ Object 



167
168
169
 
# File 'opal/opal/corelib/complex.rb', line 167

def denominator
  @real.denominator.lcm(@imag.denominator)
end

#eql?(other) ⇒ Boolean

Returns:



173
174
175
 
# File 'opal/opal/corelib/complex.rb', line 173

def eql?(other)
  Complex === other && @real.class == @imag.class && self == other
end

#fdiv(other) ⇒ Object 



177
178
179
180
181
182
183
 
# File 'opal/opal/corelib/complex.rb', line 177

def fdiv(other)
  unless Numeric === other
    raise TypeError, "#{other.class} can't be coerced into Complex"
  end

  self / other
end

#finite? ⇒ Boolean

Returns:



185
186
187
 
# File 'opal/opal/corelib/complex.rb', line 185

def finite?
  @real.finite? && @imag.finite?
end

#hash ⇒ Object 



189
190
191
 
# File 'opal/opal/corelib/complex.rb', line 189

def hash
  "Complex:#{@real}:#{@imag}"
end

#infinite? ⇒ Boolean

Returns:



195
196
197
 
# File 'opal/opal/corelib/complex.rb', line 195

def infinite?
  @real.infinite? || @imag.infinite?
end

#inspect ⇒ Object 



199
200
201
 
# File 'opal/opal/corelib/complex.rb', line 199

def inspect
  "(#{self})"
end

#numerator ⇒ Object 



207
208
209
210
211
212
213
 
# File 'opal/opal/corelib/complex.rb', line 207

def numerator
  d = denominator

  Complex(@real.numerator * (d / @real.denominator),
    @imag.numerator * (d / @imag.denominator),
  )
end

#polar ⇒ Object 



217
218
219
 
# File 'opal/opal/corelib/complex.rb', line 217

def polar
  [abs, arg]
end

#rationalize(eps = undefined) ⇒ Object 



225
226
227
228
229
230
231
232
233
234
235
236
237
 
# File 'opal/opal/corelib/complex.rb', line 225

def rationalize(eps = undefined)
  %x{
    if (arguments.length > 1) {
      #{raise ArgumentError, "wrong number of arguments (#{`arguments.length`} for 0..1)"};
    }
  }

  if @imag != 0
    raise RangeError, "can't' convert #{self} into Rational"
  end

  real.rationalize(eps)
end

#real? ⇒ Boolean

Returns:



239
240
241
 
# File 'opal/opal/corelib/complex.rb', line 239

def real?
  false
end

#rect ⇒ Object Also known as: rectangular 



243
244
245
 
# File 'opal/opal/corelib/complex.rb', line 243

def rect
  [@real, @imag]
end

#to_f ⇒ Object 



251
252
253
254
255
256
257
 
# File 'opal/opal/corelib/complex.rb', line 251

def to_f
  unless @imag == 0
    raise RangeError, "can't convert #{self} into Float"
  end

  @real.to_f
end

#to_i ⇒ Object 



259
260
261
262
263
264
265
 
# File 'opal/opal/corelib/complex.rb', line 259

def to_i
  unless @imag == 0
    raise RangeError, "can't convert #{self} into Integer"
  end

  @real.to_i
end

#to_r ⇒ Object 



267
268
269
270
271
272
273
 
# File 'opal/opal/corelib/complex.rb', line 267

def to_r
  unless @imag == 0
    raise RangeError, "can't convert #{self} into Rational"
  end

  @real.to_r
end

#to_s ⇒ Object 



275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
 
# File 'opal/opal/corelib/complex.rb', line 275

def to_s
  result = @real.inspect

  result +=
    if (Number === @imag && @imag.nan?) || @imag.positive? || @imag.zero?
      '+'
    else
      '-'
    end

  result += @imag.abs.inspect

  if Number === @imag && (@imag.nan? || @imag.infinite?)
    result += '*'
  end

  result + 'i'
end