Logo by kameelian - Contribute your own Logo!

END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

it was a great time but no longer maintainable by c.Kleinhuis contact him for any data retrieval,
thanks and see you perhaps in 10 years again

this forum will stay online for reference
News: Support us via Flattr FLATTR Link
 
*
Welcome, Guest. Please login or register. April 20, 2024, 02:25:21 AM


Login with username, password and session length


The All New FractalForums is now in Public Beta Testing! Visit FractalForums.org and check it out!


Pages: [1]   Go Down
  Print  
Share this topic on DiggShare this topic on FacebookShare this topic on GoogleShare this topic on RedditShare this topic on StumbleUponShare this topic on Twitter
Author Topic: Nested Divergence, Rotation, Nabla  (Read 1095 times)
0 Members and 1 Guest are viewing this topic.
kram1032
Fractal Senior
******
Posts: 1863


« on: December 12, 2011, 10:07:14 PM »

Hi there. LOOONG time no see. smiley

Now that I started studying, I've learned a thing or two about multidimensional derivatives.
The case where you take a generic differentiable nested function's derivative is quite easy and was found by one of you guys quite a time ago.
Back then I unnecessarily started a thread of that topic, deriving the result myself.

It basically was, as the chain rule quickly shows, f'(x) f'(f(x)) f'(f(f(x))) f'(...

Now I wonder what happens if you take the complex numbers as vector.
E.g what's Nabla*(f(...,...)+i g(...,...))

By quickly trying it out I noticed that in case of quaternions, the real part of such an operation ends up looking a lot like divergence and the imaginary part appears to be the curl.
This I was able to verify. Wikipedia actually has that listed in the article about Quaternions.
However, in that case, the original quaternions where "pure", so they didn't have any real part.
(0,1/dx,1/dy,1/dz)(0,f(x,y,z),g(x,y,z),h(x,y,z))=div(f,g,h)+I curl(f,g,h)

So the complex case, not being useful when getting rid of the real part, should behave a little differently.
Same goes for a full quaternion equation set.

Now the problem is, and I checked on it with Mathematica which I now have, what's the pattern of a nested two argument function (alternatively four arguments for quaternions) ?
There obviously is one and Mathematica doesn't have a problem with finding that but I failed to derive a simple running derivative type way of calculating  that.

I'd be interested in
1) how does the shape directly change?
2) could this be used for any kind of interesting colouring?

I did do a search but nothing like this popped up accordingly so I hope this isn't as pointless as the scalar running derivative turned out to be.

In case nothing interesting shows up by doing so, maybe doing a transpose multiplication returns something more interesting...
« Last Edit: December 12, 2011, 11:45:56 PM by kram1032 » Logged
kram1032
Fractal Senior
******
Posts: 1863


« Reply #1 on: December 13, 2011, 03:17:32 PM »

That would be an example for a derivative with respect to x, of two two-argument functions, nested four times.
Any obvious patterns for the Eagle-eyed programmer wizzes of you? smiley

\begin{doublespace}
\noindent\(\(f^{(1,0)}(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))) \\
\left(f^{(1,0)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))) \right.\\
\left(g^{(1,0)}(x,y) f^{(0,1)}(f(x,y),g(x,y))+f^{(1,0)}(x,y) f^{(1,0)}(f(x,y),g(x,y))\right)+\left.f^{(0,1)}\right(f(f(x,y),g(x,y)),\\
\left.g(f(x,y),g(x,y))) \left(f^{(1,0)}(x,y) g^{(1,0)}(f(x,y),g(x,y))+g^{(1,0)}(x,y) g^{(0,1)}(f(x,y),g(x,y))\right)\right)+\\
f^{(0,1)}(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))) \\
\left(g^{(0,1)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))) \right.\\
\left(f^{(1,0)}(x,y) g^{(1,0)}(f(x,y),g(x,y))+g^{(1,0)}(x,y) g^{(0,1)}(f(x,y),g(x,y))\right)+\\
\left(g^{(1,0)}(x,y) f^{(0,1)}(f(x,y),g(x,y))+f^{(1,0)}(x,y) f^{(1,0)}(f(x,y),g(x,y))\right) \\
\left.g^{(1,0)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))\right)\)\)
\end{doublespace}



