MaplePrimes - Questions, Posts, help, musings, answers ... all things Maple and math

Ask a Question

Create a Post

MaplePrimes Announcement

The Maple Conference is Back!

by: Kat 147 Products: Maple Maple Toolboxes

February 15 2019

5 0

It is my pleasure to announce the return of the Maple Conference! On October 15-17^th, in Waterloo, Ontario, Canada, we will gather a group of Maple enthusiasts, product experts, and customers, to explore and celebrate the different aspects of Maple.

Specifically, this conference will be dedicated to exploring Maple’s impact on education, new symbolic computation algorithms and techniques, and the wide range of Maple applications. Attendees will have the opportunity to learn about the latest research, share experiences, and interact with Maple developers.

In preparation for the conference we are welcoming paper and extended abstract submissions. We are looking for presentations which fall into the broad categories of “Maple in Education”, “Algorithms and Software”, and “Applications of Maple” (a more extensive list of topics can be found here).

You can learn more about the event, plus find our call-for-papers and abstracts, here: https://www.maplesoft.com/mapleconference/

Read Full Post

Featured Post

acer

15618

2D contour plot and legend

There have been several posts, over the years, related to visual cues about the values associated with particular 2D contours in a plot.

Some people ask or post about color-bars [1]. Some people ask or post about inlined labelling of the curves [1, 2, 3, 4, 5, 6, 7]. And some post about mouse popup/hover-over functionality [1]., which got added as general new 2D plot annotation functionality in Maple 2017 and is available for the plots:-contourplot command via its contourlabels option.

