Class: NumRu::VArrayComposite

Inherits:
VArray
  • Object
show all
Defined in:
lib/numru/gphys/varraycomposite.rb

Constant Summary

Constants inherited from VArray

NumRu::VArray::Binary_operators, NumRu::VArray::Binary_operatorsL, NumRu::VArray::Binary_operatorsL_comp, NumRu::VArray::Binary_operatorsL_other, NumRu::VArray::Binary_operators_Uconv, NumRu::VArray::Binary_operators_Unone, NumRu::VArray::Binary_operators_Uop, NumRu::VArray::Math_funcs, NumRu::VArray::Math_funcs_nondim, NumRu::VArray::Math_funcs_radian, NumRu::VArray::NArray_type1_methods, NumRu::VArray::NArray_type2_methods, NumRu::VArray::NArray_type2_methods1, NumRu::VArray::NArray_type2_methods2, NumRu::VArray::NArray_type2_methods3, NumRu::VArray::NArray_type3_methods, NumRu::VArray::Unary_operators

Instance Method Summary collapse

Methods inherited from VArray

#+@, #att_names, #attr_copy, #axis_cyclic?, #axis_cyclic_extendible?, #axis_draw_positive, #axis_modulo, #bin_mean, #bin_sum, #coerce, #convert_units, #copy, #del_att, #get_att, #long_name, #long_name=, #mapped?, #marshal_dump, #marshal_load, #name, new2, #rename!, #replace_val, #reshape, #reshape!, #sqrt, #to_g1D, #transpose, #units, #units=

Constructor Details

#initialize(varrays) ⇒ VArrayComposite

Returns a new instance of VArrayComposite.



72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
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
145
146
147
148
149
150
 
# File 'lib/numru/gphys/varraycomposite.rb', line 72

def initialize( varrays )

  if !varrays.is_a?(NArray) || varrays.typecode != NArray::OBJECT
  raise ArgumentError, "argument must be a NArray of VArray (or nil)"
  end

  nvas = varrays.shape
  vrank = nvas.length
  varrays.each{|va|
  if !va.nil?
    @first_vary = va
    @attr = @first_vary.attr_copy
    @name = @first_vary.name
    @crank = va.rank   # rank of the component VArrays
    break
  end
  }
  if vrank > @crank
    @rank = vrank       #=>  @rank > @crank
  else
  @rank = @crank
    if vrank < @crank
      (@crank - vrank).times{
        varrays = varrays.newdim(vrank)
        nvas.push(1)
      }
      vrank = @crank   # rank of varrays is increased to @crank (==@rank)
    end
  end
  un0 = varrays[0].units
  for i in 1...varrays.length
    va = varrays[i]
    varrays[i] = va.convert_units(un0) if va && va.units != un0
  end

  @bound_idx = Array.new     # will be Array of Array
  for dim in 0...vrank
  @bound_idx[dim] = [0]    # the frst element is always 0
  for i in 0...nvas[dim]
    idx = [true]*dim + [i..i,false]
    set=false
    len=0
    varrays[*idx].each{|va|
 if !set && va
          ### commented out for the duck typing:
   #if !va.is_a?(VArray)
          #  raise ArgumentError,"Not a VArray: #{va.inspect}"
   #end
          if dim < @crank
            len = va.shape_current[dim]
          else
            len = 1
          end
   @bound_idx[dim][i+1] = @bound_idx[dim][i] + len
   set=true
 elsif va
          if dim < @crank
            lc = va.shape_current[dim]
          else
            lc = 1
          end
   if lc != len
    raise ArgumentError,"Non-uniformity in the #{i}th element"+
               " of the #{dim}th dimension (#{va.shape_current[dim]} for #{len})"
   end
 end
    }
    if !set
 raise "No VArray is found in the #{i}th elem of the #{dim}th dim"
    end
  end
  end
  #< finish >
  @varrays = varrays.dup
  @shape = @bound_idx.collect{|a| a[-1]}
  @length = 1
  @shape.each{|l| @length *= l}
  # p @bound_idx, @varrays
end

Instance Method Details

#[](*slicer) ⇒ Object 



232
233
234
235
236
237
238
239
240
 
# File 'lib/numru/gphys/varraycomposite.rb', line 232

def [](*slicer)
  slicer = __rubber_expansion(slicer)
  varrays = _div_idx(*slicer)
  if varrays.length == 1
  varrays[0]
  else
  VArrayComposite.new( varrays )
  end
