| Bernstein Tensor product polynomials | 
These polynomials are represented as an array of coefficients in the tensor product Bernstein basis.
The corresponding type is Polynomial Bernstein C, where C is the type of coefficients. The coefficients C should be a field. Here we describe the main functionnalities available for these polynomials.
               
                    Mmx]  
                   
                     
                      use "realroot"
                    
                
            
                   
              
               
                    Mmx]  
                   
                     
                      R := QQ['x,'y,'z ]
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      B := bernstein R
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      b := B << "x^2+x+1"
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      type(b:> Generic)
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      coefficients b
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      b2 := b*b
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      coefficients(b2)
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      diff(b2,0)
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      b3 := B<<"x*y^2+x-5"
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      coefficients b3
                    
                
            
                   
              
               
                       
                    
               
                    Mmx]  
                   
                     
                      coefficients (b3, 0)
                    
                
            
                   
              
               
                    Mmx]