Class: NArrayMiss

Inherits:
Object
  • Object
show all
Includes:
NMMath
Defined in:
lib/narray_miss.rb,
lib/narray_miss.rb

Constant Summary collapse

BYTE =

Class Constants

— NArrayMiss::BYTE

type code for 1 byte unsigned integer.

— NArrayMiss::SINT

type code for 2 byte signed integer.

— NArrayMiss::INT

type code for 4 byte signed integer.

— NArrayMiss::SFLOAT

type code for single precision float.

— NArrayMiss::FLOAT

type code for double precision float.

— NArrayMiss::SCOMPLEX

type code for single precision complex.

— NArrayMiss::COMPLEX

type code for double precision complex.

— NArrayMiss::OBJECT

type code for Ruby object.

go back to ((<Index>))

NArray::BYTE
SINT = 
NArray::SINT
INT = 
NArray::INT
SFLOAT = 
NArray::SFLOAT
FLOAT = 
NArray::FLOAT
SCOMPLEX = 
NArray::SCOMPLEX
COMPLEX = 
NArray::COMPLEX
OBJECT = 
NArray::OBJECT

Class Method Summary collapse

Instance Method Summary collapse

Class Method Details

.[](*arg) ⇒ Object 



147
148
149
 
# File 'lib/narray_miss.rb', line 147

def self.[](*arg)
  NArrayMiss.to_nam(NArray[*arg])
end

.__new__ ⇒ Object 



73
 
# File 'lib/narray_miss.rb', line 73

alias :__new__ :new

._load(o) ⇒ Object 



1267
1268
1269
1270
1271
1272
 
# File 'lib/narray_miss.rb', line 1267

def self._load(o)
  ary, mask = Marshal::load(o)
  ary = NArray._load(ary)
  mask = NArray._load(mask)
  NArrayMiss.to_nam_no_dup(ary,mask)
end

.byte(*arg) ⇒ Object 



123
124
125
 
# File 'lib/narray_miss.rb', line 123

def self.byte(*arg)
  NArrayMiss.new(BYTE,*arg)
end

.complex(*arg) ⇒ Object 



141
142
143
 
# File 'lib/narray_miss.rb', line 141

def self.complex(*arg)
  NArrayMiss.new(COMPLEX,*arg)
end

.float(*arg) ⇒ Object 



135
136
137
 
# File 'lib/narray_miss.rb', line 135

def self.float(*arg)
  NArrayMiss.new(FLOAT,*arg)
end

.int(*arg) ⇒ Object 



129
130
131
 
# File 'lib/narray_miss.rb', line 129

def self.int(*arg)
  NArrayMiss.new(INT,*arg)
end

.new(*arg) ⇒ Object

Class Methods

— NArrayMiss.new(typecode, size, …)

create (({NArrayMiss})) of ((|typecode|)).
All elements are initialized with 0.

— NArrayMiss.byte(size, …)

same as NArrayMiss.new(NArrayMiss::BYTE, ((|size|)), ...).

— NArrayMiss.sint(size, …)

same as NArrayMiss.new(NArrayMiss::SINT, ((|size|)), ...).

— NArrayMiss.int(size, …)

same as NArrayMiss.new(NArrayMiss::INT, ((|size|)), ...).

— NArrayMiss.sfloat(size, …)

same as NArrayMiss.new(NArrayMiss::SFLOAT, ((|size|)), ...).

— NArrayMiss.float(size, …)

same as NArrayMiss.new(NArrayMiss::FLOAT, ((|size|)), ...).

— NArrayMiss.scomplex(size, …)

same as NArrayMiss.new(NArrayMiss::SCOMPLEX, ((|size|)), ...).

— NArrayMiss.complex(size, …)

same as NArrayMiss.new(NArrayMiss::COMPLEX, ((|size|)), ...).

— NArrayMiss.object(size, …)

same as NArrayMiss.new(NArrayMiss::OBJECT, ((|size|)), ...).

— NArrayMiss[](value, …)

create (({NArrayMiss})) form [((|value|)), ...].

— NArrayMiss.to_nam(array [,mask])

create (({NArrayMiss})) from ((|array|)).
((|array|)) must be (({Numeric})) (({Array})) or (({NArray})).

— NArrayMiss.to_nam_no_dup(array [,mask])

convert from ((|array|)) to (({NArrayMiss})).

go back to ((<Index>))



118
119
120
121
122
 
# File 'lib/narray_miss.rb', line 118

def self.new(*arg)
  array = NArray.new(*arg)
  mask = NArray.byte(*arg[1..-1])
  __new__(array, mask)
end

.object(*arg) ⇒ Object 



144
145
146
 
# File 'lib/narray_miss.rb', line 144

