here is the big orbitals and all there combinations, I also would like to
show you that $n is important( I noticed you took it out)
in that it can give you all the colors...
$Best,
-Mark Baker
use PDL;
use PDL::Complex;
use PDL::Graphics::TriD;
use PDL::NiceSlice;
use PDL::Math; keeptwiddling3d;
$PDL::BIGPDL=1;
## remember this can crash the computer if you dont have enough RAM and
swap in linux or pagefile in windows
for (;;){
$c = (40000) ; ## speed
for $phase(-360..360) { ## 360
$phases = $phase*0.0000001;
$frequency = ($c**-1);
$n= (6.28*$c*$frequency+($phases));
$photons = 4178; ### [Density] 2444 4178
$t=rvals(exp(zeros($photons)));
$cz=-1**$t*$c; # -1**$t*$c
$cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
$cx=-1**$c*sin(rvals($t))*$c; #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
$w=$cz-$cy-$cx;
$g=sin($w); #sin
$r=cos($cy+$c+$cz); #cos
$b=cos($w); #cos
$i=($cz-$cx-$cy);
$q=$c*($i*$n);
### altitude longitude latitude
$x = $b*sin($q); ### cos change the spin of the electron here
$y = $r*cos($q); ### sin
$z = $g*sin($q); ### cos
$rad=sum(sqrt(($x**2)+($y**2)+($z**2)));
$theta= sin($y*$x*$z)**-1;
$phi=cos($z*$rad)->rvals;
$x = $rad*cos($theta)*sin($phi); ### cos change the spin of the electron
here
$y = $rad*sin($theta)*sin($phi); ### sin
$z = $rad*cos($phi); ### cos
$xx = $b*sin($b); ### cos change the spin of the electron here
$yy = $r*($q); ### sin
$zz = $g*cos($g); ### cos
$dd++;
## torus
#$xxx = $xx*sin($yy)*cos($zz); ### cos change the spin of the electron
here
#$yyy = $xx*cos($yy)*cos($zz); ### sin
#$zzz = $xx*sin($zz); ### cos
# Draw a line
### d orbital i think
$xxxx = $xx*sin($yy)*cos($zz)+$xx; ### cos change the spin of the
electron here
$yyyy = $xx*cos($yy)*cos($zz)+$xx; ### sin
$zzzz = $xx*sin($zz); ### cos
# Draw a line
#### p orbital
$xxx = $xx*cos($yy)*sin($zz); ### cos change the spin of the electron
here
$yyy = $xx*sin($yy)*sin($zz); ### sin
$zzz = $xx*cos($zz);
$vv=sequence(4178);
### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
#$xxx=rotate($xxx,$phase);
#$yyy=rotate($yyy,$phase);
#$zzz=rotate($zzz,$phase);
#$xxx=$xxx*$vv;
#$yyy=$yyy*$vv;
#$zzz=$zzz*$vv;
#$xxxx=$xxxx*$vv;
#$yyyy=$yyyy*$vv;
#$zzzz=$zzzz*$vv;
$matrixa = ones(1000000);
$matrixb = ones(1000000);
$matrixc = ones(1000000);
$colora = ones(1000000);
$colorb = ones(1000000);
$colorc = ones(1000000);
#####################################################################
############# changing the colors can show the second electron in orbit
$matrixa->range(5000,[4178]) .=(($zzzz))-.2 ;
$matrixb->range(5000,[4178]) .=((-$xxxx))-.2 ;
$matrixc->range(5000,[4178]) .=(($yyyy))-.2 ;
$colora->range(5000,[4178]) .= $r+$n ;
$colorb->range(5000,[4178]) .= $g+$n ;
$colorc->range(5000,[4178]) .= $b+$n ;
$matrixa->range(10000,[4178]) .=((-$zzzz))-.2 ;
$matrixb->range(10000,[4178]) .=(($xxxx))-.2 ;
$matrixc->range(10000,[4178]) .=((-$yyyy))-.2 ;
$colora->range(10000,[4178]) .= $r+$n ;
$colorb->range(10000,[4178]) .= $g+$n ;
$colorc->range(10000,[4178]) .= $b+$n ;
$matrixa->range(15000,[4178]) .=((-$zzzz))-.2 ;
$matrixb->range(15000,[4178]) .=(($xxxx))-.2 ;
$matrixc->range(15000,[4178]) .=(($yyyy))-.2 ;
$colora->range(15000,[4178]) .= $r-$n ;
$colorb->range(15000,[4178]) .= $g-$n ;
$colorc->range(15000,[4178]) .= $b+$n ;
$matrixa->range(20000,[4178]) .=(($zzzz))-.2 ;
$matrixb->range(20000,[4178]) .=((-$xxxx))-.2 ;
$matrixc->range(20000,[4178]) .=((-$yyyy))-.2 ;
$colora->range(20000,[4178]) .= $r-$n ;
$colorb->range(20000,[4178]) .= $g-$n ;
$colorc->range(20000,[4178]) .= $b+$n ;
$matrixa->range(25000,[4178]) .=((-$yyyy))-.2 ;
$matrixb->range(25000,[4178]) .=(($zzzz))-.2 ;
$matrixc->range(25000,[4178]) .=(($xxxx))-.2 ;
$colora->range(25000,[4178]) .= $r+$n ;
$colorb->range(25000,[4178]) .= $g-$n ;
$colorc->range(25000,[4178]) .= $b-$n ;
$matrixa->range(30000,[4178]) .=(($yyyy))-.2 ;
$matrixb->range(30000,[4178]) .=((-$zzzz))-.2 ;
$matrixc->range(30000,[4178]) .=((-$xxxx))-.2 ;
$colora->range(30000,[4178]) .= $r+$n ;
$colorb->range(30000,[4178]) .= $g-$n ;
$colorc->range(30000,[4178]) .= $b-$n ;
$matrixa->range(35000,[4178]) .=(($xxxx))-.2 ;
$matrixb->range(35000,[4178]) .=((-$yyyy))-.2 ;
$matrixc->range(35000,[4178]) .=((-$zzzz))-.2 ;
$colora->range(35000,[4178]) .= $r-$n ;
$colorb->range(35000,[4178]) .= $g+$n ;
$colorc->range(35000,[4178]) .= $b+$n ;
$matrixa->range(40000,[4178]) .=((-$xxxx))-.2 ;
$matrixb->range(40000,[4178]) .=(($yyyy))-.2 ;
$matrixc->range(40000,[4178]) .=(($zzzz))-.2 ;
$colora->range(40000,[4178]) .= $r-$n ;
$colorb->range(40000,[4178]) .= $g+$n ;
$colorc->range(40000,[4178]) .= $b+$n ;
$matrixa->range(45000,[4178]) .=(($xxxx))-.2 ;
$matrixb->range(45000,[4178]) .=(($yyyy))-.2 ;
$matrixc->range(45000,[4178]) .=(($zzzz))-.2 ;
$colora->range(45000,[4178]) .= $r+$n ;
$colorb->range(45000,[4178]) .= $g+$n ;
$colorc->range(45000,[4178]) .= $b-$n ;
$matrixa->range(50000,[4178]) .=((-$xxxx))-.2 ;
$matrixb->range(50000,[4178]) .=((-$yyyy))-.2 ;
$matrixc->range(50000,[4178]) .=((-$zzzz))-.2 ;
$colora->range(50000,[4178]) .= $r+$n ;
$colorb->range(50000,[4178]) .= $g+$n ;
$colorc->range(50000,[4178]) .= $b-$n ;
$matrixa->range(55000,[4178]) .=(($zzzz))-.2 ;
$matrixb->range(55000,[4178]) .=(($xxxx))-.2 ;
$matrixc->range(55000,[4178]) .=(($yyyy))-.2 ;
$colora->range(55000,[4178]) .= $r+$n ;
$colorb->range(55000,[4178]) .= $g-$n ;
$colorc->range(55000,[4178]) .= $b+$n ;
$matrixa->range(60000,[4178]) .=((-$zzzz))-.2 ;
$matrixb->range(60000,[4178]) .=((-$xxxx))-.2 ;
$matrixc->range(60000,[4178]) .=((-$yyyy))-.2 ;
$colora->range(60000,[4178]) .= $r+$n ;
$colorb->range(60000,[4178]) .= $g-$n ;
$colorc->range(60000,[4178]) .= $b+$n ;
#
$matrixa->range(65000,[4178]) .=(($yyyy))-.2 ;
$matrixb->range(65000,[4178]) .=(($zzzz))-.2 ;
$matrixc->range(65000,[4178]) .=(($xxxx))-.2 ;
$colora->range(65000,[4178]) .= $r-$n ;
$colorb->range(65000,[4178]) .= $g+$n ;
$colorc->range(65000,[4178]) .= $b-$n ;
$matrixa->range(70000,[4178]) .=((-$yyyy))-.2 ;
$matrixb->range(70000,[4178]) .=((-$zzzz))-.2 ;
$matrixc->range(70000,[4178]) .=((-$xxxx))-.2 ;
$colora->range(70000,[4178]) .= $r-$n ;
$colorb->range(70000,[4178]) .= $g+$n ;
$colorc->range(70000,[4178]) .= $b-$n ;
# $matrixa->range(75000,[4178]) .=.0009*($x)+(-.2) ;
# $matrixb->range(75000,[4178]) .=.0009*($y)+(-.2) ;
# $matrixc->range(75000,[4178]) .=.0009*($z)+(-.2) ;
# $colora->range(75000,[4178]) .= $r ;
# $colorb->range(75000,[4178]) .= $g ;
# $colorc->range(75000,[4178]) .= $b ;
# $matrixa=rotate($matrixa,$dd);
#$matrixb=rotate($matrixb,$dd);
#$matrixc=rotate($matrixc,$dd);
+points3d[$matrixa+=$matrixa,$matrixb+=$matrixb,$matrixc+=$matrixc],
[$colora,$colorb,$colorc];#{PointSize=>1};
}} #}
On Sun, Oct 15, 2023 at 9:15 PM Mark NanoNebulas <nanonebu...@gmail.com>
wrote:
> i found the torus matrix here
>
> $xxx = $xx*sin($yy)*cos($zz);
> $yyy = $xx*cos($yy)*cos($zz);
> $zzz = $xx*sin($zz);
>
> and if you add the $xx to $xxx and $yyy you get bigger orbs
> like
>
> $xxx = $xx*sin($yy)*cos($zz)+$xx;
> $yyy = $xx*cos($yy)*cos($zz)+$xx;
> $zzz = $xx*sin($zz);
>
> would that be all the orbitals i need ???
>
> Best,
>
> -Mark Baker
>
>
>
>
>
>
> On Sun, Oct 15, 2023 at 5:46 PM Mark NanoNebulas <nanonebu...@gmail.com>
> wrote:
>
>> the s orbital is x y z the p orbital is xxxx yyyy zzzz the d orbitals are
>> xxx yyy zzz
>> is that right and can they be used to describe all atoms ???
>>
>> -Mark Baker
>>
>> On Sun, Oct 15, 2023 at 5:26 PM Mark NanoNebulas <nanonebu...@gmail.com>
>> wrote:
>>
>>> have you used the range function ??? here i have all orbitals that can
>>> be without using hold3d();
>>>
>>> could you put the Rotxz into 4178 dimensions I dont understand how your
>>> doing it in 3
>>> here are all the orbitals i got 14 of them total in your opinion can i
>>> use these to describe all the atoms
>>>
>>> check out the range function it lets you use a lot of piddles in one
>>> matrix...
>>>
>>> Best,
>>> -Mark
>>> #####################
>>>
>>> use PDL;
>>> use PDL::Complex;
>>> use PDL::Graphics::TriD;
>>> use PDL::NiceSlice;
>>> use PDL::Math; keeptwiddling3d;
>>> $PDL::BIGPDL=1;
>>> ## remember this can crash the computer if you dont have enough RAM and
>>> swap in linux or pagefile in windows
>>>
>>>
>>>
>>> for (;;){
>>> $c = (4000) ; ## speed
>>> for $phase(-360..360) { ## 360
>>> $phases = $phase*0.0000001;
>>> $frequency = ($c**-1);
>>> $n= (6.28*$c*$frequency+($phases));
>>>
>>> $photons = 4178; ### [Density] 2444 4178
>>>
>>>
>>> $t=rvals(exp(zeros($photons)));
>>> $cz=-1**$t*$c; # -1**$t*$c
>>> $cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
>>> $cx=-1**$c*sin(rvals($t))*$c;
>>> #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
>>>
>>>
>>> $w=$cz-$cy-$cx;
>>> $g=sin($w); #sin
>>> $r=cos($cy+$c+$cz); #cos
>>> $b=cos($w); #cos
>>> $i=($cz-$cx-$cy);
>>> $q=$c*($i*$n);
>>>
>>> ### altitude longitude latitude
>>> $x = $b*sin($q); ### cos change the spin of the electron here
>>> $y = $r*cos($q); ### sin
>>> $z = $g*sin($q); ### cos
>>>
>>>
>>> $xx = $b*sin($b); ### cos change the spin of the electron here
>>> $yy = $r*($q); ### sin
>>> $zz = $g*cos($g); ### cos
>>> $dd++;
>>> $xxx = $xx*cos($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### cos change the
>>> spin of the electron here
>>> $yyy = $xx*sin($yy)*sin($zz)+sin($xx*(1/sqrt 2)); ### sin
>>> $zzz = $xx*cos($zz); ### cos
>>> # Draw a line
>>>
>>> $xxxx = $xx*cos($yy)*sin($zz); ### cos change the spin of the
>>> electron here
>>> $yyyy = $xx*sin($yy)*sin($zz); ### sin
>>> $zzzz = $xx*cos($zz);
>>>
>>>
>>>
>>> $vv=sequence(4178);
>>> ### wavefunction #exp(-sqrt(($x**2)+($y**2)+($z**2)))
>>>
>>> #$xxx=rotate($xxx,$phase);
>>> #$yyy=rotate($yyy,$phase);
>>> #$zzz=rotate($zzz,$phase);
>>>
>>> #$xxx=$xxx*$vv;
>>> #$yyy=$yyy*$vv;
>>> #$zzz=$zzz*$vv;
>>>
>>> #$xxxx=$xxxx*$vv;
>>> #$yyyy=$yyyy*$vv;
>>> #$zzzz=$zzzz*$vv;
>>>
>>>
>>>
>>>
>>> $matrixa = ones(1000000);
>>> $matrixb = ones(1000000);
>>> $matrixc = ones(1000000);
>>> $colora = ones(1000000);
>>> $colorb = ones(1000000);
>>> $colorc = ones(1000000);
>>>
>>> $matrixa->range(5000,[4178]) .=0.8*(($xxx))-.2 ;
>>> $matrixb->range(5000,[4178]) .=0.8*(($yyy))-.2 ;
>>> $matrixc->range(5000,[4178]) .=0.8*(($zzz))-.2 ;
>>> $colora->range(5000,[4178]) .= $r ;
>>> $colorb->range(5000,[4178]) .= $g ;
>>> $colorc->range(5000,[4178]) .= $b ;
>>>
>>> $matrixa->range(10000,[4178]) .=0.8*((-$xxx))-.2 ;
>>> $matrixb->range(10000,[4178]) .=0.8*((-$yyy))-.2 ;
>>> $matrixc->range(10000,[4178]) .=0.8*((-$zzz))-.2 ;
>>> $colora->range(10000,[4178]) .= $r ;
>>> $colorb->range(10000,[4178]) .= $g ;
>>> $colorc->range(10000,[4178]) .= $b ;
>>>
>>> $matrixa->range(15000,[4178]) .=0.8*((-$xxx))-.2 ;
>>> $matrixb->range(15000,[4178]) .=0.8*((-$yyy))-.2 ;
>>> $matrixc->range(15000,[4178]) .=0.8*(($zzz))-.2 ;
>>> $colora->range(15000,[4178]) .= $r ;
>>> $colorb->range(15000,[4178]) .= $g ;
>>> $colorc->range(15000,[4178]) .= $b ;
>>>
>>> $matrixa->range(20000,[4178]) .=0.8*(($xxx))-.2 ;
>>> $matrixb->range(20000,[4178]) .=0.8*((-$yyy))-.2 ;
>>> $matrixc->range(20000,[4178]) .=0.8*(($zzz))-.2 ;
>>> $colora->range(20000,[4178]) .= $r ;
>>> $colorb->range(20000,[4178]) .= $g ;
>>> $colorc->range(20000,[4178]) .= $b ;
>>>
>>> $matrixa->range(25000,[4178]) .=0.8*((-$xxx))-.2 ;
>>> $matrixb->range(25000,[4178]) .=0.8*(($yyy))-.2 ;
>>> $matrixc->range(25000,[4178]) .=0.8*((-$zzz))-.2 ;
>>> $colora->range(25000,[4178]) .= $r ;
>>> $colorb->range(25000,[4178]) .= $g ;
>>> $colorc->range(25000,[4178]) .= $b ;
>>>
>>> $matrixa->range(30000,[4178]) .=0.8*(($xxx))-.2 ;
>>> $matrixb->range(30000,[4178]) .=0.8*((-$yyy))-.2 ;
>>> $matrixc->range(30000,[4178]) .=0.8*((-$zzz))-.2 ;
>>> $colora->range(30000,[4178]) .= $r ;
>>> $colorb->range(30000,[4178]) .= $g ;
>>> $colorc->range(30000,[4178]) .= $b ;
>>>
>>> $matrixa->range(35000,[4178]) .=0.8*(($xxx))-.2 ;
>>> $matrixb->range(35000,[4178]) .=0.8*(($yyy))-.2 ;
>>> $matrixc->range(35000,[4178]) .=0.8*((-$zzz))-.2 ;
>>> $colora->range(35000,[4178]) .= $r ;
>>> $colorb->range(35000,[4178]) .= $g ;
>>> $colorc->range(35000,[4178]) .= $b ;
>>>
>>> $matrixa->range(40000,[4178]) .=0.8*((-$xxx))-.2 ;
>>> $matrixb->range(40000,[4178]) .=0.8*(($yyy))-.2 ;
>>> $matrixc->range(40000,[4178]) .=0.8*(($zzz))-.2 ;
>>> $colora->range(40000,[4178]) .= $r ;
>>> $colorb->range(40000,[4178]) .= $g ;
>>> $colorc->range(40000,[4178]) .= $b ;
>>>
>>> $matrixa->range(45000,[4178]) .=(($xxxx))-.2 ;
>>> $matrixb->range(45000,[4178]) .=(($yyyy))-.2 ;
>>> $matrixc->range(45000,[4178]) .=(($zzzz))-.2 ;
>>> $colora->range(45000,[4178]) .= $r ;
>>> $colorb->range(45000,[4178]) .= $g ;
>>> $colorc->range(45000,[4178]) .= $b ;
>>>
>>> $matrixa->range(50000,[4178]) .=((-$xxxx))-.2 ;
>>> $matrixb->range(50000,[4178]) .=((-$yyyy))-.2 ;
>>> $matrixc->range(50000,[4178]) .=((-$zzzz))-.2 ;
>>> $colora->range(50000,[4178]) .= $r ;
>>> $colorb->range(50000,[4178]) .= $g ;
>>> $colorc->range(50000,[4178]) .= $b ;
>>>
>>> $matrixa->range(55000,[4178]) .=(($zzzz))-.2 ;
>>> $matrixb->range(55000,[4178]) .=(($xxxx))-.2 ;
>>> $matrixc->range(55000,[4178]) .=(($yyyy))-.2 ;
>>> $colora->range(55000,[4178]) .= $r ;
>>> $colorb->range(55000,[4178]) .= $g ;
>>> $colorc->range(55000,[4178]) .= $b ;
>>>
>>> $matrixa->range(60000,[4178]) .=((-$zzzz))-.2 ;
>>> $matrixb->range(60000,[4178]) .=((-$xxxx))-.2 ;
>>> $matrixc->range(60000,[4178]) .=((-$yyyy))-.2 ;
>>> $colora->range(60000,[4178]) .= $r ;
>>> $colorb->range(60000,[4178]) .= $g ;
>>> $colorc->range(60000,[4178]) .= $b ;
>>> #
>>> $matrixa->range(65000,[4178]) .=(($yyyy))-.2 ;
>>> $matrixb->range(65000,[4178]) .=(($zzzz))-.2 ;
>>> $matrixc->range(65000,[4178]) .=(($xxxx))-.2 ;
>>> $colora->range(65000,[4178]) .= $r ;
>>> $colorb->range(65000,[4178]) .= $g ;
>>> $colorc->range(65000,[4178]) .= $b ;
>>>
>>> $matrixa->range(70000,[4178]) .=((-$yyyy))-.2 ;
>>> $matrixb->range(70000,[4178]) .=((-$zzzz))-.2 ;
>>> $matrixc->range(70000,[4178]) .=((-$xxxx))-.2 ;
>>> $colora->range(70000,[4178]) .= $r ;
>>> $colorb->range(70000,[4178]) .= $g ;
>>> $colorc->range(70000,[4178]) .= $b ;
>>>
>>> # $matrixa=rotate($matrixa,$dd);
>>> #$matrixb=rotate($matrixb,$dd);
>>> #$matrixc=rotate($matrixc,$dd);
>>>
>>>
>>>
>>> +points3d[$matrixa,$matrixb,$matrixc],
>>> [$colora,$colorb,$colorc];#{PointSize=>1};
>>>
>>>
>>> }} #}
>>>
>>>
>>>
>>> On Sat, Oct 14, 2023 at 10:05 PM Luis Mochan <moc...@icf.unam.mx> wrote:
>>>
>>>> This worked for me to generate x-z d orbitals. The other planes are
>>>> similar.
>>>> I simply applied 4 different 45 degree rotations to your points and
>>>> plotted them all together. I added two additional points to fix the
>>>> scale.
>>>>
>>>> ##################################################
>>>> use PDL;
>>>> use PDL::Complex;
>>>> use PDL::Graphics::TriD;
>>>> use PDL::Math;
>>>> my $c = 10; ## speed
>>>> my $Rotxz=pdl([[1/sqrt 2, 0, -1/sqrt 2],[0,1,0], [1/sqrt 2, 0, 1/sqrt
>>>> 2]]);
>>>> nokeeptwiddling3d;
>>>> while(1){
>>>> for my $phase(-360..360) { ## 360
>>>> my $photons = 4178; ### [Density] 2444 4178
>>>> my $t=rvals(zeros($photons));
>>>> my $cz=$c; # -1**$t*$c
>>>> my $cy=sin($t*$c);
>>>> my $cx=sin($t)*$c;
>>>> my $w=$cz-$cy-$cx;
>>>> my $color=(pdl[my $r=cos($cy+$c+$cz), my $g=sin($w), my
>>>> $b=cos($w)])->mv(-1,0);
>>>> my $i=($cz-$cx-$cy);
>>>> my $q=$c*($i*$phase);
>>>> my $R = $b*sin($b);
>>>> my $phi = $r*$q;
>>>> my $theta = $g*cos($g);
>>>> my $points=($R*pdl[sin($theta)*cos($phi),
>>>> sin($phi)*sin($theta),cos($theta)]);
>>>> my $allpoints=pdl(
>>>> $Rotxz x $points, $Rotxz->transpose x $points, -$Rotxz x
>>>> $points, -$Rotxz->transpose x $points
>>>> )->mv(1,0);
>>>> points3d(pdl[[1,1,1],[-1,-1,-1]]);
>>>> hold3d();
>>>> points3d($allpoints, $color->dummy(2));
>>>> release3d();
>>>> }
>>>> }
>>>>
>>>>
>>>> On Sat, Oct 14, 2023 at 06:01:34PM -0700, Mark NanoNebulas wrote:
>>>> > i found out how to move the point in the y axis here
>>>> > $dd++;
>>>> > $xxx = $xx*cos($yy)*sin($zz)+sin($xx); ### cos change the spin of
>>>> the
>>>> > electron here
>>>> > $yyy = $xx*sin($yy)*sin($zz)+sin($xx+$dd); ### sin
>>>> > $zzz = $xx*cos($zz); ### cos
>>>> >
>>>> > this seems to work but can it make all the d orbitals i need ???
>>>> >
>>>> > -Mark Baker
>>>> >
>>>> >
>>>> > On Sat, Oct 14, 2023 at 7:18 AM Luis Mochan <moc...@icf.unam.mx>
>>>> wrote:
>>>> >
>>>> > > Hi,
>>>> > > I didn't quite understand the code. I simplified it a bit, removed
>>>> > > many noop instructions, changed the speed, negated the twiddling and
>>>> > > the result looks like a nice lobe of a p orbital, made up of points
>>>> > > along some spiral, rotating about its axis. Is this the rotation you
>>>> > > want or do you want to rotate the symmetry axis from the z
>>>> direction to an
>>>> > > arbitrary direction? You could apply a 3x3 rotation matrix to the
>>>> > > $points array in the code below. If you want a point to remain
>>>> fixed,
>>>> > > first translate it to the origin, the rotate, and then rotate back.
>>>> > > Regards,
>>>> > > Luis
>>>> > > ----------------
>>>> > > ##################################################
>>>> > > use PDL;
>>>> > > use PDL::Complex;
>>>> > > use PDL::Graphics::TriD;
>>>> > > use PDL::Math;
>>>> > > nokeeptwiddling3d;
>>>> > > my $c = 10; ## speed
>>>> > > while(1){
>>>> > > for my $phase(-360..360) { ## 360
>>>> > > my $photons = 4178; ### [Density] 2444 4178
>>>> > > my $t=rvals(zeros($photons));
>>>> > > my $cz=$c; # -1**$t*$c
>>>> > > my $cy=sin($t*$c);
>>>> > > my $cx=sin($t)*$c;
>>>> > > my $w=$cz-$cy-$cx;
>>>> > > my $color=(pdl[my $r=cos($cy+$c+$cz), my $g=sin($w), my
>>>> > > $b=cos($w)])->mv(-1,0);
>>>> > > my $i=($cz-$cx-$cy);
>>>> > > my $q=$c*($i*$phase);
>>>> > > my $R = $b*sin($b);
>>>> > > my $phi = $r*$q;
>>>> > > my $theta = $g*cos($g);
>>>> > > my $points=($R*pdl[sin($theta)*cos($phi),
>>>> > > sin($phi)*sin($theta),cos($theta)])->mv(-1,0);
>>>> > > points3d $points,$color;
>>>> > > }
>>>> > > }
>>>> > >
>>>> > >
>>>> > >
>>>> > > On Fri, Oct 13, 2023 at 10:36:17AM -0700, Mark NanoNebulas wrote:
>>>> > > > Hi im trying to rotate this obital so that i can place them all
>>>> over in
>>>> > > > balanced arrangements of orbs like the chemistry book says but i
>>>> dont
>>>> > > know
>>>> > > > how to rotate a matrix in all directions, can anybody help
>>>> > > > heres the code
>>>> > > >
>>>> > > > ##################################################
>>>> > > > use PDL;
>>>> > > > use PDL::Complex;
>>>> > > > use PDL::Graphics::TriD;
>>>> > > > use PDL::Math; keeptwiddling3d;
>>>> > > > $PDL::BIGPDL=1;
>>>> > > > ## remember this can crash the computer if you dont have enough
>>>> RAM and
>>>> > > > swap in linux or pagefile in windows
>>>> > > >
>>>> > > >
>>>> > > >
>>>> > > > for (;;){
>>>> > > > $c = (400) ; ## speed
>>>> > > > for $phase(-360..360) { ## 360
>>>> > > > $phases = $phase*0.0000001;
>>>> > > > $frequency = ($c**-1);
>>>> > > > $n= (6.28*$c*$frequency+($phases));
>>>> > > >
>>>> > > > $photons = 4178; ### [Density] 2444 4178
>>>> > > >
>>>> > > >
>>>> > > > $t=rvals(exp(zeros($photons)));
>>>> > > > $cz=-1**$t*$c; # -1**$t*$c
>>>> > > > $cy=-1**sin($t*$c); #-1**$c*sin($t*$c)
>>>> > > > $cx=-1**$c*sin(rvals($t))*$c;
>>>> > > #-1**$c*bessj0(rvals($t,{Centre=>[0]}))*$c;
>>>> > > >
>>>> > > >
>>>> > > > $w=$cz-$cy-$cx;
>>>> > > > $g=sin($w); #sin
>>>> > > > $r=cos($cy+$c+$cz); #cos
>>>> > > > $b=cos($w); #cos
>>>> > > > $i=($cz-$cx-$cy);
>>>> > > > $q=$c*($i*$n);
>>>> > > > $xx = $b*sin($b); ###
>>>> > > > $yy = $r*($q); ###
>>>> > > > $zz = $g*cos($g); ###
>>>> > > >
>>>> > > > $xxx = $xx*cos($yy)*sin($zz); ###
>>>> > > > $yyy = $xx*sin($yy)*sin($zz); ###
>>>> > > > $zzz = $xx*cos($zz);
>>>> > > >
>>>> > > > points3d[$xxx,$yyy,$zzz],[$r,$g,$b];
>>>> > > >
>>>> > > > ##########################################################
>>>> > > >
>>>> > > > if i can get this orbital to rotate in all directions, then i can
>>>> > > represent
>>>> > > > atoms
>>>> > > > better cause they are balanced by the number of electrons so the
>>>> orbs
>>>> > > > repeal each other, please any help would be appreciated !!!
>>>> > > >
>>>> > > > Best,
>>>> > > > -Mark Baker
>>>> > >
>>>> > >
>>>> > > > _______________________________________________
>>>> > > > pdl-general mailing list
>>>> > > > pdl-general@lists.sourceforge.net
>>>> > > > https://lists.sourceforge.net/lists/listinfo/pdl-general
>>>> > >
>>>> > >
>>>> > > --
>>>> > >
>>>> > > o
>>>> > > W. Luis Mochán, | tel:(52)(777)329-1734
>>>> /<(*)
>>>> > > Instituto de Ciencias Físicas, UNAM | fax:(52)(777)317-5388
>>>> `>/ /\
>>>> > > Av. Universidad s/n CP 62210 |
>>>> (*)/\/ \
>>>> > > Cuernavaca, Morelos, México | moc...@fis.unam.mx
>>>> /\_/\__/
>>>> > > GPG: 791EB9EB, C949 3F81 6D9B 1191 9A16 C2DF 5F0A C52B 791E B9EB
>>>> > >
>>>> > >
>>>> > > _______________________________________________
>>>> > > pdl-general mailing list
>>>> > > pdl-general@lists.sourceforge.net
>>>> > > https://lists.sourceforge.net/lists/listinfo/pdl-general
>>>> > >
>>>>
>>>>
>>>> > _______________________________________________
>>>> > pdl-general mailing list
>>>> > pdl-general@lists.sourceforge.net
>>>> > https://lists.sourceforge.net/lists/listinfo/pdl-general
>>>>
>>>>
>>>> --
>>>>
>>>> o
>>>> W. Luis Mochán, | tel:(52)(777)329-1734 /<(*)
>>>> Instituto de Ciencias Físicas, UNAM | fax:(52)(777)317-5388 `>/
>>>> /\
>>>> Av. Universidad s/n CP 62210 |
>>>> (*)/\/ \
>>>> Cuernavaca, Morelos, México | moc...@fis.unam.mx /\_/\__/
>>>> GPG: 791EB9EB, C949 3F81 6D9B 1191 9A16 C2DF 5F0A C52B 791E B9EB
>>>>
>>>
_______________________________________________
pdl-general mailing list
pdl-general@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pdl-general