Edit: Does the Latex work right now? Seems to be only an unknown image on my end... Writing it all by hand is pretty... long and error prone, so I'd like to avoid that <.<

Edit: Here is the actual final result to be viewed. It's the Nabla operator on a complex function, viewed as vector after four iterations:
\begin{doublespace}
\noindent\(\(\left(f^{(1,0)}(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))) \right.\\
\left(f^{(1,0)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))) \left(\left(g^{(1,0)}(x,y)+i g^{(0,1)}(x,y)\right) f^{(0,1)}(f(x,y),g(x,y))+\left(f^{(1,0)}(x,y)+i
f^{(0,1)}(x,y)\right) f^{(1,0)}(f(x,y),g(x,y))\right)+\right.\\
\left.f^{(0,1)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))) \left(\left(f^{(1,0)}(x,y)+i f^{(0,1)}(x,y)\right) g^{(1,0)}(f(x,y),g(x,y))+\left(g^{(1,0)}(x,y)+i
g^{(0,1)}(x,y)\right) g^{(0,1)}(f(x,y),g(x,y))\right)\right)+\left.f^{(0,1)}\right(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),\\
g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))) \left(g^{(0,1)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))) \left(\left(f^{(1,0)}(x,y)+i f^{(0,1)}(x,y)\right) g^{(1,0)}(f(x,y),g(x,y))+\left(g^{(1,0)}(x,y)+i
g^{(0,1)}(x,y)\right) g^{(0,1)}(f(x,y),g(x,y))\right)+\right.\\
\left.\left.\left(\left(g^{(1,0)}(x,y)+i g^{(0,1)}(x,y)\right) f^{(0,1)}(f(x,y),g(x,y))+\left(f^{(1,0)}(x,y)+i f^{(0,1)}(x,y)\right) f^{(1,0)}(f(x,y),g(x,y))\right)
g^{(1,0)}(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))\right)\right) \\
\left(f'(f(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))))+i g'(f(f(f(f(x,y),g(x,y)),g(f(x,y),g(x,y))),g(f(f(x,y),g(x,y)),g(f(x,y),g(x,y)))))\right)\)\)
\end{doublespace}


« Last Edit: December 13, 2011, 04:51:57 PM by kram1032 » Logged
cKleinhuis
Administrator
Fractal Senior
*******
Posts: 7044


formerly known as 'Trifox'


WWW
« Reply #2 on: December 13, 2011, 03:56:30 PM »

ehrm, kram, you need to put [latex] around your latex ...
Logged

---

divide and conquer - iterate and rule - chaos is No random!
cKleinhuis
Administrator
Fractal Senior
*******
Posts: 7044


formerly known as 'Trifox'


WWW
« Reply #3 on: December 13, 2011, 03:57:56 PM »

there seems to be a problem, my test: [latex]x^2[/latex]
Logged

---

divide and conquer - iterate and rule - chaos is No random!
cKleinhuis
Administrator
Fractal Senior
*******
Posts: 7044


formerly known as 'Trifox'


WWW
« Reply #4 on: December 13, 2011, 04:05:39 PM »

damn, dunno what going wrong, i just updated the late plugin, but somehow it doesnt work, neither i got any error messages ... shit!
Logged

---

divide and conquer - iterate and rule - chaos is No random!
cKleinhuis
Administrator
Fractal Senior
*******
Posts: 7044


formerly known as 'Trifox'


WWW
« Reply #5 on: December 13, 2011, 04:30:39 PM »

\lim_{n \to \infty}\sum_{k=1}^n \frac{1}{k^2}= \frac{\pi^2}{6}



try using "jstex" as tag, it is behaving very weird right now, i installed a js latex parser plugin....
Logged

---

divide and conquer - iterate and rule - chaos is No random!
kram1032
Fractal Senior
******
Posts: 1863


« Reply #6 on: December 13, 2011, 04:57:50 PM »

Well, as you can see, that worked in the sense that it now shows anything at all. Though it does NOT format the text :-/
Any obvious flaw? I directly used the Mathematica Output, saved as a .tex file.

Btw, if we're at it, the symbol in the editor is broken too.

Edit: Oh, ok... eqn too long... interesting... <.<
Great. I gotta write it by hand when I have time...
Logged
kram1032
Fractal Senior
******
Posts: 1863