def self.object(*arg)
  NArrayMiss.new(OBJECT,*arg)
end

.scomplex(*arg) ⇒ Object 



138
139
140
 
# File 'lib/narray_miss.rb', line 138

def self.scomplex(*arg)
  NArrayMiss.new(SCOMPLEX,*arg)
end

.sfloat(*arg) ⇒ Object 



132
133
134
 
# File 'lib/narray_miss.rb', line 132

def self.sfloat(*arg)
  NArrayMiss.new(SFLOAT,*arg)
end

.sint(*arg) ⇒ Object 



126
127
128
 
# File 'lib/narray_miss.rb', line 126

def self.sint(*arg)
  NArrayMiss.new(SINT,*arg)
end

.to_nam(*arg) ⇒ Object 



185
186
187
188
189
190
191
192
193
194
 
# File 'lib/narray_miss.rb', line 185

def self.to_nam(*arg)
  if !(Numeric===arg[0]) && !(Array===arg[0]) && !arg[0].is_a?(NArray)
    raise "first argument must be Numeric, NArray or Array"
  end
  arg[0] = arg[0].dup if !(Numeric===arg[0])
  if arg.length==2 && !(Numeric===arg[1]) && arg[1].class!=TrueClass && arg[1].class!=FalseClass then
    arg[1] = arg[1].dup
  end
  NArrayMiss.to_nam_no_dup(*arg)
end

.to_nam_no_dup(*arg) ⇒ Object 



150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
 
# File 'lib/narray_miss.rb', line 150

def self.to_nam_no_dup(*arg)
  if arg.length > 2 || arg.length==0 then
    raise("NArrayMiss.to_nar( array [,mask]] )")
  end

  array = arg[0]
  if Numeric===array then array = NArray[array] end
  if Array===array then array = NArray.to_na(array) end
  if !array.is_a?(NArray) then
    raise("argument must be Numeric, NArray or Array")
  end

  if arg.length==2 then
    mask = arg[1]
    if Numeric===mask then mask = array.ne(mask) end
    if Array===mask then
	mask = NArray.to_na(mask).ne(0)
    end
    if mask.class == FalseClass then
	mask = NArray.byte(*array.shape)
    end
    if mask.class == TrueClass then
	mask = NArray.byte(*array.shape).fill(1)
    end
    if !(NArray===mask && mask.typecode==BYTE) then
 raise("mask must be Numeric, Array, true, false or NArray(byte)")
    end
    if mask.length!=array.length
      raise "mask.length must be same as array.length"
    end
  else
    mask = NArray.byte(*array.shape).fill(1)
  end
  __new__(array,mask)
end

Instance Method Details

#-@ ⇒ Object

Arithmetic operator

— NArrayMiss#-@ — NArrayMiss#+(other) — NArrayMiss#-(other) — NArrayMiss#*(other) — NArrayMiss#/(other) — NArrayMiss#%(other) — NArrayMiss#**(other) — NArrayMiss#abs — NArrayMiss#add!(other) — NArrayMiss#sbt!(other) — NArrayMiss#mul!(other) — NArrayMiss#div!(other) — NArrayMiss#mod!(other) — NArrayMiss#mul_add(other, dim, …)



384
385
386
387
388
 
# File 'lib/narray_miss.rb', line 384

def -@
  array = @array.dup
  array[@mask] = -@array[@mask]
  NArrayMiss.to_nam_no_dup(array, @mask.dup)
end

#[](*arg) ⇒ Object

Slicing Array

— NArrayMiss#[](index)

return the value at [((|index|))].
((|index|)) must be (({Integer, Range, Array, true})).
Index order is FORTRAN type.

— NArrayMiss#slice(index)

