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 



235
236
237
238
239
240
241
242
243
 
# File 'lib/numru/gphys/varraycomposite.rb', line 235

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

#[]=(*args) ⇒ Object 



245
246
247
248
249
250
251
252
253
254
 
# File 'lib/numru/gphys/varraycomposite.rb', line 245

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 



267
268
269
 
# File 'lib/numru/gphys/varraycomposite.rb', line 267

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:



271
272
273
274
275
276
 
# File 'lib/numru/gphys/varraycomposite.rb', line 271

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



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

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 



277
278
279
280
 
# File 'lib/numru/gphys/varraycomposite.rb', line 277

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 



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

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 



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

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
210
211
212
 
# 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
    if val.is_a?(NArray) && v.is_a?(NArrayMiss)
      val = NArrayMiss.to_nam(val)
    end
  val[*vidx] = v
  }
  val
end

#val=(narray) ⇒ Object 



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

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