« Reply #7 on: December 17, 2011, 02:01:13 PM »

Here is the derivative for a single iteration. It should be some kind of simple pattern that most likely has a nice loop-like algorithmic form for additional iterations. But it quickly becomes so cluttered that it's hard to tell.
I guess this kind of stuff is at the edge of the usefulness of this kind of notation. Maybe we could develop some kind of nice IFS notation that makes it easy to grasp these kinds of patterns and possibly other things in the theory of fractals in general...

df(f,g)/dx=f(1,0)(x,y) f(1,0)(f(x,y),g(x,y))+f(0,1)(f(x,y),g(x,y)) g(1,0)(x,y)
 df(f,g)/dy=f(0,1)(f(x,y),g(x,y)) g(0,1)(x,y)+f(0,1)(x,y) f(1,0)(f(x,y),g(x,y))
 dg(f,g)/dx=g(0,1)(f(x,y),g(x,y)) g(1,0)(x,y)+f(1,0)(x,y) g(1,0)(f(x,y),g(x,y))
 dg(f,g)/dy=g(0,1)(x,y) g(0,1)(f(x,y),g(x,y))+f(0,1)(x,y) g(1,0)(f(x,y),g(x,y))

Let's see...
The first line roughly says:
d/dx f(0,x,y) * d/dx f(1,x,y) + d/dx g(0,x,y) d/dy f(1,x,y)

The first half of that looks a lot like the expansion we already have for the scalar running derivative. So I guess, this extends to:
f'(x) f'(f(x)) f'(f(f(x))) ...

The second half looks similar, except that at first you have the derivative of the second function, followed by the derivative of the first one by the other variable. A single iteration is too little to find a pattern in that... So let's look at a deeper nesting level...

An additional iteration gives:

d/dx f(0,x,y) d/dx f(1,x,y) d/dx f(2,x,y) +
d/dx g(0,x,y) d/dy g(1,x,y) d/dy f(2,x,y) +
d/dx g(0,x,y) d/dy f(1,x,y) d/dx f(2,x,y) +
d/dx f(0,x,y) d/dx g(1,x,y) d/dy f(2,x,y)

The first number represents the iteration.
d/dx here isn't used in the usual sense. It's just to see, which variable is derived.

d/dx f(0,x,y) here would mean
f(1,0)(x,y)
or alternatively
f(x)(x,y)

I'm detecting some kind of binomial/pascal-triangle expansion/permutation thing.

fdx->fdx->fdx is what we already found for scalar derivatives.
Here we have three more terms in case of the second iteration... and I have the feeling that they'll EXPlode for higher iterations, as permutations tend to do.
fdx -> fdx -> fdx
gdx -> gdy -> fdy
gdx -> fdy -> fdx
fdx -> gdx -> fdy

aaa
cdb
cba
acb

Still not a clear pattern...
The last two lines are anti-cyclic.
If the second one was bac, it would fit better into the picture of this.
« Last Edit: December 17, 2011, 03:25:42 PM by kram1032 » Logged
Pages: [1]   Go Down
  Print  
 
Jump to:  

Related Topics
Subject Started by Replies Views Last post
Rotation Animation Let's collaborate on something! « 1 2 » Wel lEnTaoed 15 4036 Last post November 22, 2010, 01:53:15 PM
by Bent-Winged Angel
EJS outer-para nested-7 FractalForums.com Banner Logos mrob 0 1849 Last post October 22, 2010, 02:36:35 PM
by mrob
Too much nested child boards Discuss Fractal Forums PhotoComix 8 1194 Last post June 25, 2012, 06:24:17 PM
by Alef
Best way to make these tree iterations? Arrays? location marikers? nested etc Programming Mrz00m 3 3729 Last post October 01, 2012, 11:03:21 PM
by cKleinhuis
Worms Crawl Into My Brain Divergence Fragmentarium Gallery 3dickulus 0 604 Last post March 04, 2016, 02:35:53 AM
by 3dickulus

Powered by MySQL Powered by PHP Powered by SMF 1.1.21 | SMF © 2015, Simple Machines

Valid XHTML 1.0! Valid CSS! Dilber MC Theme by HarzeM
Page created in 0.153 seconds with 24 queries. (Pretty URLs adds 0.01s, 2q)