end

#[]=(*args) ⇒ Object 



242
243
244
245
246
247
248
249
250
251
 
# File 'lib/numru/gphys/varraycomposite.rb', line 242

def []=(*args)
  args = __rubber_expansion(args)
  val = args.pop
  slicer = args
  slicer.collect!{|i| (i.is_a?(Numeric)) ? i..i : i }
  subva = self[*slicer]
  val = val.val if val.is_a?(VArray)
  subva.val= val
  val
end

#file ⇒ Object 



264
265
266
 
# File 'lib/numru/gphys/varraycomposite.rb', line 264

def file
  @varrays.collect{|va| va.file}
end

#inspect ⇒ Object 



163
164
165
166
 
# File 'lib/numru/gphys/varraycomposite.rb', line 163

def inspect
  "<#{self.class.to_s} shape=#{shape.inspect} #_of_tiles=#{@bound_idx.collect{|b| b.length-1}.inspect}>"
  # "bounds=#{@bound_idx.collect{|b| b[1..-2]}.inspect}"
end

#length ⇒ Object Also known as: total 



175
176
177
 
# File 'lib/numru/gphys/varraycomposite.rb', line 175

def length
  @length
end

#name=(nm) ⇒ Object

Raises:

  • (ArgumentError) 


268
269
270
271
272
273
 
# File 'lib/numru/gphys/varraycomposite.rb', line 268

def name=(nm)
  raise ArgumentError, "name should be a String" if ! nm.is_a?(String)
  @name = nm
  @varrays.each{|va| va.name=nm}
  nm
end

#ntype ⇒ Object

def ntype

  @first_vary.ntype
end



256
257
258
 
# File 'lib/numru/gphys/varraycomposite.rb', line 256

def ntype
  __ntype(typecode)
end

#rank ⇒ Object 



172
173
174
 
# File 'lib/numru/gphys/varraycomposite.rb', line 172

def rank
  @rank
end

#rename(nm) ⇒ Object 



274
275
276
277
 
# File 'lib/numru/gphys/varraycomposite.rb', line 274

def rename(nm)
  raise "method rename is not vailable. (But rename! and name= are available)"
  ## self.dup.set_name_shallow(nm)
end

#set_att(name, val) ⇒ Object 



154
155
156
157
158
159
 
# File 'lib/numru/gphys/varraycomposite.rb', line 154

def set_att(name,val)
  @varrays.each{|va|
  va.set_att(name,val)
  }
  self
end

#shape ⇒ Object Also known as: shape_current 



168
169
170
 
# File 'lib/numru/gphys/varraycomposite.rb', line 168

def shape
  @shape
end

#typecode ⇒ Object 



260
261
262
 
# File 'lib/numru/gphys/varraycomposite.rb', line 260

def typecode
  @first_vary.typecode
end

#val ⇒ Object 



180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
 
# File 'lib/numru/gphys/varraycomposite.rb', line 180

def val
  val = nil
  loop_multi_dim_index( @varrays.shape ){|index|
  vidx = Array.new
  for d in 0...@rank
      if d < @crank
        vidx[d] = (@bound_idx[d][index[d]])..(@bound_idx[d][index[d]+1]-1)
      else
        vidx[d] = @bound_idx[d][index[d]]
      end
  end
    units = @attr["units"]
    if units
      v = @varrays[*index].convert_units(units).val
    else
      v = @varrays[*index].val
    end
  if !val
    if NArray===v
 val = NArray.new(v.typecode,*@shape)
    elsif NArrayMiss===v
 val = NArrayMiss.new(v.typecode,*@shape)
    else
 raise TypeError, "Unexpected type #{v.class}"
    end
  end
  val[*vidx] = v
  }
  val
end

#val=(narray) ⇒ Object 



211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
 
# File 'lib/numru/gphys/varraycomposite.rb', line 211

def val=(narray)
  loop_multi_dim_index( @varrays.shape ){|index|
  if narray.is_a?(Numeric)
    sub = narray
  else
    __check_ary_class(narray)
    vidx = Array.new
    for d in 0...@rank
        if d < @crank
          vidx[d] = (@bound_idx[d][index[d]])..(@bound_idx[d][index[d]+1]-1)
        else
          vidx[d] = @bound_idx[d][index[d]]
        end
    end
    sub=narray[*vidx]
  end
  @varrays[*index].val = sub
  }
  narray
end