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.



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

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

Instance Attribute Details

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

Returns the value of attribute imag.



23
24
25
 
# File 'opal/opal/corelib/complex.rb', line 23

def imag
  @imag
end

#real ⇒ Object (readonly)

Returns the value of attribute real.



23
24
25
 
# File 'opal/opal/corelib/complex.rb', line 23

def real
  @real
end

Class Method Details

.from_string(str) ⇒ Object 



273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
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
 
# File 'opal/opal/corelib/complex.rb', line 273

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 #{::Kernel.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 #{::Kernel.Complex(0, 1)};
      }
      if (str[0] === '-' && str[1] === 'i') {
        // -i => Complex(0, -1)
        return #{::Kernel.Complex(0, -1)};
      }
      if (str[0] === '+' && str[1] === 'i') {
        // +i => Complex(0, 1)
        return #{::Kernel.Complex(0, 1)};
      }
      // anything => Complex(0, 0)
      return #{::Kernel.Complex(0, 0)};
    }

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

.polar(r, theta = 0) ⇒ Object 



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

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

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

.rect(real, imag = 0) ⇒ Object Also known as: rectangular 



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

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

  new(real, imag)
end

Instance Method Details

#*(other) ⇒ Object 



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

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

#**(other) ⇒ Object 



102
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
 
# File 'opal/opal/corelib/complex.rb', line 102

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 = ::Kernel.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 



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

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

#-(other) ⇒ Object 



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

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

#-@ ⇒ Object 



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

def -@
  ::Kernel.Complex(-@real, -@imag)
end

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



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

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?
    ::Kernel.Complex(@real.quo(other), @imag.quo(other))
  else
    __coerced__ :/, other
  end
end

#==(other) ⇒ Object 



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

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 



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

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

#abs2 ⇒ Object 



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

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

#angle ⇒ Object Also known as: arg 



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

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

#coerce(other) ⇒ Object 



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

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

#conj ⇒ Object Also known as: conjugate 



158
159
160
 
# File 'opal/opal/corelib/complex.rb', line 158

def conj
  ::Kernel.Complex(@real, -@imag)
end

#denominator ⇒ Object 



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

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

#eql?(other) ⇒ Boolean

Returns:



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

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

#fdiv(other) ⇒ Object 



170
171
172
173
174
175
176
 
# File 'opal/opal/corelib/complex.rb', line 170

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

  self / other
end

#finite? ⇒ Boolean

Returns:



178
179
180
 
# File 'opal/opal/corelib/complex.rb', line 178

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

#hash ⇒ Object 



182
183
184
 
# File 'opal/opal/corelib/complex.rb', line 182

def hash
  [::Complex, @real, @imag].hash
end

#infinite? ⇒ Boolean

Returns:



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

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

#inspect ⇒ Object 



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

def inspect
  "(#{self})"
end

#numerator ⇒ Object 



194
195
196
197
198
199
200
 
# File 'opal/opal/corelib/complex.rb', line 194

def numerator
  d = denominator

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

#polar ⇒ Object 



202
203
204
 
# File 'opal/opal/corelib/complex.rb', line 202

def polar
  [abs, arg]
end

#rationalize(eps = undefined) ⇒ Object 



206
207
208
209
210
211
212
213
214
215
216
217
218
 
# File 'opal/opal/corelib/complex.rb', line 206

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

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

  real.rationalize(eps)
end

#real? ⇒ Boolean

Returns:



220
221
222
 
# File 'opal/opal/corelib/complex.rb', line 220

def real?
  false
end

#rect ⇒ Object Also known as: rectangular 



224
225
226
 
# File 'opal/opal/corelib/complex.rb', line 224

def rect
  [@real, @imag]
end

#to_f ⇒ Object 



228
229
230
231
232
233
234
 
# File 'opal/opal/corelib/complex.rb', line 228

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

  @real.to_f
end

#to_i ⇒ Object 



236
237
238
239
240
241
242
 
# File 'opal/opal/corelib/complex.rb', line 236

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

  @real.to_i
end

#to_r ⇒ Object 



244
245
246
247
248
249
250
 
# File 'opal/opal/corelib/complex.rb', line 244

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

  @real.to_r
end

#to_s ⇒ Object 



252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
 
# File 'opal/opal/corelib/complex.rb', line 252

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