Results=Point/@Table[Flatten[{Pts[[i]],v[i]*10}],{i,nPts}] //Simplify;(*グラフを適当な縮尺にするためにv[i]に係数をかける*)
Show[Graphics3D[{PointSize[0.05],Results} ],
ViewPoint->{2.980,1.502,0.558},Axes->True,FaceGrids->All]
(*パラメータによる位置リストの作成*)
ParameterPosition=Table[
Flatten[{Apply[Plus,{ Pts[[ ele[[i]] ]] [[1]],
sP*tP*( Pts[[ ele[[i]] ]] [[2]] - Pts[[ ele[[i]] ]] [[1]] ),
sP*(1-tP)*( Pts[[ ele[[i]] ]] [[3]] - Pts[[ ele[[i]] ]] [[1]] )}] }],
{i,1,nElements}];
(*高速化のために重み係数を近似計算*)
kA=Rationalize[ N[k,7],10^-5 ];
ResultV=Table[
Flatten[{
1/10*kA[[i]].Flatten[Append[{1},ParameterPosition[[i]] ]]}],{i,1,nElements}];
ResultVSimp=Rationalize[ Simplify[ N[ResultV,5] ], 10^-5];
ResultVC=Table[Hue[Flatten[ResultVSimp][[i]]/70],{i,1,nElements}]; //Simplify (*グラフ表示のために適当な係数をかける*)
ResultVO=Table[Flatten[{
ParameterPosition[[i]],
ResultVSimp[[i]]*50,
ResultVC[[i]] }],
{i,1,nElements}];(*グラフ表示のために適当な係数をかける*)
ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->5,FaceGrids->All,ViewPoint->{2.446, 1.299, 1.945},
Lighting->False,Shading -> True, PlotRange->All]
Out of memory. Exiting.
ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->10,FaceGrids->All,Shading -> True,
ViewPoint->{0.000,0.000,9.000},AxesLabel->{"X","Y","Z"},Lighting->False,PlotRange->All]
ParametricPlot3D[Evaluate[ResultVO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->10,FaceGrids->All,Lighting->False,
AxesLabel->{"X","Y","Z"},Lighting->False,PlotRange->All,ViewPoint->{-2.392, 2.081, 1.184}]
ResultEC=Table[Hue[N[(Part[ kA[[i]],2]^2+Part[ kA[[i]],3]^2)^(1/2)]/8],{i,1,nElements}]; //Simplify
ResultEO=Table[Flatten[{
ParameterPosition[[i]],
ResultVSimp[[i]]*50,
ResultEC[[i]] }],
{i,1,nElements}];(*グラフ表示のために適当な係数をかける*)
ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-1.199,2.909,1.245},
Lighting->False,PlotRange->All]
ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-1.312, 2.713, 1.540},
Lighting->False,PlotRange->All]
ParametricPlot3D[Evaluate[ResultEO],{sP,0,1},{tP,0,1},Compiled->False,
PlotPoints->6,FaceGrids->All,ViewPoint->{-0.017,9.000,-0.009},
Lighting->False,PlotRange->All]