Another possibility consists of a legend for 2D contour plots, with distinct entries for each contour value. That is not currently available from the plots:-contourplot command as documented. That may be possible, possibly with effort, from contour-plotting approaches based on individual plots:-implicitplot calls for each contour level [eg. Kitonum's procedure, or an undocumented, alternate internal driver for plots:-contourplot].

Since I like the functionality provided by the contourlabels option I thought that I'd highjack that (and the _HOVERCONTENT plotting substructure that plot-annotations now generate) and get a relatively convenient way to get a color-key via the 2D plotting legend. This is not supposed to be super-efficient.

Here below are some examples. I hope that it illustrates some useful functionality that could be added to the contourplot command. It can also be used to get a color-key for use with densityplot.

restart;

contplot:=proc(ee, rng1, rng2)
  local clabels, clegend, i, ncrvs, newP, otherdat, others, tcrvs, tempP;
  (clegend,others):=selectremove(type,[_rest],identical(:-legend)=anything);
  (clabels,others):= selectremove(type,others,identical(:-contourlabels)=anything);
  if nops(clegend)>0 then
    tempP:=:-plots:-contourplot(ee,rng1,rng2,others[],
                                ':-contourlabels'=rhs(clegend[-1]));
    tempP:=subsindets(tempP,'specfunc(:-_HOVERCONTENT)',
                      u->`if`(has(u,"null"),NULL,':-LEGEND'(op(u))));
    if nops(clabels)>0 then
      newP:=plots:-contourplot(ee,rng1,rng2,others[],
                              ':-contourlabels'=rhs(clabels[-1]));
      tcrvs:=select(type,[op(tempP)],'specfunc(CURVES)');
      (ncrvs,otherdat):=selectremove(type,[op(newP)],'specfunc(CURVES)');
      return ':-PLOT'(seq(':-CURVES'(op(ncrvs[i]),op(indets(tcrvs[i],'specfunc(:-LEGEND)'))),
                          i=1..nops(ncrvs)),
                      op(otherdat));
    else
      return tempP;
    end if;
  elif nops(clabels)>0 then
    return plots:-contourplot(ee,rng1,rng2,others[],
                              ':-contourlabels'=rhs(clabels[-1]));
  else
    return plots:-contourplot(ee,rng1,rng2,others[]);
  end if;
end proc:

>	contplot(x^2+y^2, x=-2..2, y=-2..2, coloring=["Yellow","Blue"], contours = 9, size=[500,400], legendstyle = [location = right], legend=true, contourlabels=true, view=[-2.1..2.1,-2.1..2.1] );

>	contplot(x^2+y^2, x=-2..2, y=-2..2, coloring=["Yellow","Blue"], contours = 17, size=[500,400], legendstyle = [location = right], legend=['contourvalue',$("null",7),'contourvalue',$("null",7),'contourvalue'], contourlabels=true, view=[-2.1..2.1,-2.1..2.1] );

# Apparently legend items must be unique, to persist on document re-open.

contplot(x^2+y^2, x=-2..2, y=-2..2,
      coloring=["Yellow","Blue"],
      contours = 11,
      size=[500,400],
      legendstyle = [location = right],
      legend=['contourvalue',seq(cat($(` `,i)),i=2..5),
              'contourvalue',seq(cat($(` `,i)),i=6..9),
              'contourvalue'],
      contourlabels=true,
      view=[-2.1..2.1,-2.1..2.1]
);

>	contplot(x^2+y^2, x=-2..2, y=-2..2, coloring=["Green","Red"], contours = 8, size=[400,450], legend=true, contourlabels=true );

>	contplot(x^2+y^2, x=-2..2, y=-2..2, coloring=["Yellow","Blue"], contours = 13, legend=['contourvalue',$("null",5),'contourvalue',$("null",5),'contourvalue'], contourlabels=true );

(low,high,N):=0.1,7.6,23:
conts:=[seq(low..high*1.01, (high-low)/(N-1))]:
contplot(x^2+y^2, x=-2..2, y=-2..2,
      coloring=["Yellow","Blue"],
      contours = conts,
      legend=['contourvalue',$("null",floor((N-3)/2)),'contourvalue',$("null",ceil((N-3)/2)),'contourvalue'],
      contourlabels=true
);

plots:-display(
  subsindets(contplot((x^2+y^2)^(1/2), x=-2..2, y=-2..2,
                      coloring=["Yellow","Blue"],
                      contours = 7,
                      filledregions),
             specfunc(CURVES),u->NULL),
  contplot((x^2+y^2)^(1/2), x=-2..2, y=-2..2,
      coloring=["Yellow","Blue"],
      contours = 7, #grid=[50,50],
      thickness=0,
      legendstyle = [location=right],
      legend=true),
  size=[600,500],
  view=[-2.1..2.1,-2.1..2.1]
);

plots:-display(
  contplot(x^2+y^2, x=-2..2, y=-2..2,
      coloring=["Yellow","Blue"],
      contours = 5,
      thickness=0, filledregions),
  contplot(x^2+y^2, x=-2..2, y=-2..2,
      coloring=["Yellow","Blue"],
      contours = 5,
      thickness=3,
      legendstyle = [location=right],
      legend=typeset("<=",contourvalue)),
  size=[700,600],
  view=[-2.1..2.1,-2.1..2.1]
);

N:=11:
plots:-display(
  contplot(sin(x)*y, x=-2*Pi..2*Pi, y=-1..1,
      coloring=["Yellow","Blue"],
      contours = [seq(-1+(i-1)*(1-(-1))/(N-1),i=1..N)],
      thickness=3,
      legendstyle = [location=right],
      legend=true),
   plots:-densityplot(sin(x)*y, x=-2*Pi..2*Pi, y=-1..1,
      colorscheme=["zgradient",["Yellow","Blue"],colorspace="RGB"],
      grid=[100,100],
      style=surface, restricttoranges),
   plottools:-line([-2*Pi,-1],[-2*Pi,1],thickness=3,color=white),
   plottools:-line([2*Pi,-1],[2*Pi,1],thickness=3,color=white),
   plottools:-line([-2*Pi,1],[2*Pi,1],thickness=3,color=white),
   plottools:-line([-2*Pi,-1],[2*Pi,-1],thickness=3,color=white),
   size=[600,500]
);

N:=13:
plots:-display(
  contplot(sin(x)*y, x=-2*Pi..2*Pi, y=-1..1,
      coloring=["Yellow","Blue"],
      contours = [seq(-1+(i-1)*(1-(-1))/(N-1),i=1..N)],
      thickness=6,
      legendstyle = [location=right],
      legend=['contourvalue',seq(cat($(` `,i)),i=2..3),
              'contourvalue',seq(cat($(` `,i)),i=5..6),
              'contourvalue',seq(cat($(` `,i)),i=8..9),
              'contourvalue',seq(cat($(` `,i)),i=11..12),
              'contourvalue']),
   plots:-densityplot(sin(x)*y, x=-2*Pi..2*Pi, y=-1..1,
      colorscheme=["zgradient",["Yellow","Blue"],colorspace="RGB"],
      grid=[100,100],
      style=surface, restricttoranges),
   plottools:-line([-2*Pi,-1],[-2*Pi,1],thickness=6,color=white),
   plottools:-line([2*Pi,-1],[2*Pi,1],thickness=6,color=white),
   plottools:-line([-2*Pi,1],[2*Pi,1],thickness=6,color=white),
   plottools:-line([-2*Pi,-1],[2*Pi,-1],thickness=6,color=white),
  size=[600,500]
);

Download contour_legend_post.mw

February 12

9 1

Read Full Post

Featured Post

ecterrab

6543

A Complete Guide for...

A Complete Guide for performing Tensors computations using Physics

This is an old request, a complete guide for using Physics to perform tensor computations. This guide, shown below with Sections closed, is linked at the end of this post as a pdf file with all the sections open, and also as a Maple worksheet that allows for reproducing its contents. Most of the computations shown are reproducible in Maple 2018.2.1, and a significant part also in previous releases, but to reproduce everything you need to have the Maplesoft Physics Updates version 283 or higher installed. Feedback one how to improve this presentation is welcome.

Physics is a package developed by Maplesoft, an integral part of the Maple system. In addition to its commands for Quantum Mechanics, Classical Field Theory and General Relativity, Physics includes 5 other subpackages, three of them also related to General Relativity: Tetrads , ThreePlusOne and NumericalRelativity (work in progress), plus one to compute with Vectors and another related to the Standard Model (this one too work in progress).

The presentation is organized as follows. Section I is complete regarding the functionality provided with the Physics package for computing with tensors in Classical and Quantum Mechanics (so including Euclidean spaces), Electrodynamics and Special Relativity. The material of section I is also relevant in General Relativity, for which section II is all devoted to curved spacetimes. (The sub-section on the Newman-Penrose formalism needs to be filled with more material and a new section devoted to the EnergyMomentum tensor is appropriate. I will complete these two things as time permits.) Section III is about transformations of coordinates, relevant in general.

For an alphabetical list of the Physics commands with a brief one-line description and a link to the corresponding help page see Physics: Brief description of each command .

I. Spacetime and tensors in Physics

This section contains all what is necessary for working with tensors in Classical and Quantum Mechanics, Electrodynamics and Special Relativity. This material is also relevant for computing with tensors in General Relativity, for which there is a dedicated Section II. Curved spacetimes .

	Default metric and signature, coordinate systems

Tensors, their definition, symmetries and operations

Physics comes with a set of predefined tensors, mainly the spacetime metric , the space metric , and all the standard tensors of General Relativity. In addition, one of the strengths of Physics is that you can define tensors, in natural ways, by indicating a matrix or array with its components, or indicating any generic tensorial expression involving other tensors.

In Maple, tensor indices are letters, as when computing with paper and pencil, lowercase or upper case, latin or greek, entered using indexation, as in , and are displayed as subscripts as in . Contravariant indices are entered preceding the letter with ~, as in , and are displayed as superscripts as in . You can work with two or more kinds of indices at the same time, e.g., spacetime and space indices.

To input greek letters, you can spell them, as in mu for , or simpler: use the shortcuts for entering Greek characters . Right-click your input and choose Convert To → 2-D Math input to give to your input spelled tensorial expression a textbook high quality typesetting.

Not every indexed object or function is, however, automatically a tensor. You first need to define it as such using the Define command. You can do that in two ways:

1.	Passing the tensor being defined, say , possibly indicating symmetries and/or antisymmetries for its indices.

2.	Passing a tensorial equation where the left-hand side is the tensor being defined as in 1. and the right-hand side is a tensorial expression - or an Array or Matrix - such that the components of the tensor being defined are equal to the components of the tensorial expression.

After defining a tensor - say or - you can perform the following operations on algebraic expressions involving them

•	Automatic formatting of repeated indices, one covariant the other contravariant

•	Automatic handling of collisions of repeated indices in products of tensors

•	Simplify products using Einstein's sum rule for repeated indices.

•	SumOverRepeatedIndices of the tensorial expression.

•	Use TensorArray to compute the expression's components

•	TransformCoordinates .

If you define a tensor using a tensorial equation, in addition to the items above you can:

•	Get each tensor component by indexing, say as in or

•	Get all the covariant and contravariant components by respectively using the shortcut notation and .

•	Use any of the special indexing keywords valid for the pre-defined tensors of Physics; they are: definition, nonzero, and in the case of tensors of 2 indices also trace, and determinant.

•	No need to specify the tensor dependency for differentiation purposes - it is inferred automatically from its definition.

•	Redefine any particular tensor component using Library:-RedefineTensorComponent

•	Minimizing the number of independent tensor components using Library:-MinimizeTensorComponent

•	Compute the number of independent tensor components - relevant for tensors with several indices and different symmetries - using Library:-NumberOfTensorComponents .

The first two sections illustrate these two ways of defining a tensor and the features described. The next sections present the existing functionality of the Physics package to compute with tensors.

	Defining a tensor passing the tensor itself

	Defining a tensor passing a tensorial equation

	Automatic formatting of repeated tensor indices and handling of their collisions in products

	Tensor symmetries

	Substituting tensors and tensor indices

	Simplifying tensorial expressions

	SumOverRepeatedIndices

	Visualizing tensor components - Library:-TensorComponents and TensorArray

	Modifying tensor components - Library:-RedefineTensorComponent

	Enhancing the display of tensorial expressions involving tensor functions and derivatives using CompactDisplay

	The LeviCivita tensor and KroneckerDelta

	The 3D space metric and decomposing 4D tensors into their 3D space part and the rest

	Total differentials, the d_[mu] and dAlembertian operators

	Tensorial differential operators in algebraic expressions

	Inert tensors

	Functional differentiation of tensorial expressions with respect to tensor functions

	The Pauli matrices and the spacetime Psigma[mu] 4-vector

	The Dirac matrices and the spacetime Dgamma[mu] 4-vector

	Quantum not-commutative operators using tensor notation

II. Curved spacetimes

Physics comes with a set of predefined tensors, mainly the spacetime metric , the space metric , and all the standard tensors of general relativity, respectively entered and displayed as: Einstein[mu,nu] = , Ricci[mu,nu] = , Riemann[alpha, beta, mu, nu] = , Weyl[alpha, beta, mu, nu], = , and the Christoffel symbols Christoffel[alpha, mu, nu] = and Christoffel[~alpha, mu, nu] = respectively of first and second kinds. The Tetrads and ThreePlusOne subpackages have other predefined related tensors. This section is thus all about computing with tensors in General Relativity.

	Loading metrics from the database of solutions to Einstein's equations

	Setting the spacetime metric indicating the line element or a Matrix

	Covariant differentiation: the D_[mu] operator and the Christoffel symbols

	The Einstein, Ricci, Riemann and Weyl tensors of General Relativity

	A conversion network for the tensors of General Relativity

	Tetrads and the local system of references - the Newman-Penrose formalism

	The ThreePlusOne package and the 3+1 splitting of Einstein's equations

	III. Transformations of coordinates

	See Also
	Physics , Conventions used in the Physics package , Physics examples , Physics Updates

Download Tensors_-_A_Complete_Guide.mw, or the pdf version with sections open: Tensors_-_A_Complete_Guide.pdf

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

January 23

8 1

Read Full Post

How MapleSim Helped Canada’s Paralympic Curling...

by: John McPhee 20 Products: Maple MapleSim

January 23 2019

4 0

MaplePrimes Questions

Recent

Unanswered

Maple

MapleSim

Share via:

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

MaplePrimes Announcement

The Maple Conference is Back!

Featured Post

2D contour plot and legend

Featured Post

A Complete Guide for...

How MapleSim Helped Canada’s Paralympic Curling...

MaplePrimes Questions

Recent

Unanswered

Maple

MapleSim

What not to enter in 2d math ?

simplifying a complex symbolic expression

What's wrong with the spherical version of Bessel...

Generating random number from a regression model

Closest fit position of line

convert Surface data to region

Problem when iterating a discrete map

Procedure for Sort "by proxy" (sort list using...

seq for sparse rtables

How do I create a vector / array of time dependent...

Put constraint on symbolic calculation

Getting an error when using pointplot command

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again

Ask a Question

Create a Post

MaplePrimes Announcement

Featured Post

Featured Post

MaplePrimes Questions Recent Unanswered Maple MapleSim

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

MaplePrimes Questions

Recent

Unanswered

Maple

MapleSim