动态数组继续赋值问题

快乐枫

一个动态数组 a ,开始分配内存 allocate(a(128))
之后为a赋值，如果超过了128个，继续给a 分配128，也就是a变成 a（256）,之前的赋值也保存

这样可以实现吗？

Mid_Farmer

按楼主的说法来看 楼主已经先确定了算法 算法肯定能实现 只是看怎么编程和程序的效率高不高而已 如果你把数据形式拿来看看 说不定有更多的算法告诉你

快乐枫

Mid_Farmer 发表于 2014-8-10 10:39
按楼主的说法来看 楼主已经先确定了算法 算法肯定能实现 只是看怎么编程和程序的效率高不高而已 如果你把数 ...

可以实现是吗？你指的数据形式是？

斥鷃

声明的时候就是a(:)这种的，没有128，你都要开始赋值了你还不知道数组维数？应该可以用某种方式确定要赋值的数是否是超过128的，要是超过就直接allocate(a(256))就好了。你论述的那种方式fortran实现感觉是不是要转存一个变量~因为上面你已经用allocate声明过了，重新改维数不知道行不行。

快乐枫

斥鷃 发表于 2014-8-10 11:22
声明的时候就是a(:)这种的，没有128，你都要开始赋值了你还不知道数组维数？应该可以用某种方式确定要赋值 ...

提前是多少不知道，是一边计算一边赋值，超了128再扩展数组的范围，维数知道，是否超过256也不确定，超过了继续扩，所以直接allocate(a(256))不能解决问题

斥鷃

你计算的东西贴上来，才好想法子确定维数，要是计算不复杂的话让它算一遍，记一下数，然后再重定义变量，再算一遍赋值，你的计算是写在while循环里面？

Mid_Farmer

快乐枫 发表于 2014-8-10 10:50
可以实现是吗？你指的数据形式是？

你还是贴数据吧 楼上的说得很清楚了

nuist2007

学以致用

lysx

快乐枫 发表于 2014-8-10 14:25
提前是多少不知道，是一边计算一边赋值，超了128再扩展数组的范围，维数知道，是否超过256也不确定，超过 ...

用指针和链表啊

快乐枫

斥鷃 发表于 2014-8-10 15:47
你计算的东西贴上来，才好想法子确定维数，要是计算不复杂的话让它算一遍，记一下数，然后再重定义变量，再 ...

是用的while循环,如下，在update the solution部分，是每循环一步，更新一下Z 和Y 的数据，提前给Z和Y一个范围，用if判断n 大于它们的维数了，如果超了，再往后扩展，不知如何实现。

因为不止这一个地方用，先记一下数再重新算一遍的方法不太可取，不知道 有没有直接解决它的方法呢

do while((t<tfinal).and.(t+h>t))
   
    if(t+h>tfinal)then

      h=tfinal-t

    end if

!----------------------------------------------------------------------!
! Compute the slopes

    call iceocean(t,y1,temp)
    f(:,1)=temp

    do j=1,5
      
      call iceocean(t+alpha(j)*h,y1+matmul(h*f,beta(:,j)),temp)
      f(:,j+1)=temp
   
    end do

!----------------------------------------------------------------------!
! Estimate the error and the acceptable error

    delta=sqrt(maxval(matmul(h*f,gama(:,2))*matmul(h*f,gama(:,2))))
    tau=tol*max(sqrt(maxval(y1*y1)),1.0)

!----------------------------------------------------------------------!
! Update the solution only if the error is acceptable

    if(delta<=tau)then

      t=t+h
      y1=y1+h*matmul(f,gama(:,1))
      call iceocean(t,y1,temp)
      n=n+1

      if(n>size(Z))
          ?
      end if

      Z(n)=t
      Y(n,:)=y1

    end if

    if(delta/=0.0)then

      h=min(hmax, 0.8*h*(tau/delta)**pow)

    end if

!----------------------------------------------------------------------!
! end the main loop

  end do