Changeset 174f4b6 in tol


Ignore:
Timestamp:
Dec 11, 2009, 9:44:17 AM (17 years ago)
Author:
Víctor de Buen <vdebuen@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
db53538
Parents:
a1affcf
Message:

Counting number of nodal, stress and border points for each reference point

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@1792 be4846ae-f051-483a-831f-b64e7d08527f

File:
1 edited

Legend:

Unmodified
Added
Removed

  • tol_pkg/MulVarFunAprox/DesignC2.tol

    ra1affcf r174f4b6  
    3030  Real _.Z.maxNeighbour = ?;
    3131  VMatrix _.Y = Rand(0,2,0,0);
     32  VMatrix _.Yj = Rand(0,1,0,0);
     33  VMatrix _.Yt = Rand(0,1,0,0);
     34  VMatrix _.Yb = Rand(0,1,0,0);
     35  VMatrix _.Ym = Rand(0,1,0,0);
     36  VMatrix _.Yn = Rand(0,1,0,0);
    3237  VMatrix _.Z = Rand(0,2,0,0);
    3338  VMatrix _.Z.h = Rand(0,1,0,0);
     39  VMatrix _.Z.Y.neighbour= Rand(0,1,0,0);
    3440  VMatrix _.B = Rand(0,2,0,0);
    3541  VMatrix _.B.h = Rand(0,0,0,0);
    … …  
    5056  VMatrix _.LQ = Rand(0,0,0,0);
    5157  VMatrix _.Q.nonNullCellsByColumns = Rand(0,0,0,0);
     58
    5259  //////////////////////////////////////////////////////////////////////////////
    5360  NameBlock get.base(Real sel)
    … …  
    7885    Real build.stress.points(0);
    7986    Real build.border.points(0);
     87    VMatrix _.Ym := _.Yj+((_.Yt+_.Yb)*Case(
     88      prior.Z.order< 0, 0,
     89      prior.Z.order==0, 1,
     90      prior.Z.order==1, 1*(1+_.n),
     91      prior.Z.order==2, 1*(1+_.n*3/2+(1/2)*_.n^2)));
    8092    Real _.M := _.J+_.Mb+_.Mz;
    8193    Real build.base(0);
    … …  
    100112    If(!m, Rand(0,_.N,0,0), {
    101113    WriteLn("TRACE "+_MID+"2");
    102     WriteLn("TRACE "+_MID+"2");
     114    WriteLn("TRACE "+_MID+"3");
    103115    VMatrix X.j=DifEq(1/(1-B),Rand(m,1,1,1));
    104     WriteLn("TRACE "+_MID+"2");
     116    WriteLn("TRACE "+_MID+"4");
    105117    Set X.dif = For(1,_.X.maxNeighbour,VMatrix(Real i)
    106118    {
    … …  
    108120      X - SubRow(_.Y,X.h.i)
    109121    });
    110     WriteLn("TRACE "+_MID+"2");
     122    WriteLn("TRACE "+_MID+"5");
    111123    Set X.weight = get.X.weight(X.dif);
    112     WriteLn("TRACE "+_MID+"2");
     124    WriteLn("TRACE "+_MID+"6");
    113125    Matrix QX.triplet = Group("ConcatRows", For(1,_.X.maxNeighbour,Matrix(Real i)
    114126    {
    … …  
    127139      }))
    128140    }));
     141    WriteLn("TRACE "+_MID+"7");
    129142    Pack(Convert(Triplet(QX.triplet, m, _.N),"Cholmod.R.Sparse"))
    130143  })};
    … …  
    274287    WriteLn("TRACE "+_MID+"4");
    275288    VMatrix _.Q := _.QX << _.QB << _.PZ;
    276     WriteLn("TRACE "+_MID+"7");
     289    WriteLn("TRACE "+_MID+"5");
    277290    VMatrix _.Q.nonNullCellsByColumns :=
    278291    {
    … …  
    283296      }))
    284297    };
    285     WriteLn("TRACE "+_MID+"8");
     298    WriteLn("TRACE "+_MID+"6");
    286299    VMatrix _.LQ := CholeskiFactor(_.Q,"XtX");
    287     WriteLn("TRACE "+_MID+"9");
     300    WriteLn("TRACE "+_MID+"7");
    288301    True
    289302  }
Note: See TracChangeset for help on using the changeset viewer.