same as (({NArrayMiss#[]})) but keeps the rank of original array by not elimiting dimensions whose length became equal to 1 (which (({NArrayMiss#[]})) dose).
This is not the case with the 1-dimensional indexing and masking.

— NArrayMiss#set_without_validation(index,value)

replace elements at ((|index|)) by ((|value|)).

— NArrayMiss#[]=(index, value)

replace elements at ((|index|)) by ((|value|)) and
make replaced elements valid.



262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
 
# File 'lib/narray_miss.rb', line 262

def [](*arg)
  if arg[0].class == NArrayMiss && arg[0].typecode == BYTE
    obj = @array[arg[0].to_na(0)]
    if Numeric===obj
      return obj
    else
      return NArrayMiss.to_nam_no_dup(obj)
    end
  else
    obj = @array[*arg]
    if Numeric===obj
      return obj
    else
      return NArrayMiss.to_nam_no_dup(obj,@mask[*arg])
    end
  end
end

#[]=(*arg) ⇒ Object 



306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
 
# File 'lib/narray_miss.rb', line 306

def []=(*arg)
  if arg.length == 2 && arg[0].class == NArrayMiss && arg[0].typecode == BYTE
    idx = arg[0].to_na(0)
    self.set_without_validation(idx,arg[-1])
    if arg[-1].class != NArrayMiss && arg[-1] then
      @mask[idx] = 1
    end
  else
    self.set_without_validation(*arg)
    if arg[-1].class != NArrayMiss && arg[-1] then
      if arg.length==1 then
        @mask=1
      else
        @mask[*arg[0..-2]] = 1
      end
    end
  end
  return self
end

#__clone__ ⇒ Object 



1181
 
# File 'lib/narray_miss.rb', line 1181

alias __clone__ clone

#_dump(limit) ⇒ Object 



1264
1265
1266
 
# File 'lib/narray_miss.rb', line 1264

def _dump(limit)
  Marshal::dump([@array._dump(nil),@mask._dump(nil)])
end

#abs ⇒ Object 



404
405
406
407
408
 
# File 'lib/narray_miss.rb', line 404

def abs
  array = @array.dup
  array[@mask] = @array[@mask].abs
  NArrayMiss.to_nam_no_dup(array, @mask.dup)
end

#accum(*arg) ⇒ Object

Statistics

— NArrayMiss#sum(dim, … [“min_count”=>i])

return summation of elements in specified dimensions.
Elements at which the number of elements for summation is less than ((|i|)) is invalid.

— NArrayMiss#accum(dim, …)

same as (({NArrayMiss#sum})) but not elimiting dimensions whose length became equal to 1.

— NArrayMiss#min(dim, …)

return minimum in specified dimensions.
Elements at which the number of valid elements in the dimension is less than ((|i|)) is invalid.

— NArrayMiss#max(dim, …)

return maximum in specified dimensions.
Elements at which the number of valid elements in the dimension is less than ((|i|)) is invalid.

— NArrayMiss#median(dim, …)

return median in specified dimensions.
Elements at which the number of valid elements in the dimension is less than ((|i|)) is invalid.

— NArrayMiss#mean(dim, …)

return mean of elements in specified dimensions.
Elements at which the number of elements for then mean is less than ((|i|)) is invalid.

— NArrayMiss#stddev(dim, …)

return standard deviation of elements in specified dimensions.
Elements at which the number of elements for calculation the standard deviation is less than ((|i|)) is invalid.



536
537
538
539
540
541
542
543
544
545
 
# File 'lib/narray_miss.rb', line 536

def accum(*arg)
  if @mask.count_true == 0 then
    return nil
  else
    array = @array.dup
    array[@mask.not] = 0
    return NArrayMiss.to_nam_no_dup(array.accum(*arg),
	       @mask.to_type(NArray::INT).accum(*arg).ne(0))
  end
end

#all? ⇒ Boolean

Returns:

  • (Boolean) 


929
930
931
 
# File 'lib/narray_miss.rb', line 929

def all?
  @array[@mask].all?
end

#all_invalid ⇒ Object 



1068
1069
1070
1071
 
# File 'lib/narray_miss.rb', line 1068

def all_invalid
  @mask[true]=0
  self
end

#all_valid ⇒ Object 



1064
1065
1066
1067
 
# File 'lib/narray_miss.rb', line 1064

def all_valid
  @mask[true]=1
  self
end

#all_valid? ⇒ Boolean

Returns:

  • (Boolean) 


1128
1129
1130
 
# File 'lib/narray_miss.rb', line 1128

def all_valid?
  @mask.all?
end

#angle ⇒ Object 



968
969
970
 
# File 'lib/narray_miss.rb', line 968

def angle
  NArrayMiss.to_nam_no_dup(@array.angle,@mask)
end

#any? ⇒ Boolean

Returns:

  • (Boolean) 


932
933
934
 
# File 'lib/narray_miss.rb', line 932

def any?
  @array[@mask].any?
end

#any_valid? ⇒ Boolean

Returns:

  • (Boolean) 


1135
1136
1137
 
# File 'lib/narray_miss.rb', line 1135

def any_valid?
  @mask.any?
end

#clone ⇒ Object 



1182
1183
1184
1185
1186
1187
 
# File 'lib/narray_miss.rb', line 1182

def clone
  obj = __clone__
  obj.set_array(@array.clone)
  obj.set_mask(@mask.clone)
  return obj
end

#coerce(x) ⇒ Object 



1189
1190
1191
1192
1193
1194
1195
1196
1197
 
# File 'lib/narray_miss.rb', line 1189

def coerce(x)
  if Numeric===x then
    return [NArrayMiss.new(NArray[x].typecode,*self.shape).fill(x),self]
  elsif x.class==Array || x.class==NArray then
    return [NArrayMiss.to_nam(x), self]
  else
    raise("donnot know how to cange #{x.class} to NArrayMiss")
  end
end

#collect(&blk) ⇒ Object 



881
882
883
 
# File 'lib/narray_miss.rb', line 881

def collect(&blk)
  self.dup.collect!(&blk)
end

#collect! ⇒ Object 



875
876
877
878
879
880
 
# File 'lib/narray_miss.rb', line 875

def collect!
  for i in 0..self.total-1
    self[i] = yield(self[i])
  end
  self
end

#complex? ⇒ Boolean

Returns:

  • (Boolean) 


1172
1173
1174
 
# File 'lib/narray_miss.rb', line 1172

def complex?
  @array.complex?
end

#conj ⇒ Object 



965
966
967
 
# File 'lib/narray_miss.rb', line 965

def conj
  NArrayMiss.to_nam_no_dup(@array.conj,@mask)
end

#count_false ⇒ Object

— NArrayMiss#count_false

return the number of elements whose value==0 and valid.

— NArrayMiss#count_true

return the number of elements whose value!=0 and valid.

— NArrayMiss#mask(mask)

return (({NArrayMiss#get_mask!&((|mask|))})).

— NArrayMiss#all?

return true if all the valid elements are not 0, else false.

— NArrayMiss#any?

return true if any valid element is not 0, else false.

— NArrayMiss#none?

return true if none of the valid elements is not 0, else false.

— NArrayMiss#where

return (({NArray})) of indices where valid elements are not 0.

— NArrayMiss#where2

return (({Array})) including two (({NArray}))s of indices,
where valid elements are not 0, and 0, respectively.



907
908
909
910
911
912
913
 
# File 'lib/narray_miss.rb', line 907

def count_false
  if @array.typecode==BYTE then
    return @array.count_false-@mask.count_false
  else
    raise("cannot count_true NArrayMiss except BYTE type")
  end
end

#count_invalid(*arg) ⇒ Object 



1146
1147
1148
1149
1150
1151
1152
1153
 
# File 'lib/narray_miss.rb', line 1146

def count_invalid(*arg)
  if arg.length==0 then
    return @mask.count_false
  else
    return NArray.int(*@mask.shape).fill(1).sum(*arg)-
    @mask.to_type(NArray::INT).sum(*arg)
  end
end

#count_true ⇒ Object 



914
915
916
917
918
919
920
 
# File 'lib/narray_miss.rb', line 914

def count_true
  if @array.typecode==BYTE then
    return (@array&@mask).count_true
  else
    raise("cannot count_true NArrayMiss except BYTE type")
  end
end

#count_valid(*arg) ⇒ Object 



1139
1140
1141
1142
1143
1144
1145
 
# File 'lib/narray_miss.rb', line 1139

def count_valid(*arg)
  if arg.length==0 then
    return @mask.count_true
  else
    return @mask.to_type(NArray::INT).sum(*arg)
  end    
end

#dim ⇒ Object

NArrayMiss information

— NArrayMiss#dim

return the dimension which is the number of indices.

— NArrayMiss#rank

same as (({NArrayMiss#dim})).

— NArrayMiss#shape

return the (({Array})) of sizes of each index.

— NArrayMiss#size

return the number of total elements.

— NArrayMiss#total

alias to size

— NArrayMiss#length

alias to size

— NArrayMiss#rank_total

return the number of total of the shape.

— NArrayMiss#typecode

return the typecode.



222
223
224
 
# File 'lib/narray_miss.rb', line 222

def dim
  @array.dim
end

#dup ⇒ Object 



1177
1178
1179
 
# File 'lib/narray_miss.rb', line 1177

def dup
  NArrayMiss.to_nam(@array,@mask)
end

#each ⇒ Object

Iteration

— NArrayMiss#each{|x| …} — NArrayMiss#each_valid{|x| …} — NArrayMiss#each_valid_with_index{|x,i| …} — NArrayMiss#collect{|x| …} — NArrayMiss#collect{|x| …}



860
861
862
863
864
 
# File 'lib/narray_miss.rb', line 860

def each
  for i in 0..self.total-1
    yield(@array[i])
  end
end

#each_valid ⇒ Object 



865
866
867
868
869
 
# File 'lib/narray_miss.rb', line 865

def each_valid
  for i in 0..self.total-1
    yield(@array[i]) if @mask[i]
  end
end

#each_valid_with_index ⇒ Object 



870
871
872
873
874
 
# File 'lib/narray_miss.rb', line 870

def each_valid_with_index
  for i in 0..self.total-1
    yield(@array[i],i) if @mask[i]
  end
end

#get_array ⇒ Object 



1113
1114
1115
 
# File 'lib/narray_miss.rb', line 1113

def get_array
  @array.dup
end

#get_array! ⇒ Object 



1110
1111
1112
 
# File 'lib/narray_miss.rb', line 1110

def get_array!
  @array
end

#get_mask ⇒ Object 



1107
1108
1109
 
# File 'lib/narray_miss.rb', line 1107

def get_mask
  @mask.dup
end

#get_mask! ⇒ Object 



1104
1105
1106
 
# File 'lib/narray_miss.rb', line 1104

def get_mask!
  @mask
end

#hton ⇒ Object Also known as: ntoh 



999
1000
1001
 
# File 'lib/narray_miss.rb', line 999

def hton
  NArrayMiss.to_nam(@array.hton,@mask.hton)
end

#htov ⇒ Object Also known as: vtoh 



1003
1004
1005
 
# File 'lib/narray_miss.rb', line 1003

def htov
  NArrayMiss.to_nam(@array.htov,@mask.htov)
end

#im ⇒ Object 



975
976
977
 
# File 'lib/narray_miss.rb', line 975

def im
  NArrayMiss.to_nam_no_dup(@array.im,@mask)
end

#imag ⇒ Object 



962
963
964
 
# File 'lib/narray_miss.rb', line 962

def imag
  NArrayMiss.to_nam_no_dup(@array.imag,@mask)
end

#imag=(arg) ⇒ Object 



971
972
973
974
 
# File 'lib/narray_miss.rb', line 971

def imag=(arg)
  @array.image=(arg)
  self
end

#inspect ⇒ Object 



1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
 
# File 'lib/narray_miss.rb', line 1200

def inspect
#    "array -> " + @array.inspect + "\nmask  -> " + @mask.inspect
  count_line = 0
  max_line = 10
  max_col = 80
  sep = ", "
  const = Hash.new
  NArray.constants.each{|c| const[NArray.const_get(c)] = c}
  str_ret = "NArrayMiss."+const[typecode].to_s.downcase+"("+shape.join(",")+"):"
  if rank == 0 then
    str_ret += " []"
    return str_ret
  else
     str_ret += "\n"
  end
  str = ""
  index = Array.new(rank,0)
  index[0] = true
  i = 1
  (rank-1).times{ str_ret += "[ " }
  while(true)
    i.times{ str_ret += "[ " }

    str = @array[*index].inspect
    ary = str[str.index("[")+1..str.index("]")-1].strip.split(/\s*,\s*/)
    miss = @mask[*index].where2[1]
    miss = miss[miss<ary.length].to_a
    if ary[-1]=="..." && miss[-1]==ary.length-1 then miss.pop end
    for j in miss
	ary[j] = "-"
    end
    while ( rank*4+ary.join(", ").length > max_col )
      ary.pop
      ary[-1] = "..."
    end
    str_ret += ary.join(", ")
    i = 1
    while (i<rank)
	if index[i]<shape[i]-1 then
 str_ret += " ]"+sep+"\n"
 count_line += 1
 index[i] += 1
 break
	else
 str_ret += " ]"
 index[i] = 0
 i += 1
	end
    end

    if i>=rank then
	str_ret += " ]"
	return str_ret
    elsif count_line>=max_line then
	str_ret += " ..."
	return str_ret
    end

    (rank-i).times{ print("  ") }
  end
  return str_ret
end

#integer? ⇒ Boolean

Others

— NArrayMiss#integer?

return true if (({NArrayMiss})) is byte, sint or int, else false.

— NArrayMiss#complex?

return true if (({NArrayMiss})) is scomplex or complex, else false.

— NArrayMiss#dup — NArrayMiss#coerce(object) — NArrayMiss#inspect

go back to ((<Index>))

Returns:

  • (Boolean) 


1169
1170
1171
 
# File 'lib/narray_miss.rb', line 1169

def integer?
  @array.integer?
end

#mask(arg) ⇒ Object 



921
922
923
924
925
926
927
 
# File 'lib/narray_miss.rb', line 921

def mask(arg)
  obj = self.dup
  if arg.class==NArrayMiss then
    arg = arg.get_array!&arg.get_mask!
  end
  obj.set_mask(@mask&arg)
end

#mean(*arg) ⇒ Object 



587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
 
# File 'lib/narray_miss.rb', line 587

def mean(*arg)
  if @mask.count_true == 0 then
    return nil
  end
  min_count = 1
  options=["min_count"]
  if arg.length!=0 && arg[-1].class==Hash then
    option = arg[-1]
    arg2=arg[0..-2]
    option.each_key{|key|
	if !options.index(key) then
 raise(ArgumentError,key+" option is not exist")
	end
    }
    min_count = option["min_count"]
  else
    arg2=arg
  end
  count = @mask.to_type(NArray::INT).sum(*arg2)
  if count.class == NArray then
    count = NArrayMiss.to_nam(count,count.ge(min_count))
    if count.ge(min_count).count_true == 0
      return nil
    else
      return self.sum(*arg2)/count
    end
  else
    if count < min_count then
      return nil
    else
	return self.sum(*arg2)/count
    end
  end
end

#median(*arg) ⇒ Object 



662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
 
# File 'lib/narray_miss.rb', line 662

def median(*arg)
  if arg.length==0 then
    return @array[@mask].median
  else
    nshape = NArray.to_na(@array.shape)
    nshape[arg]=1
    nslice = nshape[nshape.ne(1).where]
    index = NArray.object(@mask.rank)
    index[nshape.eq(1).where] = true
    obj = NArrayMiss.new(@array.typecode,*nslice.to_a)
    total = 1
    nslice.each{|n| total *= n}
    for i in 0...total
      index[nshape.ne(1).where] = pos(i,nslice)
      mask = NArray.byte(*@array.shape).fill(0)
      mask[*index] = 1
      mask = @mask&mask
      if mask.count_true != 0 then
        obj[*pos(i,nslice)] = @array[mask].median
      end
    end
    return obj
  end
end

#mul_add(*arg) ⇒ Object 



420
421
422
423
424
425
426
 
# File 'lib/narray_miss.rb', line 420

def mul_add(*arg)
  if arg.length==1 then
    return (self*arg[0]).sum
  else
    return (self*arg[0]).sum(*arg[1..-1])
  end
end

#newdim(*arg) ⇒ Object Also known as: rewrank 



779
780
781
782
 
# File 'lib/narray_miss.rb', line 779

def newdim(*arg)
  obj = self.dup
  obj.newdim!(*arg)
end

#newdim!(*arg) ⇒ Object Also known as: rewrank!, rewrank= 



772
773
774
775
776
 
# File 'lib/narray_miss.rb', line 772

def newdim!(*arg)
  @array = @array.newdim!(*arg)
  @mask = @mask.newdim!(*arg)
  self
end

#none? ⇒ Boolean

Returns:

  • (Boolean) 


935
936
937
 
# File 'lib/narray_miss.rb', line 935

def none?
  @array[@mask].none?
end

#none_valid? ⇒ Boolean Also known as: all_invalid?

Returns:

  • (Boolean) 


1131
1132
1133
 
# File 'lib/narray_miss.rb', line 1131

def none_valid?
  @mask.none?
end

#not ⇒ Object 



499
500
501
 
# File 'lib/narray_miss.rb', line 499

def not
  NArrayMiss.to_nam_no_dup(@array.not, @mask.dup)
end

#randomn ⇒ Object 



359
360
361
362
363
 
# File 'lib/narray_miss.rb', line 359

def randomn
  obj = self.dup
  obj[@mask] = @array[@mask].randomn
  obj
end

#rank ⇒ Object 



225
226
227
 
# File 'lib/narray_miss.rb', line 225

def rank
  @array.rank
end

#rank_total(*arg) ⇒ Object 



237
238
239
 
# File 'lib/narray_miss.rb', line 237

def rank_total(*arg)
  @array.rank_total(*arg)
end

#real ⇒ Object

Complex compound number (only for scomplex and complex)

— NArrayMiss#real — NArrayMiss#imag — NArrayMiss#conj — NArrayMiss#angle — NArrayMiss#imag=(other) — NArrayMiss#im



959
960
961
 
# File 'lib/narray_miss.rb', line 959

def real
  NArrayMiss.to_nam_no_dup(@array.real,@mask)
end

#reshape(*arg) ⇒ Object 



767
768
769
770
 
# File 'lib/narray_miss.rb', line 767

def reshape(*arg)
  obj = self.dup
  obj.reshape!(*arg)
end

#reshape!(*arg) ⇒ Object Also known as: shape=

Changing Shapes of indices

— NArrayMiss#reshape!(size, …)

change shape of array.

— NArrayMiss#reshape(size, …)

same as (({NArrayMiss#reshape!})) but create new object.

— NArrayMiss#shape=(size, …)

same as (({NArrayMiss#reshape!})).

— NArrayMiss#newdim!(dim)

insert new dimension with size=1

— NArrayMiss#newdim(dim)

same as (({NArrayMiss#newdim!})) but create new object.

— NArrayMiss#rewrank!(dim)

same as (({NArrayMiss#newdim!})).

— NArrayMiss#rewrank=(dim)

same as (({NArrayMiss#newdim!})).



762
763
764
765
766
 
# File 'lib/narray_miss.rb', line 762

def reshape!(*arg)
  @array = @array.reshape!(*arg)
  @mask = @mask.reshape!(*arg)
  self
end

#set_invalid(*pos) ⇒ Object Also known as: invalidation 



1059
1060
1061
1062
 
# File 'lib/narray_miss.rb', line 1059

def set_invalid(*pos)
  @mask[*pos] = 0
  self
end

#set_mask(mask) ⇒ Object 



1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
 
# File 'lib/narray_miss.rb', line 1072

def set_mask(mask)
  if mask.class == Array then
    tmp = NArray.byte(*@mask.shape)
    tmp[true] = mask
    mask = tmp
  end
  if mask.class == NArrayMiss then
    mask = mask.to_na(0)
  end
  if mask.class == NArray then
    if mask.typecode != BYTE then
	raise("mask must be NArrayMiss.byte, NArray.byte or Array")
    end
    if @array.shape != mask.shape then
	raise("mask.shape must be same as array")
    end
    @mask = mask.dup
    return self
  else
    raise("mask must be NArray.byte or Array")
  end
end

#set_missing_value(val) ⇒ Object 



1099
1100
1101
1102
 
# File 'lib/narray_miss.rb', line 1099

def set_missing_value(val)
  obj = self.dup
  obj.set_missing_value!(val)
end

#set_missing_value!(val) ⇒ Object 



1095
1096
1097
1098
 
# File 'lib/narray_miss.rb', line 1095

def set_missing_value!(val)
  @array[@mask.not] = val
  self
end

#set_valid(*pos) ⇒ Object Also known as: validation

Mask and missing value

— NArrayMiss#set_valid(index)

validate element at ((|index|)).
((|index|)) must be (({Integer, Range, Array, or ture})).

— NArrayMiss#validation(index)

alias to set_valid

— NArrayMiss#set_invalid(index)

invaliadate element at ((|index|)).
((|index|)) must be (({Integer, Range, Array, or ture})).

— NArrayMiss#invalidation(index)

alias to set_invalid

— NArrayMiss#all_valid

set all elements valid

— NArrayMiss#all_invalid

set all elements invalid

— NArrayMiss#set_mask(mask)

masking by ((|mask|))

— NArrayMiss#set_missing_value(value)

replace invalid elements by ((|value|)).

— NArrayMiss#get_mask!

return (({NArray})) of byte as mask.

— NArrayMiss#get_mask

return (({NArray})) of byte as mask.

— NArrayMiss#get_array!

return (({NArray})) as data.

— NArrayMiss#get_array

return (({NArray})) as data.

— NArrayMiss#valid?

return (({Array})) whose elements are true or false corresponding to valid or invalid of elements, respectively.

— NArrayMiss#all_valid?

return true if all elements are valid, else false.

— NArrayMiss#none_valid?

return true if all elements are invalid, else false.

— NArrayMiss#all_invalid?

alias to none_valid?

— NArrayMiss#any_valid?

return true if any elements are valid, else false.

— NArrayMiss#count_valid

return the number of valid elements.

— NArrayMiss#count_invalid

return the number of invalid elements.



1054
1055
1056
1057
 
# File 'lib/narray_miss.rb', line 1054

def set_valid(*pos)
  @mask[*pos] = 1
  self
end

#set_without_validation(*arg) ⇒ Object 



283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
 
# File 'lib/narray_miss.rb', line 283

def set_without_validation(*arg)
  if arg.length==1 then
    if !arg[0] then
	@mask[] = 0
    elsif arg[0].class == NArrayMiss then
	@array[] = arg[0].get_array!
	@mask[] = arg[0].get_mask!
    else
	@array[] = arg[0]
    end
  else
    if !arg[-1] then
	@mask[*arg[0..-2]] = 0
    elsif arg[-1].class == NArrayMiss then
	@array[*arg[0..-2]] = arg[-1].get_array!
	@mask[*arg[0..-2]] = arg[-1].get_mask!
    else
	@array[*arg[0..-2]] = arg[-1]
    end
  end
  return self
end

#shape ⇒ Object 



228
229
230
 
# File 'lib/narray_miss.rb', line 228

def shape
  @array.shape
end

#size ⇒ Object Also known as: total, length 



231
232
233
 
# File 'lib/narray_miss.rb', line 231

def size
  @array.size
end

#slice(*arg) ⇒ Object 



279
280
281
 
# File 'lib/narray_miss.rb', line 279

def slice(*arg)
  NArrayMiss.to_nam_no_dup(@array.slice(*arg),@mask.slice(*arg))
end

#stddev(*arg) ⇒ Object 



621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
 
# File 'lib/narray_miss.rb', line 621

def stddev(*arg)
  if @mask.count_true == 0 then
    return nil
  end
  min_count=2
  options=["min_count"]
  if arg.length!=0 && arg[-1].class==Hash then
    option = arg[-1]
    arg = arg[0..-2]
    option.each_key{|key|
	if !options.index(key) then
 raise(ArgumentError,key+" option is not exist")
	end
    }
    min_count = option["min_count"].to_i
    if min_count<2
      raise(ArgumentError, "min_count must be >= 2")
    end
  end
  count = @mask.to_type(NArray::INT).sum(*arg)
  count2 = @mask.to_type(NArray::INT).accum(*arg)
  if count.class==NArray then
    count = NArrayMiss.to_nam_no_dup(count,count.ge(min_count))
    if count.get_mask!.count_true == 0
      return nil
    end
  else
    if count < min_count then
	return nil
    end
  end
  if self.integer? then
    a = self.to_f
  else
    a = self
  end
  var = ( (a-a.accum(*arg)/count2)**2 ).sum(*arg)/(count-1)
  obj = NMMath::sqrt(var)
  return obj
end

#swap_byte ⇒ Object

Byte swap

— NArrayMiss#swap_byte

swap byte order.

— NArrayMiss#hton

convert to network byte order.

— NArrayMiss#ntoh

convert from network byte order.

— NArrayMiss#htov

convert to VAX byte order.

— NArrayMiss#vtoh

convert from VAX byte order.



994
995
996
997
998
 
# File 'lib/narray_miss.rb', line 994

def swap_byte
  obj = self.dup
  obj.set_without_validation(@array.swap_byte)
  obj
end

#to_a ⇒ Object 



824
825
826
 
# File 'lib/narray_miss.rb', line 824

def to_a
  @array.to_a
end

#to_na(*arg) ⇒ Object 



837
838
839
 
# File 'lib/narray_miss.rb', line 837

def to_na(*arg)
  return self.dup.to_na!(*arg)
end

#to_na!(*arg) ⇒ Object 



827
828
829
830
831
832
833
834
835
836
 
# File 'lib/narray_miss.rb', line 827

def to_na!(*arg)
  if arg.length==0
    return @array
  elsif arg.length==1 then
    self.set_missing_value!(arg[0])
    return @array
  else
    raise(ArgumentError, "Usage: NArray#to_na([missing_value])")
  end
end

#to_s ⇒ Object 



840
841
842
 
# File 'lib/narray_miss.rb', line 840

def to_s
  @array.to_s
end

#to_string ⇒ Object 



843
844
845
846
847
848
 
# File 'lib/narray_miss.rb', line 843

def to_string
  obj = NArrayMiss.obj(*@array.shape)
  obj.set_without_validation( @array.to_string )
  obh.set_mask(@mask)
  obj
end

#to_type(typecode) ⇒ Object 



821
822
823
 
# File 'lib/narray_miss.rb', line 821

def to_type(typecode)
  NArrayMiss.to_nam_no_dup(@array.to_type(typecode), @mask.dup)
end

#transpose(*arg) ⇒ Object

Transpose

— NArrayMiss#transpose(dim0, dim1, …)

transpose array. The 0-th dimension goes to the ((|dim0|))-th dimension of new array.



734
735
736
737
738
739
740
741
 
# File 'lib/narray_miss.rb', line 734

def transpose(*arg)
  obj = self.dup
  shape = arg.collect{|i| obj.shape[i]}
  obj.reshape!(*shape)
  obj.set_without_validation( @array.transpose(*arg) )
  obj.set_mask(@mask.transpose(*arg))
  obj
end

#typecode ⇒ Object 



241
242
243
 
# File 'lib/narray_miss.rb', line 241

def typecode
  @array.typecode
end

#valid? ⇒ Boolean

Returns:

  • (Boolean) 


1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
 
# File 'lib/narray_miss.rb', line 1117

def valid?
  where = self.get_mask!.where2
  tf = Array.new(self.total)
  for i in where[0]
    tf[i] = true
  end
  for i in where[1]
    tf[i] = false
  end
  tf
end

#where ⇒ Object 



939
940
941
 
# File 'lib/narray_miss.rb', line 939

def where
  (@array&@mask).where
end

#where2 ⇒ Object 



942
943
944
 
# File 'lib/narray_miss.rb', line 942

def where2
  self.where-@mask.where
end

#~@ ⇒ Object

Bitwise operator (only for byte, sint and int)

— NArrayMiss#~@ — NArrayMiss#&(other) — NArrayMiss#|(other) — NArrayMiss#^(other)



437
438
439
 
# File 'lib/narray_miss.rb', line 437

def ~@
  NArrayMiss.to_nam_to_dup(~@array, @mask.dup)
end